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block  20  Continued 

Phase  I  contained  both  theoretical  development  and  experimental  tests 
designed  to  answer  questions  in  the  application  of  thermocouples  for  brake 
temperature  measurement.  A  theoretical  model  was  developed  containing  conduc¬ 
tion,  convection,  and  radiation  heat  transfer  to  a  thermocouple  bead  placed  in 
a  hole  at  the  bottom  of  which  the  temperature  was  desired.  The  effects  of 
thermocouple  lead  wire  diameter,  bead  diameter  and  distance  from  the  hole  bot¬ 
tom  were  explored.  Experimentally,  nineteen  tests  were  conducted  on  ten  type 
K  thermocouple  configurations  in  a  bench  test  apparatus  consisting  of  a  block 
of  carbon  approximately  .45  x  2.  x  3.  in.  The  carbon  block  was  heated  from 
each  side  by  a  bank  of  high  intensity  tungsten-filament  quartz  lamps.  Each 
test  consisted  of  heating  the  thermocouple  to  1500°F  in  approximately  50 
seconds.  The  best  two  thermocouple  configurations  were  the  sheathed  type,  one 
of  which  had  a  grounded  junction,  and  the  other  of  which  had  an  exposed  junc¬ 
tion. 

Phase  II,  the  brake  theoretical  analysis,  contained  the  development  of  a 
transient  finite  difference  theoretical  model  to  predict  the  temperature  within 
the  stators,  rotors,  and  brake  housing  of  a  typical  disc  brake.  The  model 
contains  provisions  for  predicting  temperature  wiLhin  the  wheel,  rim,  and  tire. 
The  effort  included  development  of  one  dimensional,  two  dimensional  and  three 
I  dimensional  heat  transfer  models.  The  effort  is  unique  in  that  it  combines 
variable  thermal  conductivity,  variable  heat  capacity,  rotor-stator  interface 
contact  resistance,  experimentally  determined  radiative/convective  heat  trans¬ 
fer  coefficients,  and  allowable  temperature  variation  from  stator  to  stator 
and  rotor  to  rotor.  The  inclusion  of  variable  thermal  conductivity  is  especial¬ 
ly  important  with  the  use  of  carbon  composite  brakes  because  the  conductivity 
is  variable  not  only  with  temperature,  but  also  with  direction.  This  added 
variable  is  due  to  the  carbon  disc  build-up  by  ply  direction  or  random  fiber 
orientation.  In  addition  to  the  multidisc  theoretical  development,  a  caliper 
type  brake  is  analyzed,  and  a  finite  difference  model  including  all  required 
equations  is  developed. 

Phase  III,  the  full  scale  brake  tests,  employed  an  F-14  size,  four  rotor, 
multidisc  brake  tested  with  a  192  inch  diameter  inertial  dynamometer.  Two 
heat  sink  materials  were  tested,  beryllium  and  carbon  composite,  so  that 
different  thermal  conductivities,  densities ,  and  heat  capacities  could  be  used. 
Tests  were  conducted  at  approximately  6.7  and  11.2  million  ft-lbs  of  energy 
with  decelerations  ranging  from  about  2.  to  8.  ft/sec^  and  initial  velocities 
ranging  from  207.  to  274.  ft/sec.  Test  results  were  used  to  calculate  convec¬ 
tive/radiative  heat  transfer  coefficients,  and  for  comparison  to  theory. 

The  correlation  of  results.  Phase  IV,  showed  that  the  one  dimensional 
axial,  and  the  one  dimensional  radial  models  predict  peak  temperatures  during  a 
brake  stop  that  are  up  to  31%  and  16%  higher  than  experiment,  respectively. 

The  two  dimensional  model  predicts  peak  temperatures  within  5%  of  experiment  in 
three  out  of  four  cases  considered,  with  the  fourth  case  less  than  8%.  The 
peak  temperature  predictions  within  the  torque  tube,  fuse  plug,  drive  lug, 
bead  seat  and  contained  air  were  within  10%  of  experiment  for  most  cases. 
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CHAPTER  I 


k 


INTRODUCTION 

Background 

The  primary  function  of  a  brake  system  is  to  convert  the  kinetic 
energy  of  vehicle  motion  into  internal  energy  while  the  vehicle  is 
being  stopped.  In  addition  to  this  primary  function,  the  brake  must 
store  the  absorbed  heat  for  a  finite  amount  of  time  so  that  overheating 
of  surrounding  vehicle  hardware  will  not  occur.  Two  methods  for 
accomplishing  these  functions  are  the  use  of  drum  brakes  and  disc 
brakes.  Disc  brakes  are  used  in  aircraft,  trucks,  and  many  automobile 
braking  systems. 

A  typical  disc  brake  assembly  as  used  in  numerous  car,  truck,  and 
aircraft  applications  consists  of  two  basic  parts,  the  brake  heat  sink 
and  the  brake  housing  (Figure  1) .  The  heat  sink  is  built  up  of  alter¬ 
nating  stators  and  rotors  which  take  the  shape  of  annuli.  A  photograph 
of  actual  hardware  from  an  aircraft  is  shown  in  Figure  2.  The  rotors 
are  keyed  to  the  rotating  wheel  while  the  stators  are  keyed  to  the 
fixed  vehicle  axle.  The  number  of  keys  is  dependent  upon  the  particular 
application.  During  vehicle  stopping,  the  rotors  and  stators  are 
pushed  together  using  hydraulically  actuated  pistons  such  that  the 
vehicle  kinetic  energy  is  transformed  into  heat  due  to  the  friction 

between  the  stators  and  rotors.  The  contact  pressure  is  normally 
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maintained  at  a  relatively  constant  value  as  the  vehicle  is  brought  to 
a  stop.  This  results  in  a  time-varying  energy  input  into  the  stators 
and  rotors.  During  this  energy  input,  heat  is  lost  due  to  convection 
and  radiation  from  the  edges  of  the  stators  and  rotors  and  by  conduc¬ 
tion  at  the  key  locations  and  end  stators. 

Statement  of  the  Problem 

The  disc  brake  arrangement  described  above  must  be  designed  to 
meet  the  two  functions  previously  mentioned.  In  order  to  do  this,  the 
designer  must  have  brake  temperature  data  available  to  him.  Specifi¬ 
cally,  the  temperature  distribution  in  the  heat  sink  is  required  to 
determine  the  structural  heating  effect  on  the  stators,  rotors,  brake 
housing,  wheel,  and  tire.  The  tire  peak  temperature  is  required  to 
prevent  deterioration  of  the  rubber  properties  and  subsequent  failure 
due  to  strength  loss.  The  stator,  rotor,  brake  housing  and  wheel 
temperatures  and  gradients  are  required  so  that  thermally  induced 
stress  patterns  and  fatigue  predictions  can  be  made.  The  cooling  rate 
is  also  important  in  determining  the  vehicle  permissible  "time  between 
stops."  Since  brake  wear  in  common  heat  sinks  is  related  to  heat  sink 
rubbing  area  and  since  high  temperatures  cause  loss  of  material  through 
oxidation,  actual  heat  sink  temperatures  become  a  function  of  brake 
life.  In  order  that  the  brake  lubrication  seals  are  protected,  it  is 
necessary  that  the  hydraulic  fluid,  which  is  used  to  pressurize  the 
pistons,  remain  below  a  given  temperature.  Thus  it  is  important  to  be 
able  to  predict  the  hydraulic  fluid  temperature. 

The  objective  of  this  research  is  to  develop  a  general  model  for 
accurately  predicting  temperatures  within  the  stators,  rotors,  and 
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brake  housing  of  a  typical  disc  brake.  The  model  will  contain  pro¬ 
visions  for  the  calculation  of  temperature  within  the  wheel,  axle,  and 
tire.  One-,  two-,  and  three-dimensional  models  are  needed  for  com¬ 
parison  with  each  other,  as  well  as  with  experiment.  More  than  one  theore¬ 
tical  model  is  required  so  that  the  most  accurate  and  most  expedient 
model  for  a  particular  application  can  be  identified.  This  research 
is  to  apply  equally  well  to  caliper  type  brakes  (single-rotor)  and  to 
multidisc  type  brakes.  Separate  theoretical  models  are  required  to 
accomplish  this  generality.  Experiments  for  comparison  to  the  theory 
are  an  extremely  important  part  of  this  effort  since  the  validity  of  the 
theoretical  models  will  ultimately  be  judged  by  how  close  the  experimental 
curves  are  matched. 
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CHAPTER  II 


REVIEW  OF  THE  LITERATURE 

Overall 

Efforts  at  obtaining  the  temperature  distribution  in  a  brake  and 
adjacent  hardware  can  be  divided  into  two  classes  depending  upon  the 
particular  style  brake  analyzed.  The  first  is  a  single-rotor  type 
brake  in  which  the  stators  are  replaced  by  a  pair  of  sector-shaped 
calipers.  Energy  is  then  transferred  from  the  rotor  into  the  calipers. 
The  second  differs  from  the  first  in  that  the  stator  is  always  a 
complete  annulus.  While  class  one  is  limited  to  a  single  rotor,  class 
two  can  have  any  number  of  rotors  and  corresponding  stators.  Clearly, 
class  one  is  a  special  case  of  class  two. 

Class  One — Caliper 

•k 

The  class  one  type  brake  has  been  analyzed  by  Evans  and  Newcomb.* 
They  assume  that  temperature  is  not  a  function  of  the  radial  or  angular 
directions  within  the  brake.  Newcomb  and  El-Sherbiny2  used  a  similar 
analytical  model  but  added  in  a  liquid  cooling  boundary  condition  and 
assumed  a  constant  value  for  thermal  conductivity.  Secrist  and 
Hornbeck3  assumed  that  temperature  varied  in  the  radial  direction  but 
remained  constant  in  the  axial  and  angular  directions.  They  employed 
a  finite  difference  technique  to  obtain  a  solution.  An  experimental 
as  well  as  analytical  analysis  was  conducted  by  Limpert1*  in  which 
temperature  was  assumed  to  vary  in  the  axial  direction  only.  He 
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assumed  a  constant  thermal  conductivity  and  used  a  finite  difference 
technique.  Hartter,  Schwartz,  and  Rhee5  conducted  tests  on  a  heavy 
truck  brake  and  used  the  results  to  refine  their  finite  difference 
theoretical  analysis.  One-dimensional  conduction  in  the  axial  direc¬ 
tion  was  modeled  assuming  constant  thermal  conductivity,  Limpert8 
again  assumed  that  temperature  varied  in  the  axial  direction  only,  but 
included  within  his  finite  difference  technique  a  method  for  calcu¬ 
lating  the  stress  in  the  rotor.  Rusnak,  Schwartz,  and  Coleman7 
numerically  analyzed  the  axial  direction  allowing  for  variable  thermal 
conductivity  with  the  intent  of  rating  several  different  materials.  A 
study  was  conducted  by  Dike8  that  solved  the  heat  conduction  equation 
in  the  radial  and  axial  directions,  but  this  study  assumed  that  the 
heat  into  the  brake  at  the  sliding  surface  was  independent  of  radius. 
Santini  and  Kennedy8  conducted  an  entirely  experimental  effort  in 
which  the  caliper  part  of  the  brake  was  imbedded  with  thermocouples. 
They  pointed  out  the  lack  of  published  brake  experimental  data.  These 
two  authors,  along  with  Ling,10  conducted  constant  speed  experiments 
with  the  intent  of  evaluating  external  cooling  and  rotor  design  modifi¬ 
cations.  Schwartz11  used  the  finite  difference  technique  for  looking 
at  the  radial  and  axial  directions.  The  caliper  was  not  analyzed  and 
the  convection  heat  transfer  coefficients  were  entirely  theoretical. 
Sinclair  and  Gulick12  performed  no  analyses  but  discussed  the  advan¬ 
tages  of  using  a  dynamometer  for  testing  rather  than  a  vehicle. 

Morgan  and  Dennis13  used  the  finite  element  technique  in  the  radial 
and  axial  directions  to  predict  temperature.  Parker  and  Newcomb114 


developed  an  exact  procedure  for  predicting  the  surface  temperature  of 
a  caliper  rotor  as  a  function  of  time. 

Class  Two — Multidisc 

The  class  two  type  of  brake  (annular  stators)  has  been  studied  by 
Anderle15  with  the  assumption  of  axial  heat  transfer  only.  He 
analyzed  a  rotor-stator  pair  and  assumed  that  all  other  pairs  would 
have  the  same  temperature  profile.  A  finite  difference  approach  was 
used  in  which  half  of  the  input  energy  at  the  sliding  surface  was 
assumed  to  go  into  the  stator  with  the  remaining  half  going  into  the 
rotor.  Claridge18  developed  a  brake  model  that  considered  the  axial 
direction  only.  His  main  objective  was  to  show  how  a  computer  could 
be  used  in  the  design  of  the  complete  braking  system.  Scaringe,  Ho, 
and  Peterson17  conducted  a  primarily  experimental  effort  to  evaluate 
an  air-cooling  scheme.  The  theoretical  model  alluded  to  assumed  no 
heat  conduction  in  the  radial  or  angular  directions.  Assumptions  of 
constant  thermal  conductivity  and  constant  heat  input  were  also  made. 
Kennedy,  Wu,  and  Ling18  analyzed  a  stator-rotor  pair  employing  the 
finite  element  approach  with  assumptions  of  heat  transfer  in  the  axial 
direction  only.  The  results  were  compared  to  experiments.  Mercer18 
has  used  a  finite  difference  approach  to  analyze  a  rotor-stator  pair, 
allowing  heat  transfer  in  the  radial  as  well  as  axial  directions.  The 
energy  input  is  varied  with  radius  as  well  as  time.  The  angular 
direction  is  ignored,  as  well  as  any  losses  to  surrounding  components. 
Experimental  studies  have  been  accomplished  by  Ho  and  Peterson28, 21 
in  which  sliding  tests  were  used  to  simulate  braking.  These  tests 
were  aimed  at  materials  development.  Unfortunately,  it  has  been 


mi 


found,  sliding-type  tests  do  not  adequately  simulate  the  conditions 
within  the  brake  due  to  different  geometry  and  cooling  conditions. 

Ho,  Kennedy,  and  Peterson22  conducted  two-rotor  tests  in  which  the 
temperature  at  locations  within  the  stators  was  monitored.  Wagner23 
reports  an  experimental  comparison  of  two  materials  run  within  a  two- 
rotor  brake.  Results  consistently  show  both  radial  and  axial  (stator- 
to-stator  and  rotor-to-rotor)  temperature  variations  that  are  signifi¬ 
cant  (up  to  40%) .  Bender  and  Han24  used  a  finite  difference  model  to 
look  at  brake  temperatures  but  lumped  each  stator  and  rotor  so  that  no 
temperature  variation  was  allowed  within  the  stator  or  rotor.  Brewer25 
analyzed  one  rotor  using  a  one-dimensional  finite  difference  approach 
with  constant  properties.  Donaldson,  Stocks,  Teter,  and  Rebholz26 
also  used  finite  differences  but  looked  at  the  radial  as  well  as  the 
axial  direction.  They  analyzed  a  single  rotor  and  stator  and  assumed 
that  the  others  were  identical.  They  also  found  a  theoretical  con¬ 
vective  heat  transfer  coefficient.  Lee27  analyzed  all  three  directions 
but  performed  no  experiments,  assumed  an  energy  split  between  the 
rotor  and  stator,  used  theoretical  convection  heat  transfer  coefficients, 
and  looked  at  only  one  pair  of  discs.  Nemeth  and  Wieter2®  employed  a 
finite  difference  approach  but  used  a  lumped  system  that  resulted  in 
only  one  temperature  for  the  rotor  and  one  for  the  stator.  Stout29 
conducted  tests  on  a  brake  system  designed  for  a  space  application  in 
which  the  ambient  temperature  was  800°F. 


Summary 


In  summary,  past  analyses  have  assumed  that  all  rotors  behaved 


the  same  and  that  all  stators  behaved  the  same.  This  ignores  heat 
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loss  through  the  end  stators.  No  theoretical  model  is  currently 
available  that  accounts  for  the  experimentally  measured  axial  tempera¬ 
ture  variation  from  disc  to  disc.  With  such  an  analytical  model,  a 
brake  heat  sink  could  be  designed  such  that  it  would  have  a  minimum 
axial  temperature  gradient  at  the  most  used  stopping  conditions.  This 
would  reduce  the  high  wear  oxidation  rate  at  the  present  heat  sink 
high  temperature  area.  Past  efforts  have  assumed  constant  thermal 
properties  which  could  result  in  significant  error,  especially  at  the 
high  temperatures  that  are  being  routinely  encountered  in  present-day 
heat  sink  materials  (1500°F) .  Often  an  estimate  is  made  of  the 
fraction  of  generated  heat  that  goes  into  a  rotor  and  stator.  This 
estimate  is  not  required.  Also,  it  is  apparent  that  no  comparison  of 
model  complexities  has  been  reported.  Thus,  it  is  unclear  which  model 
should  be  used  in  a  particular  application.  Extensive  testing  is  also 
sketchy  in  the  literature.  When  tests  are  conducted,  they  are  usually 
not  accompanied  by  a  theoretical  development.  Table  1  is  a  summary 
table  that  lists  authors,  brake  type,  whether  effort  was  experimental 
or  theoretical,  and  whether  thermal  properties  were  varied  with 
temperature. 


TABLE  1 — continued 
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CHAPTER  III 


TEMPERATURE  MEASUREMENT  ANALYSIS 

Introduction 

In  any  experimental  program  it  is  important  to  know  the  limitations 
of  the  equipment  being  used.  In  this  research  the  measurement  of  tem¬ 
perature  of  brakes  and  their  surrounding  hardware  is  a  major  considera¬ 
tion.  Different  types  and  sizes  of  thermocouples  are  routinely  used 
interchangeably  by  many  organizations  with  little  thought  given  to  the 
resulting  accuracy  or  response  time  of  the  thermocouple  being  used. 
Thermocouple  manufacturers  often  give  time  constant  data,  but  this  is 
only  valid  for  a  particular  application.  The  problems  inherent  in 
thermocouple  temperature  measurement  are  discussed  in  the  literature 
along  with  recommended  methods  to  construct,  locate,  and  use  thermo¬ 
couples.  Although  this  information  is  important,  the  deviations  that 
could  result  from  compromises  that  may  have  to  be  made  as  a  result  of 
geometry,  availability,  time,  resources,  or  durability  are  not  dis¬ 
cussed  in  the  literature.  Clearly,  a  need  exists  in  a  particular 
application  for  answers  to  some  of  these  thermocouple  unknowns  in 
order  to  obtain  the  best  results. 

Theoretical  Model 

A  model  is  used  to  determine  the  temperature  reached  within  a 
thermocouple  bead  (the  welded  junction  of  the  two  dissimilar  metals  of 
the  thermocouple)  as  a  function  of  the  temperature  of  the  material 


surrounding  the  bead.  In  order  to  accomplish  this  analysis,  a  typical 
temperature  distribution  in  a  disc  of  a  brake  must  be  found,  as  well  as 
the  temperature  distribution  in  the  air  within  the  hole  into  which  the 
thermocouple  is  embedded.  This  model  is  shown  in  Figure  3. 
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Figure  3.  Temperature  Measurement  Analysis  Theoretical  Model 
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One-Dimensional  Brake  Material  Analysis 
first  step,  the  one-dimensional  temperature  distribution  is 
Figure  4  shows  the  model  used  for  obtaining  this  temperature 


H - ^ - H 


4.  One-Dimensional  Axial  Temperature  Distribution  Model 


distribution  in  the  brake  material.  The  governing  one-dimensional 
transient  equation30  in  dimensionless  form  is: 


^  «... 


where  T'  is  temperature,  y'  is  distance  and  t'  is  time.  [See 
Equation  (4)  for  definition.]  The  initial  condition  is: 

T'(y\  o)  =  0  (2) 

One  boundary  condition  is : 

t'  TP-T 

T'(o,  t')  =  — -y-^-  (3) 

sn  o 


where  t'  is  the  dimensionless  stop  time  (time  required  for  the  brake 
to  come  to  rest)  ,  Tq  is  the  initial  temperature  and  TP  is  the  peak 
temperature  during  the  stop.  Equation  (3)  was  selected  as  the  boundary 
condition  at  y'  =0  since  experiment  showed  that  a  linear  surface 
temperature  increasing  to  a  peak  at  the  stop  time  was  a  good  approxi¬ 
mation  to  the  actual  temperature  rise.  In  terms  of  dimensional  vari¬ 
ables  (depicted  without  primes),  the  above  are: 


where  a  is  the  thermal  diffusivity  and  Tq  is  the  initial  temperature. 
The  second  boundary  condition  is: 


3T* 

3y’ 


y' 


=  l 


T’(l,t’) 


(5) 


where  h^/2  is  the  disc  thickness  as  seen  in  Figure  4,  K  is  the  thermal 
conductivity  and  h  is  the  aver..>;e  convective  heat  transfer  coefficient. 
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With  the  use  of  standard  finite  difference  techniques  and  the  Euler 
time  marching  technique.  Equation  (1)  becomes: 

Tt+“(n  -  t'o)  x  [i  -  (W2-°fe]  +  [Tt<1  +  « 

+  Tt(l  - »]  *  Td-t-x-w)  «> 

where  t  indicates  time,  t  +  At  is  one  time  increment  advanced  from  t,  I 
is  the  Ith  node  in  the  finite  difference  network  shown  in  Figure  4,  DY 
is  the  distance  between  two  nodes,  and  DT  is  the  time  increment.  The 
primes  have  been  omitted  from  (6)  for  clarity.  The  computer  program 
solving  Equations  (2)  through  (6)  is  shown  in  Appendix  A.  The  material 
constants  used  in  this  program  were  for  a  carbon  composite  brake 
material.  Carbon  was  selected  at  this  point  since  full-scale  tests 
(discussed  in  Chapter  4),  employed  carbon  and  beryllium,  and  the 
carbon  material  provided  the  highest  temperature  gradients.  The  time- 
varying  boundary  condition  was  selected  because  it  was  indicative  of 
the  rubbing  surface  temperature  rise  of  past  full-scale  test  results. 
Specific  numbers  were  chosen  from  a  typical  carbon  brake  test.  Results 
showed  that  for  a  peak  temperature  of  1170°F  and  an  input  rate  of 
40°F/sec,  the  actual  temperature  at  each  node  was  lower  at  each  time 
step  by  an  amount  independent  of  the  time  step.  This  is  shown  in 
Figure  5  as  a  plot  of  temperature  lag  from  the  heat  input  surface 
versus  the  dimensionless  distance. 
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Two-Dimensional  Air  Cavity  Analysis 
The  temperature  distribution  within  the  hole  shown  in  Figure  3  is 
needed  to  evaluate  the  heat  transfer  to  the  thermocouple  bead  due  to 
convection.  The  boundary  temperatures  for  this  problem  were  obtained 
from  the  previous  one-dimensional  analysis.  Figure  6  shows  the  model 
to  be  analyzed  including  the  finite  difference  grid  points  and  the  grid 
numbering  system.  The  two-dimensional  transient  equation  in  dimension¬ 
less  form  is: 


32T' 

3x’2 


+  A2 


32T' 
3y'  2 


3T 1 
3t’ 


(7) 


where  x'  and  y'  are  coordinates  shown  in  Figure  6,  T'  is  temperature, 
t'  is  time,  and  A  is  the  ratio  of  the  length  to  height  of  the  model,  as 
shown  in  Figure  6.  In  terms  of  dimensional  variables  (depicted  without 
primes),  the  above  are: 


where  a  is  thermal  diffusivity  and  Tq  is  the  initial  temperature.  The 
initial  condition  is: 


T' (x' ,  y' ,  o)  =  0  (9) 

The  boundary  condition  of  all  nodes  on  the  right  face  is  assumed: 

T'(l,  y\  t’)  =  0  (10) 

20 
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The  boundary  condition  of  all  nodes  on  the  left  face  was  again  selected 
to  be: 

,  TP-T 

T'(0,  y',  f)  =  ^ - ^  (11) 

sn  o 

where  tg  is  the  dimensionless  stop  time  and  TP  is  the  peak  temperature 
during  the  stop.  The  remaining  three  points  on  the  upper  boundary  and 
the  remaining  three  points  on  the  lower  boundary  are  assigned  tempera¬ 
tures  based  on  values  taken  from  Figure  5,  i.e.: 


T'(2, 

5, 

t')  =  T'(2, 

1, 

t') 

=  T’(0, 

y’ , 

t') 

-27.5/T 

o 

(12) 

T'(3, 

5, 

f)  =  T’(3, 

1, 

t’) 

=  T'(0, 

y\ 

f) 

-46.1/T 

o 

(13) 

T’(4, 

5, 

rt 

V-/ 

II 

H 

1, 

t') 

=  T'(0, 

y' . 

t') 

-49.9/T 

o 

(14) 

With  the  use  of  standard  finite  difference  techniques  and  the  Euler 
time  marching  technique,  Equation  (7)  becomes: 

TC  +  At(I,  J)  =  T(I ,  J)  +  [T(I-1,  J)  -2T(I,  J) 

DX2 

+  T(I+1,  J)]  +  (DT-)-  — 

DYZ 

[T(I,  J-l)  -  2T(I,  J)  +  T(I,  J  +  1)]  (15) 

where  t  indicates  time,  t  +  At  is  one  time  interval  advanced  from  t,  I 
is  the  Ith  node  in  the  x  direction,  and  J  is  the  Jth  node  in  the  y 
direction  (Figure  6).  For  clarity,  the  primes  in  (15)  have  been 
omitted.  The  computer  program  for  solving  Equations  (9)  through  (15) 
is  given  in  Appendix  B.  The  material  of  concern  in  this  part  of  the 


analysis  was  air.  It  was  assumed  that  conduction  was  the  only  heat 
transfer  mechanism,  i.e.,  convection  and  radiation  were  neglected.  The 
results  showed  for  y'  =  .5,  that  the  first  two  nodes  to  the  right  of 
the  left-hand  boundary  differed  from  that  boundary  by  a  constant  amount 
which  was  independent  of  time.  Figure  7  gives  the  temperature  (°F)  at 
y'  =  .5  as  a  function  of  x'  and  t'.  The  x'  =  0  value  is  the  bottom  of 
the  hole  or  the  carbon  temperature.  From  Figure  7  it  can  be  seen  that 
at  x'  =  .25  the  temperature  is  27°F  less  than  the  carbon  temper;. ture. 
This  information  is  used  in  the  following: 


Thermocouple  Bead  Temperature  Analysis 
The  final  step  in  analyzing  the  thermocouple  embedded  in  a  hole 
(Figure  3)  was  to  use  the  previously  developed  carbon  temperature  and 
air  cavity  temperature  as  boundary  conditions  for  the  bead 
inserted  in  the  hole.  The  first  law  of  thermodynamics  can  be  applied 
to  the  bead  of  the  inserted  thermocouple  shown  in  Figure  3.  Since  no 
work  crosses  the  bead  boundary  and  no  kinetic  energy  or  potential 
energy  changes  are  involved,  the  only  energy  change  mechanisms  are  heat 
transfer  and  internal  energy.  The  first  law  then  takes  the  form: 


The  conduction  term  is  the  rate  of  heat  conducted  away  from  the  bead 
through  the  thermocouple  lead  wires  due  to  the  temperature  gradient 
along  the  wires.  For  this  analysis  the  worst  case  of  thermocouple 
insertion  is  assumed,  namely,  that  the  bead  does  not  touch  the  bottom 
of  the  hole.  It  is  thus  assumed  that  the  thermocouple  bead  rests 
approximately  1/16  in.  from  the  bottom  of  the  hole.  This  allows 
application  of  specific  boundary  conditions  specified  by  the  previous 
models.  Further,  it  was  assumed  that  the  temperature  gradient  in  the 
wires  was  the  same  as  that  in  the  air.  In  actual  fact,  the  temperature 
gradient  will  be  less  than  that  of  the  air  due  to  the  lead  wire  thermal 
properties  as  compared  to  those  of  air.  The  equation31  for  finding  the 
conduction  heat  transfer  rate  was: 

^conduction  "  -K\  3^ 


where  —  is  the  temperature  gradient  along  the  wires,  is  the  wire 
cross-sectional  area  and  K  is  the  thermal  conductivity  of  the  wires. 

The  thermal  conductivity  used  was  the  average  of  that  for  chromel  and 
alumel32  since  type  K.  chromel-alumel  thermocouples  were  to  be  used  in 
testing.  The  cross-sectional  area  was  written  as  a  function  of  wire 
radius  so  that  a  program  could  be  written  leaving  radius  as  a  variable. 

The  convection  heat  transfer  rate  term  arises  due  to  the  heat 
transferred  from  the  surrounding  air  to  the  thermocouple  bead.  To 
determine  the  convection  heat  transfer  rate,  Newton's  law  of  cooling 


q 


convection 


(19) 
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was  used,  where  h  is  the  average  convection  heat  transfer  coefficient, 
is  the  thermocouple  bead  surface  area,  Tgir  is  the  temperature  of 
the  air  in  the  air  cavity  at  the  location  of  the  thermocouple  bead,  and 
is  the  thermocouple  bead  temperature.  In  order  to  be  conservative, 
the  average  convection  heat  transfer  coefficient  was  taken  to  be  the 
upper  limit  free-convection  value  given  by  Kreith.33  The  radius  of  the 
bead  was  selected  to  be  twice  the  lead  wire  radius  and  the  bead  surface 
area  was  then  reduced  to  a  function  of  radius. 

The  radiation  heat  transfer  rate  contribution  in  Equation  (17) 
arises  from  the  heat  transferred  from  the  air  cavity  sides  (carbon)  to 
the  thermocouple  bead.  Assuming  blackbody  radiation,  the  equation  for 
the  radiation  contribution  is: 


q  , .  .  .  =  A  F 

radiation  c  c-b 


=  A  F  ,  o(t  4  -  T  4\ 
c  c-b  ^  c  b  J 


where  A  is  the  area  of  the  carbon  that  is  radiating  heat  to  the  bead, 
c 

F  ,  is  the  radiation  shape  factor  from  the  carbon  boundary  to  the 
c-b 

thermocouple  bead,  and  a  is  the  Stefan-Boltzmann  constant,  Tc  is  the 
carbon  temperature,  and  T,  is  the  thermocouple  bead  temperature.  The 

D 

reciprocity  theorem33  in  this  application  is: 


A  F  , 
c  c-b 


where  A^  is  the  bead  surface  area  and  F^  ^  is  the  radiation  shape 
factor  from  the  thermocouple  bead  to  the  carbon  boundary.  Substituting 


Equation  (21)  into  Equation  (20)  -gives: 
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(22) 


^radiation  \  Fb-c  °  (TC4  "  Tb4) 

This  form  of  the  equation  is  desirable  since  can  be  reduced  to  a 
function  of  radius  and  since  ^  can  be  readily  approximated  as  the 
value  1;  all  radiation  leaving  the  bead  reaches  the  carbon  boundary. 
This  relation  along  with  Equations  (17),  (18),  and  (19)  can  be  sub¬ 
stituted  into  Equation  (16)  to  yield: 


As  shown  in  Equation  (16) ,  the  time  derivative  of  the  internal  energy 
can  be  written  for  a  solid  as: 


dU 

dt 


d_ 

dt 


pVu;  du  =  C  dT, 
b 


(24) 


where  p  is  the  bead  density,  V  is  the  bead  volume  which  is  reducible  to 
a  function  of  radius,  C  is  the  bead  heat  capacity,  and  is  the  bead 

temperature.  The  values  of  p  and  C  were  taken  as  average  values  of 
chromel  and  alumel  and  obtained  from  Benedict.32  With  the  assumption 
of  constant  p,  V  and  C,  the  time  derivative  of  U  becomes: 


dU 

dt 


(25) 


Using  the  Euler  time  marching  technique  along  with  Equations  (23)  and 
(25)  yields: 
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Tt  +  At  _  Tt 

PVC  ^  I+  *  \  (Tair  -  Tb) 

r'  ' 

+  \  Fb-c  "(V  -  Tb“)  (26) 

The  carbon  temperature,  Tc»  was  represented  by  an  average  of  the  first 
node  and  second  node  temperatures  resulting  from  the  one-dimensional 
carbon  model-  The  air  temperature  was  given  as  the  midline  node 
temperature  next  to  the  heated  surface  using  the  two-dimensional  air 
cavity  model  [Node  (2,  3)  in  Figure  6].  Equation  (26)  was  programmed 
on  a  digital  computer  to  give  the  thermocouple  bead  temperature  as  a 
function  of  time.  The  resulting  program  is  shown  in  Appendix  C. 

Figure  8  shows  the  effect  of  the  variation  of  the  radius  of  the 
thermocouple  lead  wires.  The  results  shown  in  this  figure  were  obtained 
assuming  no  heat  conduction  along  thermocouple  leads.  This  would  be 
the  case  if  a  thermocouple  hole  were  placed  so  that  no  temperature 
gradient  existed  along  the  lead  wires  of  the  thermocouple.  This  situa¬ 
tion  is  approximated  by  entering  a  brake  disk  in  the  radial  direction 
with  the  thermocouple.  It  should  be  reiterated  that  this  analysis 
assumes  that  the  thermocouple  is  placed  1/16  in.  from  the  bottom  of  the 
hole. 

Figure  9  shows  the  effect  if  conduction  heat  transfer  along  the 
lead  wires  is  taken  into  account.  It  should  be  recalled,  as  discussed 
under  Equation  (19) ,  that  at  all  times  the  bead  radius  was  twice  the 
lead  wire  radius.  It  is  noted  that  at  large  lead  wire  radii,  the 
effect  of  conduction  out  of  the  lead  wires  is  not  nearly  as  important 
as  the  fact  that  the  mass  of  the  bead  has  increased  significantly.  On 
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Figure  8.  Thermocouple  Model-Radius  Effect 


Figure  9.  Thermocouple  Model-Conduction  Effect 


the  other  hand,  at  low  lead  wire  radii  the  conduction  effect  is  signifi¬ 
cant,  especially  at  the  lower  temperatures.  At  high  temperatures  the 
radiation  effect  allows  the  bead  temperature  to  asymtotically  approach 
the  temperature  at  the  bottom  of  the  hole. 

Conclusions  From  Theoretical  Model 

Conclusions  that  can  be  drawn  from  this  section  of  the  temperature 
measurement  analysis  chapter  are  summarized  below. 

1.  In  the  end  stator  of  a  brake  with  a  thermocouple  entering  in 
the  axial  direction,  a  thermocouple  placed  1/16  in.  from  the 
bottom  of  the  hole  can  have  an  error  of  370°F  out  of  840°F. 
This  occurs  when  a  .02  in.  radius  lead  wire  is  used.  This 
error  decreases  at  higher  temperature  due  to  radiation 
effects  (300°F  error  at  1050°F). 

2.  If  the  lead  wire  radius  is  reduced  to  .0015  in.,  even  with 
the  bead  not  touching  the  parent  material,  the  error  can  be 
reduced  to  40°F  at  low  temperature  (200°F)  and  to  a  minimum 
error  of  18°F  at  high  temperature  (1050°F)  for  the  model 
considered. 

3.  In  case  2,  radiation  effects  insure  that  the  bead  temperature 
is  higher  than  the  air  temperature,  but  the  conduction 
through  the  lead  wires  and  the  convection  to  the  air,  both  of 
which  remove  heat  from  the  bead,  insure  that  the  parent 
material  temperature  can  never  be  reached. 

4.  The  three  conclusions  above  are  based  on  a  properly  made 

thermocouple  junction.  If  the  junction  is  twisted  a  few 

times  at  the  end,  the  temperature  measured  is  affected  in  two 
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ways.  First,  the  EMF  passing  through  the  thermocouple  system 
does  so  partially  at  each  contact  point  of  the  twist,  the 
highest  EMF  fraction  being  at  the  point  of  contact  of  least 
resistance.  This  point  is  unknown  so  the  temperature  measured 
will  probably  not  be  at  the  thermocouple  end.  Second,  the 
thermocouple  standard  EMF  tables  are  listed  for  welded 
contact  at  the  thermocouple  wire  junction.  If  contact  occurs 
at  other  locations,  as  in  a  twist,  the  EMF  passes  through  an 
oxide  layer  on  the  lead  wires,  rather  than  through  only  pure 
metal. 

The  radial  temperature  gradient  in  the  discs  is  less  than  the 
axial  gradient  in  the  end  stator.  Thus,  the  error  in  measure¬ 
ment  is  less  in  that  case  since  the  conduction  loss  is 
negligible. 

A  standard  Type  K  thermocouple  error  at  1000°F  is  ±7.5°F  and 
2000°F  is  ±15°F.  If  this  error  is  added  to  those  discussed 
above  (twists  instead  of  beads,  improper  placement,  large 
wire  diameter),  the  overall  error  would,  of  course,  be  even 
greater.  The  .75%  error  of  above  can  be  reduced  by  purchasing 
special  wire.  This  decreases  this  type  of  error  to  no  more 
than  .375%.  Calibration  should  also  be  accomplished. 

The  fact  that  no  bead  is  formed  in  some  applications  (twists 
only  are  sometimes  used)  probably  helps  the  situation  from 
the  standpoint  of  response  time.  The  main  factor  in  decreasing 
the  lead  wire  diameter  is  to  decrease  the  bead  size;  then  the 
response  time  is  less.  Thus,  if  no  bead  is  formed  and  EMF  is 


still  transmitted,  the  results  could  be  better  than  if  a  bead 
were  formed,  if  a  large  radius  lead  wire  is  used.  This  is 
not  recommended  due  to  other  errors,  as  stated  above. 

Bench  Tests 
Objective 

Bench  tests  were  conducted  to  experimentally  determine  the  dif¬ 
ferences  that  exist  between  various  types  and  sizes  of  test  thermo¬ 
couples  when  all  were  subjected  to  the  same  heat  input  within  the  same 
amount  of  time,  utilizing  the  same  test  specimen.  The  results  were  to 
be  used  in  related  tests.  It  was  desired  to  know  the  accuracy  of 
current  brake  laboratory  thermocouple  practice  so  that  changes  could  be 
made  if  required  for  conducting  the  full-scale  brake  tests.  It  is 
desirable  to  have  the  best  experimental  data  practical  to  eliminate 
that  possible  unknown  prior  to  the  comparison  with  the  computer  model. 

Test  Apparatus 

Test  Specimen  and  Heat  Lamps.  A  test  apparatus  with  banks  of  quartz 
lamps  was  used  for  heating  the  first  specimen.  The  lamp  control  is 
such  that  the  lamps  can  be  maintained  at  a  constant  temperature  auto¬ 
matically  or  varied  manually.  For  these  tests  the  control  was  manually 
set  to  the  maximum  heat  rate  and  maintained  there  throughout  all  but 
one  of  the  48  tests  conducted.  The  material  selected  for  heating  was 
allowed  to  attain  1500°F  at  which  point  the  lamps  were  manually  turned 
off. 

In  order  to  approximate  the  heat  input  to  a  piece  of  material 


within  an  actual  brake  without  conducting  a  brake  test,  a  small  piece 
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(.45  in.  x  1.9  in.  x  3.0  in.)  of  carbon  composite  brake  materia]  was 


suspended  between  two  high-density  radiant  heaters  containing  quartz 
lamps.  This  arrangement  is  shown  in  Figure  10.  A  carbon  composite 
material  was  selected  as  the  test  specimen  since  carbon  is  being  used 
in  more  and  more  brakes  for  Air  Force  aircraft.  Since  the  carbon  com¬ 
posites  of  interest  are  about  one-fourth  the  density  of  steel,  the 
temperature  gradients  in  the  carbon  will  be  higher.  Also  shown  in 
Figure  10  is  the  rod  that  supported  the  carbon,  the  test  thermocouple, 
and  the  baseline  thermocouple.  The  thermocouple  holes  in  the  carbon 
were  1.0  in.  deep.  The  holes  were  located  in  the  middle  of  the  thick¬ 
ness  direction  and  were  0.7  in.  inward  from  the  edge  in  the  1.9  in. 
dimension  direction,  making  the  thermocouples  0.5  in.  apart.  The  lamp 
and  thermocouple  arrangement  simulates  the  brake  condition  in  which 
energy  is  input  from  each  rubbing  surface  due  to  friction,  while  the 
thermocouples  enter  from  the  outside  or  inside  radius  of  the  brake 
disc.  The  baseline  thermocouple  was  cemented  into  place  with  Omega  CC 
high  temperature  sodium  silicate  cement  that  would  withstand  tempera¬ 
tures  up  to  2000°F.  The  maximum  temperature  encountered  during  any 
test  was  1602°F.  This  baseline  thermocouple  was  installed  so  that  it 
could  be  assured,  after  each  test,  that  the  same  maximum  temperature 
occurred  in  the  carbon  block  during  each  test.  A  photograph  of  the 
carbon  block  with  the  supporting  rod  and  two  thermocouples  is  shown  in 
Figure  11. 

Data  Recording  Equipment.  The  temperatures  of  interest  were  recorded 

using  thermocouples  in  an  arrangement  as  shown  in  Figure  12.  The 

copper/alumel  junction  and  the  copper/chromel  junction  were  made  using 
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THERMOCOUPLE 


MODEL  5208  HIGH 
DENSITY  MODULAR 
RADIANT  HLATER 
CONTAINING  6 
TUNGSTEN  FILAMENT 
QUARTZ  LAMPS, 

144  VOLTS,  7.2  KILO¬ 
WATTS,  50  AM  PS,  NET 
HEAT  FLUX  666 
WATTS  /IN2,  WATER 
AND  AIR  COOLED 


TEST  SPECIMEN 
SUPPORT  ROD 


Figure  10.  Sketch  of  Thermocouple  Test  Apparatus 
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terminal  blocks.  The  thermocouple  junction  was  a  variable  and  is 
discussed  in  the  next  section.  A  Houston  Instruments  X-Y-Y'  plotter 
was  used  to  record  the  thermocouple  voltage  output.  The  plotter  was 
calibrated  using  a  Brown  Portable  Potentiometer,  Model  126W3,  from 
Minneapolis-Honeywell  Regulator  Company.  Calibration  data  points  were 
at  every  two  millivolts  to  a  full  scale  of  50  millivolts.  The  milli¬ 
volts  were  converted  to  °F  by  using  a  table  of  temperature  versus 
millivolts  published  by  the  National  Bureau  of  Standards. 


Thermocouple  Description.  Thermocouples  in  these  tests  were  all  type 
K,  chromel-alumel.  The  variations  were  in  diameter,  configuration, 
type  of  junction,  and  distance  from  bottom  of  the  hole  in  the  carbon 
block.  Table  2  shows  the  variation  in  diameter,  configuration,  and 


TABLE  2 


THERMOCOUPLE  DESCRIPTION 

Test  No. 
(Code) 

Description 

15,  47, 

48 

.01  in.  dia.  stainless  sheath,  magnesium 
oxide  insulation, 

.0015  in.  wire  dia. ,  grounded  junction. 

Omega  No.  SCASS-010G-6 

18,  19, 
(El) 

23 

,0625  in.  dia.  stainless  sheath,  magnesium 
oxide  insulation, 

.01  in.  wire  dia.,  exposed  junction.  Omega 

No.  CASS-116E-12 

24,  27 
(Ul) 

.0625  in.  dia.  stainless  sheath,  magnesium 
oxide  insulation, 

.01  in.  wire  dia.,  ungrounded  junction. 

Omega  No.  CASS-116U-12 
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TABLE  2-continued 


Test  No. 

(Code) 

Description 

29,  30 
(G2) 

.0625  in.  dia.  stainless  sheath,  magnesium 
oxide  insulation, 

.01  in.  wire  dia.,  grounded  junction.  Omega 

No.  CASS-116G-12 

32,  34 
(3) 

.01  in.  wire  dia.  unsheathed  T/C  installed  in 
.0625  in.  dia., 

.25  in.  LG  two-hole  magnesium  oxide  insulators 

36,  38,  41 

.045  in.  x  .078  in.  size,  glass  braid  insulation, 

.02  in.  wire  dia.  also  covered  with  glass 
braid  insulation,  junction  twisted  3  times. 

40B 

Same  as  above  glass  braid  except  junction 
twisted  10  times 

42 

Same  as  above  glass  braid  except  junction 
twisted  once 

44 

Same  as  above  glass  braid  except  junction  spot 
welded 

45 

.030  in.  x  .053  in.  size,  glass  braid  insulation, 
.012  in.  wire  dia.  also  covered  with  glass 
braid  insulation,  junction  spot  welded 

type  of  junction.  Five  basic  types  are  shown,  i.e.,  .01  in.  diameter 

stainless  steel  sheath,  .0625  in.  diameter  stainless  steel  sheath, 

.045  in.  x  .078  in.  diameter  glass  braid,  .030  in.  *  .053  in.  diameter 

glass  braid  and  an  unsheathed  thermocouple  installed  in  a  two-hole 

magnesium  oxide  insulator.  The  .01  in.  diameter  stainless  steel 

sheathed  thermocouple  had  a  grounded  junction  as  opposed  to  ungrounded 

and  exposed  junctions.  These  three  junction  types  are  shown  in 
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Figure  13.  The  .0625  in.  diameter  stainless  steel  sheathed  thermo¬ 
couples  were  tested  with  exposed,  grounded,  and  ungrounded  junctions. 
The  unsheathed  thermocouple  was  installed  in  a  .0625  in.  outside 
diameter,  two -hole  magnesium  oxide  insulator.  The  insulator  was 
one  inch  long  at  the  end  placed  within  the  carbon  block,  while  the 
remainder  of  the  6  in.  thermocouple  was  covered  with  .25  in.  long 
pieces  of  the  two-hole  magnesium  oxide  insulator.  The  remaining 
thermocouples  were  such  that  the  individual  wires  were  insulated  with 
glass  braid  insulation  and  the  two  insulated  wires  were  then  insulated 
and  held  together  with  the  same  type  glass  braid  insulation.  The  spot- 
welded  thermocouples  were  constructed  by  overlapping  the  lead  wires, 
using  a  hand-held  spot  welding  tool  and  cutting  off  the  remaining  wire 
ends.  Several  thermocouples  were  joined  by  twisting  the  leads  to  form 
a  junction.  One  twist  is  defined  as  turning  the  leads  through 
180  degrees  while  holding  the  main  thermocouple  stationary.  The 
twisted  area  was  void  of  insulation  but  insulation  started  immediately 
at  the  spot  where  the  twists  stopped. 

Five  of  the  thermocouples  tested  were  independently  calibrated  at 
two  temperatures,  750°F  and  1100°F.  Appendix  D  shows  the  percent 
deviation  from  a  National  Bureau  of  Standards  tested  platinum,  10% 
rhodium  versus  platinum  thermocouple.  Also  shown  are  the  thermocouple 
codes  used  and  a  list  of  tests  during  which  the  various  thermocouples 
were  used.  In  addition.  Appendix  D  shows  the  calibration  of  the  same 
ten  thermocouples  but  gives  millivolts  and  °F  instead  of  percent.  Also 
included  are  the  conditions  at  which  the  thermocouples  were  calibrated, 
as  well  as  the  equipment  used  and  the  accuracy  of  the  results. 


Test  Procedure 


The  basic  test  procedure  was  to  subject  the  carbon  block  with 
embedded  thermocouples  to  the  maximum  lamp  output,  raising  the  baseline 
thermocouple  temperature  from  100°F  to  1500°F.  The  method  of  obtaining 
the  100®F  starting  temperature  for  all  tests  was  either  to  heat  the 
block  to  about  150°F  and  allow  it  to  cool  down  to  100°F  or  to  allow  the 
block  to  cool  down  to  100°F  from  the  previous  test.  The  heat  input  of 
each  test  was  manually  terminated  at  1500°F.  Due  to  this  manual 
control  the  peak  temperature  varied  slightly.  The  temperatures  during 
the  heat  up  and  cool  down  periods  were  recorded  on  the  X-Y-Y'  plotter. 

A  detailed  list  of  the  steps  involved  in  running  a  particular  thermo¬ 
couple  test  is  given  in  Appendix  D. 

Results  and  Discussion 

The  specific  baseline  and  test  thermocouple  output  in  the  form  of 
millivolts  versus  time  is  depicted  in  Figure  14  for  a  typical  test. 

The  particular  test  from  which  this  data  was  generated  was  Test  27. 

The  tests  were  run  in  numerical  order.  Not  all  tests  run  are  reported 
since  many  tests  were  duplicates  of  those  reported  and  produced  exactly 
the  same  results.  These  tests  were  run  specifically  to  assure  repeata¬ 
bility  from  test  to  test.  Table  3  shows  test  temperatures  at  various 
baseline  temperatures.  The  particular  thermocouples  associated  with 
each  test  number  are  given  in  Table  2.  Table  3  shows  that  five  tests 
were  run  in  which  the  thermocouple  was  placed  1/16  in.  from  the  bottom 
of  the  1.0  in.  deep  hole  in  the  carbon  block.  Table  2  shows  that 
Tests  15,  47,  and  48  were  conducted  with  the  same  thermocouple.  Test  15 
utilized  a  .02  in.  diameter  hole  in  the  carbon  block.  Test  47  used  the 
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TEST  TEMPERATURE  AT  VARIOUS  BASELINE  TEMPERATURES 


\  BASELINE 

TEST  N.  TEMP 
NO  \f°F) 

300 

450 

750 

1,100 

1,350 

1,500 

15 

316 

468 

778 

1,130 

1,381 

1,521 

18a 

223 

377 

725 

1,105 

1,361 

1,505 

19 

286 

448 

777 

1,144 

1,402 

1,539 

23 

274 

448 

777 

1,146 

1,402 

1,540 

24a 

212 

388 

736 

1,112 

1,366 

1,503 

27 

221 

395 

745 

1,124 

1,381 

1,519 

29a 

225 

412 

751 

1,133 

1,386 

1,524 

30 

283 

444 

777 

1,152 

1,408 

1,538 

32a 

217 

369 

731 

1,116 

1,374 

1,517 

34 

260 

423 

762 

1,147 

1,407 

1,542 

36a 

143 

241 

555 

985 

1,257 

1,396 

38 

160 

265 

602 

1,019 

1,278 

1,429 

40B 

139 

230 

537 

935 

1,187 

1,320 

41 

145 

245 

581 

1,019 

1,281 

1,425 

42 

207 

314 

640 

1,052 

1,306 

1,473 

44 

225 

363 

697 

1,089 

1,359 

1,486 

45 

232 

384 

733 

1,120 

1,385 

1,521 

47 

318 

469 

775 

1,136 

1,377 

1,516 

48 

323 

1 

478 

778 

1,133 

1,366 

1,512 

NOTE:  a. 

.  This 

thermocouple  was 

placed  1/16 

in.  from  hole 

bottom.  All  others  were  placed  at  the  hole  bottom. 
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same  small  thermocouple  but  it  was  cemented  inside  a  .0625  in.  x 
1.0  in.  long  stainless  steel  sheath  and  cemented  into  the  bottom  of  the 
.07  in.  diameter  hole  in  the  carbon  block.  This  was  done  to  negate  the 
possible  influence  of  the  metal  sheaths  of  other  thermocouples  when 
making  comparisons.  Comparison  of  values  from  Table  3  for  Tests  15  and 

47  shows  very  little  sheath  effect.  Test  48  differed  from  47  in  that 

48  had  a  heat  input  rate  approximately  one-half  that  of  Test  47.  All 
heating  rates  were  the  same  except  for  Test  48.  To  reach  1500°F  each 
test  required  approximately  45  to  50  seconds  except  for  Test  48,  which 
required  65  to  70  seconds.  This  heating  rate  is  representative  of  that 
of  a  typical  3  or  4  rotor  aircraft  brake  with  a  carbon  heat  sink  and  a 
deceleration  of  about  9  ft/sec2  with  an  absorbed  energy  of  about 

11  million  ft-lbs.  Table  2  shows  that  Tests  18,  19,  and  23  were  con¬ 
ducted  with  the  same  thermocouple.  Tests  19  and  23  were  identical 
runs;  Test  18  was  conducted  with  the  thermocouple  1/16  in.  from  the 
bottom  of  the  hole  in  the  carbon.  Tests  36  and  38  used  the  same 
thermocouple,  while  Test  41  used  a  physically  different  thermocouple 
but  one  configured  the  same  as  for  Test  38.  This  was  accomplished  due 
to  the  large  differences  observed  between  these  tests  and  other  tests 
of  thermocouples  enclosed  in  metal  sheaths.  Table  3  shows  that  tempera¬ 
tures  in  Tests  38  and  41  are  close,  especially  at  the  higher  three 
temperatures. 

As  previously  discussed,  the  baseline  thermocouple  was  installed 
to  insure  that  the  same  temperature  was  reached  in  the  carbon  block 
from  test  to  test.  However,  due  to  (1)  room  temperature  and  humidity 
changes  from  day  to  day,  (2)  the  overall  test  apparatus  heating  up 
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slightly  during  a  test,  and  (3)  specimen  weight  loss  from  test  to  test 
(9  grams  were  lost  of  the  total  specimen  weight  of  64  grams  over  the 
last  32  tests  due  to  oxidation  of  the  carbon) ,  the  baseline  temperature 
was  not  strictly  constant  over  all  tests.  The  variation  was  as  high  as 
3 Z  near  the  peak  temperature.  To  determine  the  effect  of  this  varia¬ 
tion  on  test  thermocouple  temperatures.  Tests  19  and  23  can  be  com¬ 
pared.  The  baseline  temperature  difference  was  the  only  variable 
between  the  two  tests.  As  shown  in  Table  3  or  Table  4,  this  effect  is 
extremely  small.  The  maxiumura  difference  was  at  a  baseline  temperature 
of  300°F.  In  this  area  the  slope  of  the  temperature  versus  time  curve 
is  the  greatest,  which  adds  inaccuracy  to  reading  the  temperature  in 
this  area.  Since  all  of  the  other  five  readings  are  so  close  and  the 
possibility  of  resolution  error  at  low  temperatures  is  high,  it  is 
reasonable  to  conclude  that  the  effect  of  a  slight  baseline  temperature 
change  from  test  to  test  is  negligible. 

Since  Test  47  was  accomplished  with  the  smallest  thermocouple  wire 
combined  with  a  .0625  in.  x  1.0  in.  long  stainless  steel  sheath  and  a 
cemented  junction,  this  test  was  judged  to  be  the  most  accurate  of  all 
those  conducted.  Due  to  this,  results  are  tabulated  and  graphed 
showing  the  percent  difference  of  the  test  thermocouple  from  thermo¬ 
couple  Number  47.  Table  4  gives  such  values  for  each  test  reported. 
Also  shown  is  the  calculation  procedure.  The  numbers  range  from  -56.3% 
to  +2.3%.  A  positive  deviation  indicates  that  the  test  thermocouple 
associated  with  that  positive  deviation  responds  faster  than  the 
Test  47  thermocouple,  or  that  it  is  more  accurate  than  the  Test  47 

thermocouple  or  a  combination  of  both.  Test  47  was  selected  for  the 
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TABLE  4 

PERCENT  DIFFERENCE  OF  TEST  THERMOCOUPLE 
FROM  THERMOCOUPLE  NUMBER  47 


-29.9  -19.6 


-12.1 

-3.1 

-11.0 

-5.3 

.3 

-31.8 

-21.3 

-5.7 

- 

* 

1.8 

-.22 

.1 

NOTE(S):  a.  Table  values  calculated  by  100  *  ^Ttest  ~  ^47^^47 

b.  This  thermocouple  was  placed  1/16  in.  from  the  hole 
bottom.  All  others  were  placed  at  the  hole  bottom.  A 7 


comparison  test  because  of  its  size  and  insertion  technique,  therefore 
this  small  value  (2.3%  maximum)  is  attributed  in  most  cases  to  small 
discrepancies  in  setup,  input  energy,  or  plotter  chracteristics.  The 
exception  to  this  is  Test  48  which  used  the  same  thermocouple  and 
installation  as  Test  47.  Test  48  is  the  only  test  to  indicate  positive 
values  in  Table  4  in  the  lowest  three  temperature  columns.  Since  the 
temperature  comparisons  of  Tables  3  and  4  are  based  on  baseline  tempera¬ 
ture,  the  test  with  a  lower  heat  input  rate  (Test  48)  can  give  a  higher 
temperature  than  the  test  with  a  higher  input  rate  (Test  47).  Although 
this  is  noted,  the  highest  temperature  difference  is  only  1.88%. 

Test  15  varies  less  than  1%  from  Test  47,  showing  the  lack  of  effect  of 
inserting  the  thermocouple  in  a  .02  in.  diameter  hole  without  cement  as 
in  Test  15  compared  with  using  the  same  thermocouple  cemented  in  ide  a 
.0625  in.  diameter  sheath  which  was,  in  turn,  cemented  in  the  hole  in 
the  carbon  as  in  Test  47.  This  small  difference  can  also  be  seen  in 
Figure  15  where  six  different  thermocouples  are  compared.  The  test 
numbers  are  shown  along  with  the  thermocouple  description.  A  more 
complete  description  is  given  in  Table  2.  As  noted  in  Figure  15,  all 
of  these  thermocouples  were  placed  at  the  bottom  of  the  hole  in  the 
carbon  except  for  Test  36.  From  Table  4  it  is  noted  that  results  of 
Tests  38  and  41  are  fairly  close,  thus  only  Test  38  is  shown  on 
Figure  15.  Recall  that  38  and  41  used  physically  different  thermo¬ 
couples  but  the  same  configuration.  Also  from  Table  4,  Tests  19  and  30 
show  approximately  the  same  values.  Again,  only  Test  19  is  graphed  on 
Figure  15.  The  difference  between  the  thermocouples  of  Test  19  and 

Test  30  is  that  Test  19  had  an  exposed  junction  and  Test  30  had  a 
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Figure  15.  Accuracy  of  Different  Thermocouple  Types 


grounded  junction  (see  Figure  13).  One  obvious  and  predictable  trend 
is  that  the  difference  increases  in  absolute  value  as  the  thermocouple 
diameter  increases.  The  exposed  junction  (and  thus  the  grounded 
junction)  is  better  than  the  ungrounded  junction. 

A  graph  depicting  the  differences  of  Table  4  for  thermocouples 
with  glass  braid  insulation  is  shown  in  Figure  16.  The  graph  shows 
that  10  twists  are  less  accurate  than  3  twists  which  are,  in  turn,  less 
accurate  than  1  twist.  The  spot-welded  thermocouples  are  most  accurate 
and  the  smallest  diameter  one  is  best.  It  should  be  pointed  out  that 
the  spot  welding  was  done  with  a  hand-heJd  spot  welder,  and  as  a 
result,  it  was  impossible  to  make  a  perfect  bead.  This  could  account 
for  the  up  to  25%  difference  shown  at  low  temperatures.  Any  better 
method  to  form  the  bead  on  this  type  thermocouple  was  not  attempted 
since  it  would  be  time  consuming  if  reduced  to  everyday  practice. 

A  comparison  should  be  made  between  Test  45  and  Test  19,  both  of 
which  are  shown  in  Figure  17.  The  first  has  glass  braid  insulation  and 
the  second  has  a  stainless  steel  sheath,  but  the  wire  diameters  are 
almost  identical.  The  major  difference  is  in  the  bead  formation.  The 
thermocouple  of  Test  19  was  made  under  controlled  conditions  by  the 
manufacturer,  whereas  the  other  was  made  by  spot  welding.  Observing 
the  two  graphs  (Figures  16  and  17)  or  the  tabular  data  in  Table  4  leads 
to  the  conclusion  that  the  differences  in  test  thermocouple  45  are 
about  triple  those  of  Test  19  (or  Test  30).  The  best  thermocouple, 
other  than  the  small  diameter  one,  is  the  one  used  in  Test  19  or 
Test  30. 
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Figure  16.  Accuracy  of  Glass  Insulation  Thermocouples 


Figure  17.  Thermocouple  Comparison/Welding  Method 


The  effect  of  placing  the  thermocouple  1/16  in.  from  the  bottom  of 
the  hole  can  be  seen  in  Figure  18.  In  each  case,  as  expected,  the  most 
accurate  is  the  one  at  the  bottom  of  the  hole.  Glass  braid  insulated 
ones  are  least  accurate  followed  by  ungrounded  and  exposed  (or 
grounded,  since  results  of  Tests  19  and  30  are  almost  identical).  Note 
that  the  glass  braid  insulated  thermocouples  have  been  shown  to  have 
higher  differences  from  the  Test  47  thermocouple  not  because  of  the 
insulation  type,  which  is  merely  a  categorizing  method,  but  because  of 
the  junction  on  those  particular  thermocouples.  Also  of  interest  in 
this  area  is  the  relative  condition  of  the  thermocouples  after  testing. 
The  metal  sheathed  thermocouples  showed  no  change  but  those  with  glass 
braid  insulation  deteriorated  significantly.  After  one  exposure  to 
1500°F  during  a  test,  the  insulation  became  white  and  brittle.  Any 
handling  after  this  one  exposure  in  most  cases  led  to  the  insulation 
crumbling  and  falling  off.  After  about  five  tests  with  any  one  thermo¬ 
couple  the  insulation  would  fall  off  when  the  thermocouple  was  removed 
from  the  test  specimen.  Although  nothing  irregular  was  noted  in  the 
data,  this  phenomenon  could  cause  the  thermocouple  to  form  a  junction 
in  any  location  along  the  white  brittle  area. 

A  plot  of  temperature  versus  time  for  five  thermocouples,  including 
the  thermocouple  used  in  Test  47,  is  shown  in  Figure  19.  A  significant 
lag  is  seen  in  Tests  36,  38,  and  40B  when  they  are  compared  to  Test  47. 
The  glass  braid  insulation  thermocouple  with  ten  twists  is  the  most 
dramatic  followed  by  the  three-twist  tests.  The  best  thermocouple 
appears  to  be  that  of  Test  30.  It  is  shown  to  be  above  the  Test  47 
curve  at  the  high  temperatures.  This,  again,  is  attributed  to  small 


Figure  18.  Effect  of  Thermocouple  Distance  from  Hole  Bottom 


Temperature  vs  Time,  Thermocouple  Comparison 


discrepancies  in  setup,  input  energy,  or  plotter  characteristics.  Note 
as  before  that  Test  30  and  Test  19  show  almost  the  same  results.  This 
shows  that  either  the  exposed  (Test  19)  or  the  grounded  (Test  30)  metal 
sheathed  thermocouple  is  the  best  tested  other  than  the  Test  47  thermo¬ 
couple.  The  small  diameter  (.10  in.  sheath  diameter)  thermocouple  used 
in  Test  47  is  not  recommended  for  use  in  actual  brake  tests  due  to  its 
size  and  fragility.  The  highest  temperature  plotted  in  Figure  19  for 
each  curve  is  the  actual  peak  temperature  reached  during  the  test. 

This  is  somewhat  misleading  since  baseline  peak  temperatures  for  these 
tests  varied  slightly. 

A  comparison  of  the  peak  temperatures  reached  during  all  reported 
tests  is  shown  in  Table  5.  The  peak  temperatures  of  the  test  thermo¬ 
couples  are  not  directly  correlatable  since  the  baseline  peak  tempera¬ 
tures  are  not  identical.  To  circumvent  this  difficulty,  the  peak 
temperature  of  each  test  thermocouple  was  adjusted  by  multiplying  it  by 
the  ratio  of  the  Test  47  baseline  peak  temperature  to  the  baseline  peak 
temperature  of  the  particular  test  in  question.  This  adjustment  is 
shown  in  the  third  column  of  temperatures  in  Table  5.  The  column 
following  shows  the  percent  difference  between  each  adjusted  test  peak 
temperature  and  the  peak  temperature  of  Test  47.  The  greatest  dif¬ 
ference  is  for  Test  40B  which  shows  that  the  adjusted  peak  temperature 
of  Test  40B  is  11.1%  lower  than  that  of  Test  47.  This  is  the  glass 
braid  insulation  thermocouple  with  10  twists  as  the  junction.  The  next 
highest  difference  is  with  Test  36  which  is  5.1%  low.  This  thermo¬ 
couple  was  of  the  glass  braid  insulation  type  with  3  twists  forming  the 
junction  and  was  placed  1/16  in.  from  the  bottom  of  the  hole  in  the 
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TABLE  5 — continued 


TABLE  5 — continued 


carbon  block.  All  other  adjusted  peak  temperatures  are  within  3%  of 
that  of  Test  47.  The  last  column  in  Table  5  shows  the  difference 


between  the  time  the  adjusted  peak  temperature  of  each  test  occurred 
and  the  time  the  peak  temperature  of  Test  47  occurred.  This  informa¬ 
tion  could  only  be  read  accurately  to  the  nearest  1/4  sec.  The  largest 
lag  occurred  with  Tests  36,  41,  and  42  at  2.75  sec  followed  closely  by 
Tests  38  and  40B  at  2.5  sec.  These  five  tests  were  all  of  the  .061  in. 
diameter  glass  braid  insulation  type  with  twisted  junctions.  The 
.061  in.  diameter  glass  braid  with  the  welded  junction  (Test  44)  was 
next  at  1.25  sec.  Thermocouples  placed  1/16  in.  from  hole  bottom, 
other  than  the  .061  in.  diameter  glass  braid,  had  lag  times  of  from  .75 
to  1.0  sec.  The  three  tests  (15,  47,  48)  conducted  with  the  .01  in. 
diameter  sheath  had  no  lag  time.  Comparing  the  two  spot-welded  thermo¬ 
couples  shows  that  the  smaller  diameter  lead  wire  (.012  in.  vs  .02  in.) 
gives  a  smaller  lag  time,  0.5  vs  1.25  sec.  The  best  thermocouples  as 
far  as  supplying  the  least  lag  time  (.25  sec),  other  than  the  .01  in. 
diameter  sheath  thermocouple,  are  those  of  Tests  19,  23,  30,  34.  All 
of  these  are  either  grounded  or  exposed. 

Conclusions 

These  conclusions  are  based  on  the  particular  test  configuration 
employed.  If  other  applications  change  heating  rates  or  levels  or 
materials,  the  levels  of  error  could  change.  Each  conclusion  is 
followed  by  a  set  of  numbers  which  refers  to  the  tests  from  which  the 
conclusion  was  drawn. 


1. 


The  difference  in  temperature  caused  by  increasing  the  hole 


in  the  carbon  block  from  .02  in.  diameter  to  .07  in.  diameter  and 
cementing  a  stainless  steel  sheath  in  the  larger  hole  is  negligible 
(15,  47). 

2.  The  difference  in  temperature  caused  by  decreasing  the  heat 
input  rate  to  approximately  half  is  less  than  2%  for  the  .01  in. 
diameter  sheath  thermocouple  (47,  48). 

3.  The  method  of  creating  a  junction  by  twisting  the  thermo¬ 
couple  end  gives  temperatures  up  to  35%  low  if  1  twist  is  used  and  up 
to  55%  low  if  10  twists  are  used  (38,  40B,  41,  42). 

4.  Welding  the  thermocouple  junction  reduces  that  maximum  error 
from  the  above  value  of  35%  for  1  twist  to  29%  (42,  44). 

5.  Using  a  smaller  diameter  thermocouple  (.012  in.  vs  .02  in.) 
reduces  the  error  about  5%  (44,  45). 

6.  The  error  incurred  in  moving  the  thermocouple  1/16  in.  from 
the  bottom  of  the  hole  in  the  carbon  is  approximately  5%  (18,  19,  24, 
27,  36,  38). 

7.  The  best  two  thermocouples  are  the  .0625  in.  diameter  stain¬ 
less  steel  sheathed  thermocouples,  one  of  which  had  an  exposed  junction 
and  the  other  had  a  grounded  junction.  The  grounded  junction  is  recom¬ 
mended  in  the  brake  application  due  to  the  lack  of  junction  protection 
dn  the  exposed  junction  thermocouple.  The  smallest  diameter  thermo¬ 
couple  (.01  in.)  is  not  recommended  due  to  its  size  and  thus  inherent 

fragildty  (19,  30).  Refer  to  Table  2,  page  38,  for  dead  wire  diameters. 

8.  The  use  of  glass  braid  insulated  thermocouples  at  tempera- 
",i,i  around  or  abuve  1500°F  should  be  limited  to  short-time  usage  in 


nonmoving  environments  due  to  the  insulation  brittleness  and  probable 
crumbling  (36,  38,  40B,  41,  42,  44,  45). 

9.  As  long  as  an  excessive  (above  3)  number  of  twists  is  not 
used  and  the  thermocouple  is  located  at  the  bottom  of  the  hole,  the 
error  in  peak  temperature  is  less  than  3%  (36,  40B) . 

10.  A  lag  time  of  up  to  2.75  sec  can  be  expected  at  the  peak 
temperature  for  twisted  junction  thermocouples  (36,  38,  40B,  41,  42). 
This  lag  time  can  be  reduced  to  .25  sec  for  exposed  or  grounded  junc¬ 
tion  stainless  steel  sheathed  thermocouples  (18,  23,  30). 

11.  All  of  the  ten  thermocouples  calibrated  independently  from 
the  manufacturer  proved  to  be  more  accurate  than  the  tolerance  given  by 
the  manufacturer. 
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CHAPTER  IV 


FULL-SCALE  BRAKE  TESTS 

Introduction 

At  present,  in  the  evaluation  of  brakes  for  automobiles,  trucks, 
and  aircraft,  a  fairly  extensive  testing  program  is  essential.  It  is 
time  consuming  and  expensive  to  conduct  full-scale  brake  tests  on  the 
vehicles  they  are  designed  to  stop.  Although  the  ultimate  evaluation 
must  occur  on  the  particular  vehicle,  much  laboratory  testing  is 
necessary  first.  Two  basic  test  methods  are  generally  used  to  evaluate 
brakes  in  the  laboratory.  The  first  and  most  common  method  used  to 
test  brakes  uses  a  flywheel  to  provide  the  relative  speed  between  the 
brake-axle  assembly  and  the  flywheel.  Thus,  the  brake-axle  assembly 
remains  stationary  and  the  flywheel  turns.  The  speed  at  the  flywheel's 
outer  diameter  is  equal  to  that  speed  required  between  the  test  vehicle 
and  the  ground.  The  inertia  of  the  flywheel  is  usually  variable  to 
allow  the  vehicle  mass  to  be  simulated.  This  testing  method12  is 
universally  used  within  the  brake  industry.  In  the  second  method,  a 
machine  employing  a  flat  plate  creates  relative  motion  between  the 
plate  and  the  brake-axle  assembly.  This  test  method  is  seldom  used 
since  speeds  attainable  are  limited  due  to  power  requirements,  moving 
plate  weight  and  plate  length.  Additionally,  the  plate  dimension  in 
the  moving  direction  can  restrict  the  duration  of  the  brake  test.  If 
the  plate  is  made  longer,  available  space  becomes  a  primary 
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concern.  The  first  method  discussed  eliminates  many  of  these  problems. 

On  the  other  hand,  test  flywheels  or  dynamometers  are  expensive  to 

i 

procure  and  maintain  and  are  thus  found  mostly  in  large  corporations  or 
governmental  facilities.  This  is  the  primary  reason  that  many  efforts 

I 

dealing  with  brake  theoretical  models  are  void  of  experimental  verifi¬ 
cation  data. 

Test  Apparatus 

Dynamometer 

A  typical  dynamometer  is  shown  in  Figure  20.  The  flywheel  is  in 
three  sections.  The  center  section  rotates  and  the  outer  two  sections 
remain  stationary.  Each  of  the  outer  sections  is  a  series  of  discs 
that  can  be  transferred  to  the  center  rotating  section  prior  to  a  test. 

This  allows  different  weights  of  vehicles  to  be  simulated  from*  test  to 
test.  The  wheel  assembly,  as  shown  in  Figure  20,  is  installed  on  the 
axle  which  is,  in  turn,  fastened  rigidly  to  a  structure  reaching  to  the 
floor  and  supported  by  large  rollers.  This  structure,  as  well  as  the 
wheel  assembly  is  pneumatically  moved  toward  the  flywheel  during  a 
test.  All  of  the  full-scale  brake  tests  were  conducted  on  a  192  in. 
diameter  dynamometer  located  within  the  Landing  Gear  Test  Facility, 

Building  31,  in  Area  B  of  Wright-Patterson  Air  Force  Base.  This  is  the 
largest  dynamometer  in  the  world.  The  size  of  the  flywheel  is  important 
to  more  adequately  simulate  a  flat  road  surface  (which  has  an  infinite 
radius  of  curvature)  and  in  being  able  to  have  large  mass  capabilities. 

This  192  in.  dynamometer  has  the  capability  to  test  wheel  assemblies  on 
both  sides  and  the  small,  or  north  side,  carriage  was  used  in 
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this  research.  This  carriage  has  the  following  maximum  capabilities: 

(1)  40,000  lbf  normal  load  between  the  tire  and  flywheel 

(2)  200  miles/hour  tangential  velocity  at  the  flywheel  maximum 
radius 

(3)  205  x  106  ft-lbf  kinetic  energy 

(4)  80  in.  tire  diameter, 

(5)  inertia  equivalent  of  162,987  lbm 

The  inertia  equivalent  of  the  flywheel  is  obtained  for  a  particular 
test  by  equating  the  desired  kinetic  energy  to  the  kinetic  energy  of 
the  flywheel  in  angular  terms  and  solving  for  the  inertia  equivalent 
term.  The  following  equation  applies. 

Vf 2  - 1 IE  h2  “f2  (27) 

where  M  is  the  mass  of  the  vehicle  to  be  tested,  is  the  velocity  at 

which  the  vehicle  is  traveling,  I  is  the  flywheel  inertia,  is  the 

n  t 

flywheel  angular  velocity,  IE  is  the  inertia  equivalent  and  is  the 
flywheel  outer  radius.  Solving  for  the  inertia  equivalent  gives: 

IE  =  M Vv2y/ Rf2wf2  (28) 

Full-Scale  Brakes 

The  brake  selected  for  the  experimental  tests  is  also  used  in  the 
F-14  aircraft.  A  four-rotor  brake  was  primarily  selected  because  two 
sets  of  rotors  and  stators  were  readily  available  that  would  fit  into 
the  brake  housing.  One  complete  set  of  rotors  and  stators  for  a  brake 
is  called  a  heat  sink.  The  above  heat  sink  selection  allowed  the 
testing  of  two  heat  sink  materials  without  the  effects  of  a  different 
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housing,  tire,  and  wheel.  The  two  heat  sink  materials  used  were  carbon 
composite  and  beryllium.  This  choice  allowed  a  variation  in  density, 
heat  capacity,  and  thermal  conductivity.  A  typical  beryllium  rotor  is 
shown  in  Figure  21.  The  keyways  are  located  at  the  outside  diameter  of 
the  rotor.  The  rotor  keys  to  the  rotating  wheel  through  these  keyways 
and  thus  rotates  with  the  wheel.  Wear  pads  are  bolted  to  the  beryllium 
core  as  shown.  The  wear  pad  thickness  is  about  3%  of  the  total  disc 
thickness.  Figure  22  shows  a  typical  beryllium  stator.  Here  the 
keyways  are  at  the  inside  diameter  since  the  stator  is  keyed  to  the 
nonrotating  axle.  Note  that  the  force  must  be  greater  at  a  stator 
keyway  than  at  a  rotor  keyway  since  both  support  the  same  torque  during 
a  test.  Thus,  in  order  to  have  a  more  even  force  distribution,  the 
stator  has  ten  keyways  compared  to  eight  for  the  rotor.  Figures  23  and 
24  are  photographs  of  a  typical  carbon  composite  rotor  and  stator, 
respectively.  Other  than  the  material,  the  main  difference  between  the 
carbon  and  beryllium  is  that  the  carbon  requires  no  separate  wear  pads, 
as  does  the  beryllium.  The  separate  wear  pads  are  required  on  beryllium 
since  any  rubbing  on  beryllium  itself  creates  a  toxic  dust. 

Figure  25  shows  a  photograph  of  the  beryllium  heat  sink  ready  to 
be  assembled.  The  brake  housing  is  shown  on  the  left,  the  heat  sink 
(all  rotors  and  stators)  is  on  the  right,  and  the  torque  tube,  through 
which  the  stators  are  keyed  to  the  axle,  is  in  the  middle.  At  assembly 
the  housing  is  bolted  to  the  left  side  of  the  torque  tube,  and  the 
stator  on  the  far  right  is  bolted  to  the  right  end  of  the  torque  tube. 
This  allows  the  hydraulic  pressure  to  push  the  pressure  plate  (plate  at 
the  right  of  the  housing)  to  the  right,  engaging  the  stators  and 
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Typical  Berylli 


Figure  24.  Typical  Carbon  Composite  Stator 


rotors.  The  carbon  heat  sink  is  assembled  similarly  using  the  same 
housing  and  torque  tube.  The  assembled  brake  is  shown  in  Figure  26. 

The  torque  tube  can  be  seen  at  the  housing  inside  diameter.  The  ten 
cylindrical  objects  arranged  in  sets  of  two,  placed  symmetrically 
within  the  housing,  are  the  return  springs.  They  insure  that  the 
pressure  plate  will  be  pulled  away  from  the  stators  and  rotors  when  the 
brake  hydraulic  pressure  is  not  applied.  This  allows  the  stators  and 
rotors  to  separate  from  one  another  when  no  braking  is  desired.  The 
five  lighter  colored  cylindrical  shapes  (Figure  26)  evenly  spaced 
between  the  sets  of  return  springs  are  the  pistons  to  which  the  hydrau¬ 
lic  fluid  is  fed  to  apply  the  pressure  required  during  vehicle  braking. 
The  hydraulic  fluid  enters  the  housing  through  the  fluid  port  shown  on 
the  extreme  left  of  Figure  26.  The  fluid  is  channeled  to  each  piston 
through  the  raised  channels  shown  between  each  pair  of  pistons. 

Next,  the  brake  is  mounted  on  the  axle  or  mandrel  in  front  of  the 
dynamometer  as  shown  in  Figure  27.  The  rotating  portion  of  the  fly¬ 
wheel  is  shown  in  the  background  on  the  right  with  the  nonrotating 
portion  on  the  left.  The  structure  into  which  the  mandrel  fits  is 
mounted  on  large  rollers  and,  as  mentioned  earlier,  will  move  toward 
the  flywheel  during  a  test.  The  hydraulic  pressure  line  is  shown  at 
the  top  of  the  brake.  The  wires  with  the  connectors  are  the  thermo¬ 
couples  used  for  temperature  measurement  and  will  be  described  in  more 
detail  later. 

The  last  step  in  assembling  the  test  setup  is  to  install  the  tire 
and  wheel  over  the  brake  and  to  attach  the  end  axle  nut  as  shown  in 
Figure  28.  The  assembly  in  front  of  the  wheel  including  the  wires  is 
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Figure  27.  Beryllium  Brake  on  Test  Mandrel 


Figure  28.  Beryllium  Brake  With  Tire  and  Wheel  Assembly 


part  of  the  instrumentation  used  and  will  be  described  in  the  instru¬ 
mentation  section. 


Test  Plan 

The  basic  objective  was  to  test  both  the  carbon  and  beryllium 
brakes  to  the  same  set  of  conditions.  The  inertia  equivalent,  normal 
load,  rolling  radius  (distance  from  flywheel  outer  surface  to  the 
mandrel  centerline) ,  kinetic  energy,  and  the  inflation  pressure  were 
selected  to  be  approximately  those  of  the  F-14  aircraft.  The  inertia 
equivalent  was  10,147  lbm  and  was  calculated  as  described  earlier.  The 
normal  load  for  all  stops  was  20,000  lbf,  and  the  rolling  radius  was 
16.5  in.  The  inflation  pressure  was  not  an  independent  parameter  and 
was  selected  so  that  the  20,000  lbf  normal  load  would  result  in  a 
rolling  radius  of  16.5  in.  as  specified.  Tests  were  run  at  two  kinetic 
energies,  i.e.,  11.23  *  106  ft-lbf  and  6.74  x  106  ft- lbf.  These  are 
subsequently  referred  to  as  100%  and  60%  conditions,  respectively. 

Since  the  kinetic  energy  and  mass  were  specified,  the  initial  velocity 
was  calculated.  This  resulted  in  141  mph  for  the  60%  kinetic  energy 
tests  and  182  mph  for  the  100%  kinetic  energy  tests.  The  brake  pressure 
was  maintained  at  a  constant  value  during  a  particular  test  and  was 
selected  based  upon  a  desired  flywheel  deceleration.  Two  decelerations 
were  selected  in  order  to  have  two  heating  rates  for  the  brake.  Since 
the  deceleration  is  dependent  on  the  wearing  surface  coefficient  of 
friction  characteristics,  the  actual  deceleration  varies  sompwhat  from 
the  desired  values.  The  first  value  selected  was  3.0  ft/sec2  and  was 
selected  based  on  a  smooth  routine  F-14  landing.  The  second  was 
7.0  ft/sec2  and  was  selected  based  on  a  more  severe  braking  condition 
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but  not  at  the  maximum  capability  of  the  brake.  Based  on  previous  test 
experience,  the  brake  pressures  could  now  be  specified  since  the 
dynamometer  operator  must  have  a  pressure  value  rather  than  a  decelera¬ 
tion.  For  the  beryllium  brake  the  pressures  required  to  give  approxi¬ 
mately  3.0  and  7.0  ft/sec2  deceleration  were  270  psig  and  400  psig, 
respectively.  The  carbon  brake  required  240  and  425  psig  to  attain  the 
desired  decelerations. 

Basically  two  types  of  tests  were  conducted.  The  first  was  a 
constant  pressure  stop  at  the  various  conditions  described  above.  At 
the  completion  of  this  stop  the  brake  heat  sink  was  allowed  to  cool  to 
150°F  before  beginning  the  next  stop.  This  cooling  was  accomplished 
either  by  forced  convection  using  external  air  blowers  (shown  as 
rectangular  ducts  above  the  tire  in  Figure  20)  or  by  natural  convec¬ 
tion.  The  second  test  was  accomplished  to  allow  calculation  of 
convection-radiation  heat  transfer  coefficients  at  various  brake 
rotation  speeds  without  the  brake  being  actuated  during  the  rotation. 
This  calculation  procedure  is  discussed  later.  Thus,  the  procedure  was 
to  complete  a  brake  stop  as  described  above,  but  to  rotate  the  brake 
using  the  dynamometer  at  a  constant  speed  during  the  cooling  process. 
This  was  accomplished  with  the  blowers  as  previously  described  and  also 
without  the  blowers.  In  each  type  of  test  the  procedure  was  to  bring 
the  flywheel  to  the  desired  speed,  shut  off  the  motor,  bring  the  tire 
into  contact  with  the  flywheel,  apply  a  constant  pressure  to  the  brake 
until  the  flywheel  stopped  and  then  to  bring  the  test  wheel  away  from 
the  flywheel.  Table  6  summarizes  all  tests  conducted  as  described 


above.  It  should  be  noted  that  all  tests,  except  the 


3.1  3  stops  using  blowers  after  stators  peak 

3.2  3  stops  with  blowers  off  continuously 

4.  100%  K.E.,  182  mph,  brake  pressure  =  425  psi 

4.1  3  stops  using  blowers  after  stators  peak 

4.2  3  stops  with  blowers  off  continuously 

5.  100%  I.E.,  182  mph,  brake  pressure  =  425  psi 

Make  complete  stop  but  immediately  following  the  stop 
(no  cooling  time) ,  let  the  brake  cool  by  landing  the 
tire  at  about  1,000  lbs  normal  load,  maintaining  speed  in 
table  below  until  the  stators  reach  approximately  200°F 
(no  brake  pressure) . 


TABLE  6-continued 


Stop 

Number 

Velocity  at  which 
to  keep  flywheel 
during  brake  cooling 
(mph) 

Blowers  turned 
on  during  const 
speed  cooldown 

5.1 

50 

Yes 

5.2 

50 

No 

5.3 

150 

Yes 

5.4 

150 

No 

5.5 

100 

Yes 

5.6 

100 

No 

60%  K.E.,  141  mph,  brake  pressure  =  270  psi 

6.1  3  stops  using  blowers  after  stators  peak 

6.2  3  stops  with  blowers  off  continuously 

60%  K.E.,  141  mph,  brake  pressure  =  400  psi 

7.1  3  stops  using  blowers  after  stators  peak 

7.2  3  stops  with  blowers  off  continuously 

100%  K.E. ,  182  mph,  brake  pressure  =  270  psi 

8.1  3  stops  using  blowers  after  stators  peak 

8.2  3  stops  with  blowers  off  continuously 

100%  K.E. ,  182  mph,  brake  pressure  =  400  psi 

9.1  3  stops  using  blowers  after  stators  peak 

9.2  3  stops  with  blowers  off  continuously 

100%  K.E.,  182  mph,  brake  pressure  =  400  psi 

Make  complete  stop  but  immediately  following  the  stop 
(no  cooling  time) ,  let  the  brake  cool  by  landing  the 
tire  at  about  1,000  lbs  normal  load  maintaining  speed  in 
table  below  until  the  stators  reach  approximately  200°F 
(no  brake  pressure) . 


TABLE  6-continued 


Stop 

Number 

Velocity  at  which 
to  keep  flywheel 
during  brake  cooling 
(mph) 

Blowers  turned 
on  during  const, 
speed  cooldown 

10.1 

50 

Yes 

10.2 

50 

No 

10.3 

150 

Yes 

10.4 

150 

No 

10.5 

100 

Yes 

10.6 

100 

No 

assisted  cooling  tests,  were  repeated  three  times  to  insure  adequate 
repeatability  of  both  the  test  conditions  and  the  instrumented 
parameters. 

Instrumentation 

The  purpose  of  this  experimental  portion  of  the  investigation  was 
to  determine  temperatures  within  the  brake  as  a  function  of  time  for 
correlation  with  the  theoretical  developments.  To  accomplish  this,  14 
thermocouples  were  installed  in  the  beryllium  brake  and  24  in  the 
carbon  brake.  Based  on  the  temperature  measurement  analysis  in 
Chapter  III,  type  K  thermocouples  with  a  stainless  steel  sheath  and  a 
grounded  protected  junction  were  used.  The  location  of  the  thermo¬ 
couples  in  the  beryllium  is  shown  in  Figure  29.  The  dashed  line  indi¬ 
cates  the  direction  of  the  thermocouple  leads  and  the  darkened  circle 
indicates  the  location  of  the  thermocouple  junction.  Figure  30  shows 
the  location  of  thermocouples  for  the  carbon  heat  sink.  The  placement 
of  thermocouples,  as  well  as  the  number  of  thermocouples  in  the  carbon, 


is  somewhat  different  than  that  in  the  beryllium.  This  is  primarily 
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Figure  29.  Beryllium  Heat  Sink  Thermocouple  Locations 
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due  to  the  difficulty  of  drilling  holes  in  the  beryllium.  The  carbon 
composite  used  appeared  to  be  of  the  consistency  similar  to  ordinary 
pencil  lead.  It  was  readily  machinable  in  a  lathe  or  drill  press.  The 
excellent  machinability  of  the  carbon  facilitated  drilling  holes  at 
angles  (No.  21,  23,  25  on  Figure  30)  with  precise  depth  control  and  to 
controlled  short  depths  (No.  16-19)  and  long  depths  (No.  2-4,  9-12)  as 
well.  On  the  other  hand,  the  holes  in  the  beryllium  could  not  be 
drilled  or  controlled  nearly  as  well.  The  beryllium  material,  when 
machined,  emits  a  toxic  beryllium  dust  that  can  accumulate  within  the 
human  body  and  can  ultimately  cause  death.  Thus,  machining  of  beryllium 
must  be  conducted  within  a  special  clean  room  and  all  operators  must  be 
protected  such  that  beryllium  dust  cannot  come  in  contact  with  any  part 
of  the  body.  This  includes  the  wearing  of  gloves,  hats,  face  masks, 
goggles,  and  a  special  mask  for  filtering  the  breathing  air. 

These  machine  operator  hinderances,  plus  the  fact  that  the  holes 
were  chemically  etched  into  the  beryllium,  resulted  in  limited  control 
compared  to  that  in  the  carbon.  Therefore,  it  was  decided  to  use  only 
shallow  holes  in  the  beryllium  and  not  complicate  the  procedures  with 
critical  angles  or  depths.  The  specific  holes  in  Rotor  3  and  Stators  1 
and  2,  including  the  angled  holes,  all  of  which  are  shown  in  Figure  30, 
are  more  explicitly  defined  in  Figure  31  for  Rotor  3  and  in  Figure  32 
for  Stators  1  and  2.  All  thermocouples  other  than  those  in  the  heat 
sinks  are  shown  in  Figure  33.  Due  to  the  lower  temperature  gradients 
in  the  brake  parts  of  Figure  33,  glass  insulation-type  thermocouples 
were  selected  for  those  areas.  All  thermocouple  outputs  were  recorded 
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Figure  31.  Carbon  Rotor  3  Thermocouple  Locations 
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Figure  32.  Carbon  Stator  1  and  2  Thermocouple  Locations 
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Overall  Thermocouple  Locations 


on  digital  recorders.  To  increase  the  sampling  rate,  two  recorders 
were  used  during  each  test.  A  Digitrend  220  Doric  recorder  and  a  Kaye 
Model  8000  recorder  were  used  for  all  tests.  For  the  carbon  brake 
tests.  Thermocouples  2-13  were  recorded  on  the  Kaye  recorder  while 
Thermocouples  14-25  were  recorded  on  the  Doric.  For  the  beryllium. 
Thermocouples  2-8  were  recorded  on  the  Kaye  and  9-13  on  the  Doric. 

This  arrangement  increased  the  sampling  rate  to  12  samples  per  minute. 
To  insure  that  the  start  of  a  brake  stop  could  be  synchronized  with  the 
temperature  information,  at  least  one  channel  of  each  temperature 
recorder  was,  in  addition  to  being  recorded  on  the  temperature 
recorder,  recorded  on  magnetic  tape  along  with  brake  pressure.  This 
permitted  comparison  of  the  two  temperatures  from  the  same  channel  and 
provided  the  desired  time  reference.  Rotor  thermocouple  output  was 
channeled  through  a  commutator  (Figure  28)  in  order  to  record  the 
outputs.  The  commutator  longitudinal  axis  is  in  line  with  the  mandrel 
axis. 

Also  shown  in  Figure  28  is  the  instrument  used  to  measure  the  test 
wheel  speed.  A  light  source  is  located  just  to  the  left  of  the  com¬ 
mutator  that,  during  tire  rotation,  is  continually  interrupted  by  a 
series  of  small  bolts  at  the  outer  diameter  of  the  small  plate  in 
Figure  28.  The  corresponding  pulse  signals  are  fed  into  a  digital-to- 
analog  converter  and  the  resulting  analog  signal  is  recorded  on  mag¬ 
netic  tape.  Thus,  it  is  a  relatively  simple  matter  to  accurately 
determine  the  test  wheel  speed. 


Flywheel  speed  is  measured  similarly  and  is  also  recorded  on 
magnetic  tape.  In  this  case  the  speed  at  the  flywheel  outer  diameter 
is  found  in  feet  per  second. 

The  torque  developed  during  a  brake  stop  is  of  primary  interest  in 
determining  the  energy  absorbed  by  the  brake  as  a  function  of  time. 

The  torque  was  obtained  by  mounting  a  temperature  compensated  strain 
gauge  circuit  to  the  mandrel  and  calibrating  it  using  a  torque  arm 
applied  to  the  mandrel  when  the  wheel-brake  assembly  was  removed.  This 
signal  was  also  recorded  on  tape  for  subsequent  analysis. 

Brake  pressure  was  obtained  using  a  pressure  transducer  upstream 
of  the  hydraulic  inlet  line  as  depicted  in  Figure  27.  The  pressure 
measured  in  this  way  is  the  pressure  applied  to  the  brake  and  includes 
not  only  the  pressure  forcing  the  discs  together  but  also  that  pressure 
required  to  overcome  the  force  of  the  return  springs. 

Normal  load  was  found  by  measuring  the  pressure  within  the  pneu¬ 
matic  loading  system  and  multiplying  it  by  the  applicable  area.  This 
was  recorded  on  tape  in  pounds  force. 

The  deflection  of  the  tire  was  also  monitored  and  recorded.  A 
linear  velocity-displacement  transducer  was  used  in  this  case  and 
similarly  recorded  on  magnetic  tape. 

Each  of  the  above  mentioned  parameters,  including  the  two  tempera¬ 
ture  channels,  was  recorded  on  a  Bell  and  Howell  VR-3700B  Oatatape 
unit.  The  specific  units  used  are  shown  in  Figure  34.  Also  shown  in 
Figure  34  is  the  signal  conditioning  equipment  which  contains  the 
capability  of  recording  an  audible  signal  describing  each  brake  stop  as 

to  the  test  number  and  test  conditions.  This  was  accomplished  for  all 
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tests.  The  unit  on  the  far  right  of  Figure  34  is  the  patchboard  where 
all  signals  from  the  dynamometer  station  are  patched  into  the  signal 
conditioning  equipment.  Once  the  data  is  recorded,  or  even  while  it  is 
being  recorded,  it  may  also  be  channeled  into  a  brush  recorder  so  that 
a  hard  copy  output  of  each  parameter  can  be  displayed  as  a  function  of 
time.  Figure  35  shows  the  brush  recorders  used.  Sample  output  from 
these  recorders  is  shown  and  discussed  in  the  test  results  section. 

Results 

A  total  of  60  tests  was  conducted.  Of  these,  48  were  conducted 
for  comparison  with  theoretical  models.  The  remaining  12  were  run 
specifically  to  determine  the  convective-radiative  heat  transfer 
coefficients  for  use  in  the  theoretical  models.  Since  the  theoretical 
models  are  described  in  Chapter  V,  the  discussion  of  these  12  tests 
will  be  deferred  until  Chapter  V  where  boundary  conditions  are 
developed. 

Recorded  Parameters  Versus  Time 

Ten  tests  were  selected  for  detailed  study  and  temperature  plotting. 
Table  7  shows  the  particular  ten  stops  selected.  The  tests  can  be 
identified  by  the  test  number  which  corresponds  to  the  paragraph 
numbers  shown  in  Table  6,  the  test  plan.  The  third  digit  in  the  test 
number  is  the  particular  stop  of  the  three  conducted  at  that  paragraph. 
The  ten  stops  were  selected  to  encompass  a  range  of  kinetic  energies 
and  decelerations.  In  addition,  six  tests  were  selected  with  cooling 
blowers  on  and  four  tests  were  selected  with  cooling  blowers  off. 

Figure  36  shows  typical  data  recorded  during  each  test  as  a  function  of 


TABLE  7 


60%  of  Normal  K.E 


time.  The  graph  shown  is  for  Test  1.1.3,  the  first  listed  on  Table  7, 
but  is  typical  of  all  the  tests.  The  scale  for  each  recorded  parameter 
is  given  under  that  parameter.  The  time  divisions  are  shown  as  one- 
second  intervals  and  are  located  between  the  pressure  and  torque 
traces.  The  complete  brake  stop  can  be  followed  by  considering 
Figure  36.  The  first  thing  to  note  is  that  the  flywheel  has  been 
brought  up  to  the  required  speed.  The  tire  is  moved  into  contact  with 
the  flywheel  at  this  point.  As  noted  in  Figure  36,  the  test  wheel 
speed,  tire  deflection,  and  normal  load  all  rise  simultaneously.  About 
three  and  one-half  seconds  later  the  brake  pressure  is  applied.  The 
torque  then  increases  immediately  and  the  flywheel  begins  to  slow  down. 
At  this  point  the  brake  is  starting  to  heat  up.  This  is  evident  from 
the  rise  noted  in  the  two  temperature  traces. 

Approximately  one  minute  later,  the  flywheel  is  stopped  and  the 
pressure  and  load  are  returned  to  zero.  The  first  six  traces  are  now 
zero  and  the  cooling  period  begins.  Another  stop  is  not  made  until  the 
brake  has  cooled  to  150°F.  Table  7  lists  data  that  have  been  extracted 
from  the  traces  of  Figure  36.  The  average  deceleration  was  calculated 
by  dividing  the  initial  speed  by  the  stop  time.  The  average  torque  was 
obtained  by  approximating  the  average  torque  during  each  second,  adding 
these  values  for  the  entire  stop  and  dividing  by  the  stop  time  of  the 
stop.  Appendix  E  contains  the  data  recorded  during  each  test,  as  shown 
in  Figure  36,  for  each  of  the  ten  stops  of  Table  7.  Unfortunately,  due 
to  instrumentation  difficulties,  two  stops  do  not  have  test  wheel  speed 
shown  but  do  contain  all  of  the  other  more  important  parameters. 
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Temperature  Versus  Time 

The  temperature  in  Rotor  2,  Stator  2,  and  the  torque  tube  is  shown 
as  a  function  of  time  for  Stop  1.1-3  in  Figure  37.  The  exact  thermo¬ 
couple  locations  are  given  in  Figure  30.  The  left  side  of  Figure  37 
has  an  expanded  time  scale  to  allow  observation  of  the  initial  high 
temperature  gradients.  The  curves  on  the  right  show  the  cooling  that 
takes  place  over  the  next  70  minutes.  The  point  at  which  the  brake 
stopped  the  flywheel  is  seen  to  be  at  58  seconds  and  slightly  after  the 
peak  temperature  occurs  in  the  stator  and  rotor.  The  heat  from  the 
stators  and  rotors  flows  into  the  surrounding  brake  hardware,  including 
the  torque  tube.  The  peak  temperature  in  the  torque  tube  occurs  later 
than  the  peak  temperatures  of  the  heat  sink.  Appendix  F  contains 
curves  of  this  type  for  each  test  in  Table  7. 

Temperature  Versus  Axial  Position 

The  way  in  which  the  temperature  varies  in  the  axial  direction 
during  a  stop  is  depicted  in  Figure  38.  Here  the  temperature  plotted 
is  from  the  thermocouples  at  the  bottom  of  the  graph.  They  are  located 
in  the  disc  shown  above  the  thermocouple  number  on  Figure  38,  and  the 
exact  location  can  again  be  found  on  Figure  30.  The  data  points  were 
connected  by  straight  lines  since  the  information  between  the  data 
points  does  not  indicate  the  brake  temperatures  in  that  area.  As  can 
be  seen  in  the  graph,  the  rotors  tend  to  be  hotter  than  the  stators  at 
early  time,  but  this  difference  decreases  as  time  increases.  Also,  it 
is  evident  that  the  temperature  is  less  near  the  left  end  or  brake 
housing.  Appendix  G  gives  these  axial  position  plots  for  all  ten  of 
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the  stops  of  Table  7.  The  beryllium  curves  are  plotted  somewhat 
differently  than  the  carbon  since  the  thermocouples  in  the  beryllium 
were  not  in  a  straight  axial  line.  Due  to  this,  at  each  time  points 
are  plotted,  the  stator  points  are  connected  to  stator  points  and  the 
rotor  points  are  connected  to  rotor  points.  Figure  29  gives  the  exact 
beryllium  thermocouple  locations. 

Temperature  Versus  Radial  Position 

Figure  39  shows  a  plot  of  temperature  versus  Lime  for  the  four 
thermocouples  located  in  a  radial  line  in  the  backing  plate  of  the 
carbon  heat  sink.  The  exact  location  of  the  thermocouples  is  shown  in 
Figure  30.  The  inner  and  outer  thermocouple  locations  appear  to  be 
cooler  than  the  other  two  at  early  times,  but  all  four  cooling  curves 
coincide  at  later  times.  The  radial  position  plots  for  each  of  the 
five  carbon  brake  tests  identified  in  Table  7  are  given  in  Appendix  H. 

To  more  clearly  understand  the  temperature  pattern  in  the  radial 
direction  for  early  times.  Figure  40  was  plotted.  This  figure  shows 
temperature  versus  radial  position  at  various  times  up  until  the  stop 
time.  Here  the  particular  radial  temperature  is  seen  at  each  time 
increment.  The  thermocouple  position  is  shown  on  the  graph,  as  is  the 
wearing  surface  location.  At  early  time  the  middle  two  thermocouples 
are  higher  in  temperature  than  the  other  two.  The  curve  is  fairly  flat 
by  the  time  the  stop  time  is  reached.  Also,  the  peak  temperature 
occurs  earlier  than  at  the  stop  time.  Appendix  I  contains  similar 
temperature  versus  radial  position  plots  for  the  five  carbon  brake 
stops  of  Table  7. 
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In  order  to  observe  the  radial  temperature  variation  in  Stator  1, 
Figure  41  was  plotted.  The  two  thermocouples  plotted  are  both  at  the 
axial  center  of  Stator  1  but  are  at  different  radii,  4.4.  in.  and 
5.5  in.  The  thermocouple  at  the  largest  radius  has  a  higher  gradient 
but  reaches  about  the  same  peak  temperature  as  the  other  thermocouple. 
Figure  42  shows  the  same  characteristics;  however,  the  test  was  con¬ 
ducted  with  a  higher  energy  input  and  deceleration.  Figure  41  shows  a 
peak  temperature  of  about  900°F  while  Figure  42  shows  a  peak  of  about 
1500°F . 

Figure  43  is  a  plot  of  temperature  versus  time  for  two  thermo¬ 
couples  at  the  same  radial  position  in  Stator  1.  The  difference 
between  the  two  is  that  one  is  located  in  the  axial  center  of  the 
stator  and  the  other  is  very  close  to  the  wearing  surface  (Figure  30). 
The  temperature  measured  close  to  the  wearing  surface  is  slightly 
higher  than  that  of  the  other.  The  difference  would  probably  have 
been  higher  if  the  radius  were  increased  so  that  the  angled  thermo¬ 
couple  would  not  be  at  the  radial  edge  of  the  rubbing  surface.  At  a 
higher  energy  and  deceleration,  this  temperature  difference  is  higher 
(see  Figure  44) .  It  is  also  seen  that  during  cooling  the  curves  are 
coincident  in  the  first  case  and  almost  so  in  the  second. 

Brake  Hardware  Temperature  Versus  Time 

Four  stops  were  selected  to  show  the  temperatures  reached  in  the 
tire-wheel  assembly  other  than  in  the  heat  sink.  Figure  45  shows  such 
a  plot  for  the  carbon  brake  at  a  low  deceleration  (3.58  ft/sec'-)  and 
at  60%  of  normal  kinetic  energy.  The  thermocouple  numbers  and  loca¬ 
tions  are  specified  on  the  graph.  Figure  33  gives  a  pictorial  view  of 
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42.  Stator  1  Radial  Temperature  (Test  4.2-3) 
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Figure  45.  Temperature  vs  Time-Tire/Wheel  Assembly  (Test  1.1-3) 


these  particular  thermocouple  locations.  The  wheel  drive  lug  attains 
the  highest  temperature  while  the  hydraulic  fluid  is  the  lowest.  The 
air  temperature  in  the  tire  is  also  low  and  only  slightly  above  the 
hydraulic  fluid  temperature.  For  a  higher  deceleration  and  kinetic 
energy.  Figure  46  shows  that  the  air  temperature  in  the  tire  remains 
low  but  that  the  hydraulic  fluid  temperature  rises  significantly  over 
the  test  previously  described.  The  highest  temperature  is  again 
attained  in  the  wheel  drive  lug.  Another  difference  in  the  two  graphs 
is  the  relative  position  of  the  temperature  of  the  wheel  rim.  For  the 
low  energy  and  low  deceleration,  the  wheel  rim  temperature  is  the 
second  highest  of  the  six  curves;  in  the  other  case,  the  wheel  rim  is 
the  fourth  highest  in  temperature. 

Figures  47  and  48  show  similar  plots  for  the  beryllium  brake.  The 
large  jump  in  the  hydraulic  fluid  temperature  from  the  low  to  high 
deceleration  can  again  be  noticed.  On  the  high  deceleration  stop 
(Figure  47),  the  temperature  of  the  air  in  the  tire  did  not  rise 
significantly,  while  on  the  low  deceleration  stop  the  air  temperature 
did  rise  significantly.  Although  cooling  blowers  were  used  in  the  low 
deceleration  case,  the  air  temperature  increase  for  this  case  is 
attributed  to  the  higher  energy  input  (100%  vs  60%) .  For  the  beryllium 
brake  it  is  noted  that  although  the  wheel  rim  temperature  is  close  to 
the  highest  temperature  in  Figure  48,  it  is  the  next  to  lowest  in 
Figure  47.  The  apparent  discrepancy  in  the  results  is  probably  due  to 
two  factors.  The  first  is  that  the  combination  of  low  deceleration, 
low  energy  input  found  in  the  carbon  brake  stop  of  Figure  45  is  not 
found  in  either  of  the  two  beryllium  brake  stops.  The  conditions  of 
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Figure  48.  Temperature  vs  Time-Tire/Wheel  Assembly  (Test  8.1-2) 


Figure  46  were  also  not  duplicated  by  either  of  the  two  beryllium  brake 
stops  shown.  The  second  reason  is  that  the  carbon  brake  was  equipped 
with  heat  shields  that  were  placed  circumferentially  between  the  wheel 
rim  inside  diameter  and  the  heat  sink  outer  diameter.  These  were 
inserted  to  keep  the  wheel  rim  from  obtaining  a  high  temperature  when 
the  carbon  heat  sink  obtained  a  high  temperature.  Due  to  the  low 
temperatures  in  the  beryllium  heat  sink,  the  heat  shields  were  not  used 
for  the  beryllium  brake  tests. 
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CHAPTER  V 


THEORETICAL  DEVELOPMENT 

Introduction 

The  theoretical  models  described  in  this  chapter  are  all  concerned 
with  the  multidisc  brake  except  the  last  which  deals  with  a  single 
rotor  caliper  type  brake.  The  analysis  is  most  directed  toward  the 
multidisc  brake  because  of  the  extensive  testing  that  was  accomplished 
on  the  multidisc  brakes  as  reported  in  Chapter  IV.  Within  the  multidisc 
brake  analysis,  several  models  employing  different  assumptions  were 
developed.  These  allowed  a  gradual  buildup  to  the  most  extensive 
model.  In  addition,  the  validity  of  less  complicated  models  is  shown 
since  all  are  compared  to  experiment.  This  approach  has  seldom,  if 
ever,  been  reported  for  the  brake  application  before. 

During  the  course  of  the  theoretical  development,  the  term  "R/S 
Pair"  is  used.  This  refers  to  a  commonly  made  assumption  within  multi¬ 
disc  brake  analysis.  The  assumption  is  that,  if  the  area  between  the 
radial  centerlines  of  an  adjacent  rotor  and  stator  is  analyzed,  all 
other  such  areas  within  the  brake  are  the  same  and  need  not  be  analyzed. 
This  greatly  simplifies  the  analysis  but  has  the  drawback  of  ignoring 
the  axial  temperature  gradient  within  the  heat  sink  that  has  been  shown 
experimentally  to  exist  (see  Figure  38). 

In  several  of  the  multidisc  analyses,  a  boundary  condition 
resulting  from  convection  and  radiation  at  the  brake,  wheel,  and  tire 
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exposed  boundaries  is  utilized.  The  actual  values  of  the  convective/ 
radiative  heat  transfer  coefficients  have  been  calculated  from  the 
experimental  results.  The  values  calculated  and  the  calculation  pro¬ 
cedure  are  described  in  Appendix  J. 

Numerical  Solution  Technique  Selection 

The  equations  to  be  developed  in  subsequent  sections  of  this 
chapter  prove  to  be  nonlinear  due  to  the  inclusion  of  variable  thermal 
conductivity  and  heat  capacity.  Also,  the  boundary  condition  at  the 
rotor/stator  interface  contains  a  work  term  that  varies  with  radius  as 
well  as  time.  Due  to  these  facts,  a  numerical  solution  technique  was 
necessary.  The  objective  of  this  section  of  the  research  is  to  select 
a  simple,  one-dimensional,  theoretical  model  simulating  an  aircraft 
brake  that  will  enable  comparison  of  an  exact,  finite  difference  and 
finite  element  solution  in  the  areas  of  accuracy,  computer  storage  and 
computer  time  required.  It  should  be  noted  that  the  results  indicated 
herein  are  equally  applicable  to  any  disc  type  brake,  aircraft,  truck 
or  car. 

Theoret/ ;al  Model 

The  problem  selected  simulates  a  rotor  or  stator  disc  in  a  brake 
in  that  the  entire  piece  of  material  is  initially  at  a  constant  tem¬ 
perature,  one  edge  is  adiabatic  which  is  the  disc  centerline  and  at 
time  equal  zero,  heat  crosses  the  face  opposite  the  adiabatic  face.  In 
the  brake  disc,  heat  is  added  due  to  the  friction  between  the  discs. 

In  the  model  selected,  heat  is  taken  from  the  material  by  keeping  the 
face  at  a  constant  temperature  which  is  lower  than  the  initial 
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temperature.  The  model  selected  is  shown  in  Figure  49.  The  face  at 


ADIABATIC 

DISC 

CENTERLINE 


Figure  49.  Theoretical  Model 


x  =  L  is  insulated  and  the  entire  piece  of  material  is  initially  at 
temperature  Tq.  The  face  at  x  =  0  is  suddenly  changed  to  temperature 
Tg.  The  governing  equation  using  constant  properties  and  one  spatial 
dimension  is 


a2  t  _  i  ar 
a  x2  =  °  3t 


(29) 


where  a  is  the  thermal  diffusivity. 


Boundary  and  initial  conditions  are: 


ST 

T(x,0)  =  V  it (L>t)  =  °»  T(°.o  =  Ts 


(30) 


Define  nondimensional  quantities  as  follows: 
T  -  T 


T  -  T  ’ 

•o  s 


*'  -  r- c'  - 


a  t 


(31) 


Then  equation  (29)  becomes: 


3  2  T'  3T ' 


3  x 


,2  3t ' 


(32) 


and  the  initial  and  boundary  conditions  are: 

T’(x\0)  =  1,  T’  (0,t')  =  0,  r,  (l,t')  =  0  (33) 

where  the  x'  subscript  denotes  differentiation  with  respect  to  x’ . 

Analysis 

Exact  Solution.  Assuming  that  the  solution  to  the  above  equation 
consists  of  the  product  of  a  function  of  x'  and  a  function  of  t’,  two 
ordinary  differential  equations  can  be  obtained,  the  solutions  of  which 
can  be  combined  to  give  the  solution  to  the  original  equation.  With 
the  initial  and  boundary  conditions  applied,  the  exact  soJution^  is: 


(34) 


oo 

(x,0)  =  (2/ it)  ^  ^  exp  I  -  /n  +  ^-j2  "2  tl 
n  =  0  L  J 

coS  [(n  +  i).(l-x)] 


where  the  primes  have  been  omitted  for  convenience. 

Finite  Difference  Solution.  Equation  (32)  written  in  finite  difference 


form  is 


T  =  T^I)  1  -  — — ]  +  —■ — -  [T(I+1)  +  T(I-l)]t  (35) 


(Ax) 2  I  (Ax)2 


This  equation  can  be  programmed  directly  using  an  iterative  scheme 
within  the  computer.  It  can  also  be  written  in  matrix  form. 


Tt+At  =  A 


In  this  case,  operations  within  the  computer  would  involve  matrices. 

For  comparison  of  storage  and  time  requirements  both  methods  were  used 
and  are  referred  to  as  "Finite  Difference  Equation"  and  "Finite 
Difference  Matrix,"  respectively. 

For  stability  and  compliance  with  the  second  law  of  thermodynamics, 
the  quantity  multiplying  Tt(I)  must  be  greater  than  zero.  Since  the 
number  of  modes  was  taken  to  be  five,  the  nondimensional  Ax  value 
becomes  .25  and  the  following  equation  results: 


-«t.-  -  — ■ 


This  equation  yields  a  At  of  .03125.  This  is  a  conservative  value  and 
for  this  problem  a  value  of  .031967  was  found  by  trial  and  error  to  be 
the  maximum  value  allowable  for  a  stable  solution. 

Finite  Element  Solution.  The  first  step  in  the  Finite  Element  Solution^1* 
is  to  find  the  variational  statement  from  the  second  order  differential 
equation  given  earlier.  This  statement  is 


(38) 


This  integral  is  then  manipulated  to  obtain  equations  to  solve  for 
temperature.  This  manipulation  leads  to  the  matrix  equation 


K  T  +  C  T  =  0 


(39) 


where  the  bars  indicate  matrices  and  the  dot  indicates  differentiation 
with  respect  to  time.  K  is  the  global  conduction  matrix  and  C  is  the 
global  capacitance  matrix.  The  solution  is  advanced  in  time  by  using 
the  Euler  numerical  solution  technique: 

-t+At  _  +  At  (40) 


Evaluating  Equation  (39)  at  time  t  and  substituting  into  Equation  40 
yields : 

C  =  B  T*  (41) 


where 


B  =  C  -  AtK 


.  1^ 
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Solving  for 


zxr- 


"■ ~  ** 1  jp  *> 


^t+At 


yields 


T*  +  iC  -  C  -1  B  T4 


(42) 


This  equation  can  now  be  solved  at  each  time  step  to  obtain  the  solu¬ 
tion  to  the  differential  equation  as  long  as  initial  conditions, 
boundary  conditions,  C,  and  B  are  known.  The  number  of  nodes  in  this 
one-dimensional  problem  was  again  five,  which  results  in  four  elements. 

For  comparison  with  the  finite  difference  technique,  a  method  of 
solving  the  finite  element  equations  without  using  matrices  was  desired 
in  addition  to  the  above  described  procedure.  The  method  employed  was 
to  write  Equation  (39)  as  a  series  of  equations  and  then  solve  them 
algebraically  for  the  required  time  derivatives.  This  practice  is 
cumbersome  and  not  recommended,  but  it  does  allow  a  direct  comparison 
with  the  other  methods.  The  two  finite  element  computer  solutions  are 
referred  to  as  "Finite  Element  Matrix"  and  "Finite  Element  Equation," 
respectively.  The  one-dimensional  finite  element  critical  time  step^-3 
was  found  to  be: 


1  - 


6At 

(.25)2 


>  0 


(«) 


This  gives  the  conservative  value  of  .01042.  Trial  and  error  produced 
the  value  of  .0115  for  the  maximum  value  allowable  for  a  stable 
solution. 


119 


Comparison  of  Finite  Element  Method  to 
Finite  Difference  Method 

The  exact,  finite  difference,  and  finite  element  solutions  were 

numerically  evaluated  for  various  time  increments,  elapsed  time  values, 

and  at  five  spatial  locations.  Figures  50  and  51  show  plots  of  T  - 

T  versus  position,  where  T  is  the  temperature  predicted  by  the 

6X3C  L 

finite  element  or  finite  difference  technique  and  T  _  is  the  tern- 

exact 

perature  as  given  by  the  exact  analysis  at  the  same  time  and  position 

as  the  approximate  analyses.  As  shown  in  Figure  50  at  time  equal  .5, 

the  error  in  using  the  finite  difference  techniques  increases  with 

increasing  x  and  with  increasing  time  increment.  Recall  that  T, 

T  ,  x,  and  time  are  nondimensional  as  defined  in  the  Theoretical 

exact 

Model  Section.  It  is  interesting  to  note  that  results  with  the  finite 
difference  technique  show  increasing  error  with  increasing  x  in 
Figure  51  also,  but  show  just  the  opposite  effect  of  that  in  Figure  50 
with  respect  to  the  time  increment;  the  error  decreases  with  increasing 
time  increment.  The  difference  between  the  two  figures  other  than  the 
results  is  that  Figure  50  is  plotted  for  time  =  .5  while  Figure  51  is 
plotted  for  time  =  1.0. 

Figure  52  shows,  for  time  =  .5,  the  absolute  value  T  -  Texact 
versus  x  for  the  finite  element  technique.  Error  increases  with 
increasing  x,  but  no  trend  is  evident  of  error  change  with  time  incre¬ 
ment.  Figure  53  gives  the  same  information  as  Figure  52  but  at  time  = 
1.0.  In  this  case,  error  increases  with  increasing  x  and  increasing 
time  increment.  This  time  increment  trend  from  time  =  .5  to  time  =  1.0 
is  opposite  to  that  found  using  the  finite  difference  technique. 

Figure  54,  for  time  =  .5,  compares  finite  difference  and  finite  element 
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Figure  52.  Temperature  Difference  vs  Position 
Finite  Element ,  Time  =  . 5 


Temperature  Difference  vs  Position, 
Finite  Element,  Time  =  1.0 


EXACT 


techniques  for  two  time  increments.  As  can  be  seen,  the  larger  time 
increments  show  larger  error  than  the  smaller  time  increments  and  at  a 
given  time  increment  the  finite  difference  technique  has  greater  error 
than  the  finite  element  technique.  The  same  information  is  shown  in 
Figure  55  for  time  =  1.0.  Here  the  finite  element  technique  error  is 
larger  than  the  finite  difference  error.  For  the  finite  element 
technique  the  larger  time  increment  gives  the  higher  error  while  within 
the  finite  difference  technique,  the  higher  error  is  produced  by  the 
smaller  time  increment.  These  results  do  not  appear  conclusive  and 
must  be  looked  at  in  another  way  to  provide  meaningful  results.  A  plot 
of  T  -  T  ^  versus  time  is  given  in  Figure  56  for  DT  =  .00390625. 

This  shows  a  comparison  of  the  finite  element  and  the  finite  difference 
techniques.  Each  point  represents  an  average  of  four  T  -  values 

at  that.  time.  The  four  points  averaged  are  those  at  x  =  .25,  .50,  .75, 
and  1.0.  At  small  times  the  finite  element  technique  has  much  greater 
error  than  the  finite  difference  technique.  The  finite  element  error 
at  time  =  .03125  is  about  three  and  one-half  times  greater  than  the 
finite  difference  error  at  the  same  time.  The  maximum  error  at  time  = 
.03125  for  the  finite  element  technique  is  4.7%  compared  with  1.1%  for 
the  finite  difference  technique.  These  values  change  to  2.2%  and  1.7% 
respectively  at  time  =  1.0.  Figure  57  is  comparable  to  Figure  56  and 
is  plotted  for  DT  =  .0078125.  Here  the  maximum  error  in  the  finite 
element  technique  at  time  =  .03125  is  7.5%  compared  to  3.8%  for  the 
finite  difference  technique.  At  time  =  1.0  these  values  reduce  to  2.4% 
and  1.7%,  respectively. 
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Figure  57.  Temperature  Difference  vs  Time 
DT  =  .0078125 


In  addition  to  error  values,  the  time  required  to  perform  the 
analyses  was  recorded.  Figure  58  is  a  plot  of  the  time  increment  used 
versus  the  computer  time  used,  given  in  CPU  seconds  on  the  CDC  6600 
computer.  It  is  apparent  that  the  time  required  is  greater  for  the 
finite  element  technique  as  long  as  the  time  increment  is  the  same  and 
the  equation  formulation  within  each  technique  is  the  same.  Note  that 
the  upper  limit  in  DT  is  different  for  the  two  techniques.  The  upper 
limit  in  DT  represents  the  maximum  value  that  can  be  used  and  still 
maintain  a  stable  solution.  Thus,  a  larger  time  increment  can  be  used 
within  the  finite  difference  technique  and  this  means  that  less  com¬ 
puter  time  is  necessary  to  solve  the  problem. 

The  fact  that  the  time  required  is  higher  in  the  finite  element 
case  is  at  least  partially  caused  by  the  fact  that  more  arrays  are 
necessary  in  that  case.  Table  8  compares  the  exact,  finite  difference 
and  finite  element  techniques  in  the  areas  of  storage  requirements  and 
computer  time  necessary.  The  finite  element  method  requires  consider¬ 
ably  more  time  and  storage  than  the  finite  difference  method  when  using 
matrices.  Also,  although  array  requirements  are  the  same  when  using 
the  equation  formulation,  the  finite  element  technique  is  still  more 
time  consuming  than  the  finite  difference  technique.  Note  that 
although  the  exact  method  uses  no  arrays,  it  takes  much  more  time  than 
any  of  the  others.  This  suggests  that  even  if  the  exact  solution  to  a 
problem  exists,  it  could  be  more  economical  to  use  a  numerical  solution. 
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TABLE  8 


EXACT,  FINITE  DIFFERENCE  AND  FINITE  ELEMENT  COMPARISON 


Required 

Arrays 


One-Dimensional  0 


One-Dimensional 


Two-Dimensional  0 


Two-Dimens ional 


Computation 
Time  (CPU-SEC) 

DT  =  .031967 

DT  =  .015625 

DT  =  .0115 

DT  =  .0078125 

DT  =  .00390625 

DT  =  .001953125 


- 1 

Exact  With 

200  Terms 

Finite  Difference 
Equation,  DT  = 

.031967  C 

l! 

Finite  Element 
Equation,  DT  .  = 

.0115  Crit 

Finite  Difference 

Matrix,  DT  = 

crit  1 

.031967 

- II 

0 

1  with 

1  with 

2  with 

5  EL* 

5  EL 

4  EL 

0 

0 

0 

2  at 

(4  x  4) 

2  at 

(4  x  1) 

.811 

.010 

.013 

.010 

.019 

.011 

.013 

.025 

.017 

.015 

.029 

.024 

.027 

.057 

.041 

.050 

.111 

(U  x 

U  -H 
•H  W 
C  ^  «H 
•H  <fl  O 


2  with 


4  EL 


7  at 
(5  x  5) 
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Conclusions 


Based  on  a  comparison  of  finite  difference,  finite  element  and 
exact  solutions  of  the  one-dimensional  model  of  a  brake: 

1.  The  finite  element  method  gives  about  three  times  more  error 
at  small  times  as  does  the  finite  difference  method.  At  larger  times 
the  absolute  value  of  the  error  is  about  the  same  for  the  two  methods. 

2.  The  computer  time  required  at  the  same  time  increment  is 
greater  for  the  finite  element  method  than  for  the  finite  difference 
method. 

3.  The  critical  time  increment  as  far  as  stability  is  concerned 
is  almost  three  times  greater  in  the  finite  difference  method  as  com¬ 
pared  to  the  finite  element  method. 

4.  Storage  requirements  are  more  extensive  with  the  finite 
element  method. 

5.  The  exact  solution  requires  much  more  time  than  either  of  the 
numerical  solutions. 

6.  Whether  or  not  the  finite  difference  method  has  less  error  at 
a  given  time  and  time  increment  than  the  finite  element  method  depends 
on  both  the  particular  time  and  particular  time  increment. 

Discussion 

The  conclusions  of  above  pointed  toward  the  use  of  the  finite 
difference  technique.  Pertinent  literature  on  the  subject  was  also 
investigated.  Most  authors  seem  to  choose  a  technique  with  no  reasons 
given  for  the  selection.  Huebner,36  among  others,  has  pointed  out  that 
the  finite  element  method  is  superior  in  treating  irregular  geometries 


and  that  it  can  be  more  accurate  than  finite  difference  in  steady  state 
heat  transfer  problems.  Emery  and  Carson37  studied  the  finite  element 
method  as  used  in  heat  transfer  problems.  They  concluded: 

1.  An  oscillation  at  early  time  exists  in  transient  problems 
within  the  finite  element  method.  This  is  not  a  numerical  instability 
but  is  inherent  within  the  technique  itself. 

2.  More  severe  numerical  oscillation  limits  exist  in  transient 
problems  in  the  finite  element  technique. 

3.  Variable  property  steady  state  problems  usually  require  large 
amounts  of  computer  execution  time  when  the  finite  element  method  is 
used,  an  order  of  magnitude  more  than  the  finite  difference  technique. 

4.  Due  to  item  3,  variable  property  transient  problems  will 
compound  the  time  problem  and  thus  should  not  be  treated  by  the  finite 
element  technique. 

5.  Problems  such  as  3  and  4  of  above  are  not  only  slower  using 
the  finite  element  method  but  also  require  an  order  of  magnitude  more 
computer  storage. 

Due  to  the  results  of  this  section  of  the  research,  all  models  employ 
the  finite  difference  technique  for  solution. 

One-Dimensional  Radial 

Interior  Node  Equation  Development 

The  model  selected  for  this  radial  analysis  is  shown  in  Figure  59. 
The  nodes  and  typical  elements  of  mass  are  shown  as  well  as  applicable 
dimensions.  The  equation  for  the  interior  nodes  was  derived  from  the 
first  law  of  thermodynamics: 
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Figure  59.  One  Dimensional  Radial  Model 


dE 

qin  qout  e  dt 


(44) 


where  q.  is  the  heat  transfer  rate  into  the  elemental  mass,  q  is 
in  ’  ^out 

the  heat  transfer  rate  out  of  the  elemental  mass,  W  is  the  rate  work 

e 

crosses  the  elemental  mass  boundary,  E  is  the  total  energy  of  the 
elemental  mass  and  t  is  the  time.  If  the  kinetic  and  potential  energy 
changes  with  respect  to  time  are  discarded,  the  total  energy  becomes 
only  the  internal  energy,  U  or  the  specific  internal  energy,  u,  times 
the  mass,  pV.  Since  the  specific  internal  energy  can  be  considered  a 
function  of  temperature,  T,  and  specific  volume,  v,  the  differential  of 
u  can  be  expressed  as: 

d"  -  (|v)t  +  (!?)v  dT  <45> 

where  the  subscripts  indicate  evaluation  at  constant  temperature  and 
constant  specific  volume,  respectively.  For  a  solid,  dv  is  approxi¬ 
mately  zero  so  the  first  term  in  Equation  (45)  becomes  zero.  Also  the 
quantity  multiplying  dT  in  Equation  (45)  is  the  definition  of  specific 
heat  at  constant  volume,  C.  Thus,  Equation  (45)  substituted  into 
Equation  (44)  yields: 


qin  -  qout  -  We 


PVC  dF 


An  expression  for  in  terms  of  known  or  measurable  quantities  is 
desired.  Considering  the  bigger  of  the  two  cross-hatched  elements  in 
Figure  59,  an  expression  for  the  work  input  to  that  element  is: 


(47) 
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where  W  is  the  rate  work  crosses  the  element  boundary,  F  is  the  force 
e  e 

pushing  the  rotor  element  past  the  stator  element  (Figure  59) ,  D  is  the 


distance  through  which  the  force,  Fe>  acts  and  t  is  the  time  interval 
through  which  this  force  acts.  The  negative  sign  is  included  since  this 
work  is  done  on  the  elemental  system.  The  force  on  the  element  can  be 
written  as: 


F 

e 


V  P, 


A 

e 


(48) 


where  y  is  the  coefficient  of  friction  between  the  rotor  and  stator,  P 

e 

is  the  normal  pressure  exerted  on  the  element  and  A  is  the  area  of  the 

e 

element.  The  elemental  area,  A  ,  around  the  disc  is 

e 

A^  =  2ir  r  Ar  (49) 


where  r  is  the  radius  associated  with  the  node  in  question  and  Ar  is 
the  incremental  radius.  The  pressure  on  the  element,  P  ,  can  be 
obtained  by  using  the  condition  of  constant  brake  wear  with  radius. 
Laboratory  testing  has  shown  this  to  be  a  good  assumption.  Using  this 
assumption  the  pressure  on  the  element,  P^,  is: 


P 

e 


(50) 


where  Cj  is  a  constant. 

Combining  (48),  (49),  and  (50)  gives: 


The  distance,  D,  through  which  the  force,  F^,  acts  is  given  by: 


D  =  r  tot 


where  to  is  the  angular  velocity  of  the  rotor  with  respect  to  the 
stator.  This  angular  velocity  at  any  time,  t,  can  be  expressed  by: 


.  (V)  (i _  «.) 

\R  /  \  t  / 
\  r/  \  s  / 


where  Vq  is  the  velocity  of  the  test  wheel  at  its  point  of  contact  with 
the  fly  wheel,  is  the  tire  rolling  radius  (distance  from  axle  center- 
line  to  tire/flywheel  contact  spot)  and  tg  is  the  time  required  during 
a  particular  brake  test  for  the  rotors  to  come  to  a  complete  stop. 
Equations  (47),  (51),  (52),  and  (53)  can  now  be  combined  to  give: 


We  =  -  pCi  2irAr  r/^j  /l  -  f- 


Consider  the  torque  exerted  at  the  element  during  a  brake  stop. 


T  =  F  r 
e  e 


The  total  torque  can  be  found  by  integrating  over  the  entire  wearing 
surface.  Substituting  (51)  into  (55)  and  arranging  this  in  integral 


form  gives: 


/SOR 

p  2nCi  rdr 


where  Tfs  is  the  torque  generated  at  the  interface  of  one  rotor  and  one 
stator,  SOR  is  the  stator  outside  radius,  RIR  is  the  rotor  inside 
radius  and  dr  is  the  result  of  allowing  Ar  to  become  infinitely  small. 


Upon  integration  the  torque  is  given  by: 
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(57) 


Trs  =  MitC^SOR2  -  RIR2 

The  total  torque  transmitted  by  the  brake  is  related  to  the  torque  at 
each  rotor/stator  interface  by: 


Tt  =  NTrs  (58> 

where  T  is  the  total  torque  transmitted  by  the  brake  and  N  is  the 
number  of  rotor/stator  interfaces.  Solving  (57)  for  pirCi,  substituting 
into  (54)  and  incorporating  (58)  yields  the  following: 


W 

e 


N  ^SOR2  -  RIR2) 


(59) 


The  elemental  work  term  is  now  in  terms  of  experimentally  measurable 
quantities.  This  is  the  total  work  transmitted  into  the  stator  and  the 
rotor  element.  For  this  analysis  the  fraction  of  work  input  to  the 
stator  was  taken  to  be  the  same  amount  as  that  observed  in  the  one- 
dimensional  axial  analysis.  The  axial  analysis  is  discussed  in  another 
section. 

Referring  back  to  Equation  (46) ,  it  is  necessary  to  obtain  expres¬ 
sions  for  q.  and  q  _ .  The  heat  rate  input,  q.  ,  at  radius  r  is 
in  out  in 

expressed  by  Fourier's  Law: 

’in  '  (60> 


where  K  is  thermal  conductivity  and  is  the  area  through  which  the 

heat  conduction  takes  place.  The  subscript  r  refers  to  the  location  at 
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which  the  subscripted  term  must  be  evaluated.  Writing  as  a  function 
of  radius  and  disc  thickness,  a,  gives: 


q.  =  -  Ka2Trfr  — 
in  \  3r , 


■(■ «), 


Similarly  the  rate  of  heat  flowing  out  of  an  element  at  r  +  Ar  is: 

’out  *  -Ka2’(r  !?l«r  (6: 

The  elemental  volume,  V,  can  be  expressed  as: 


V  =  2-nrAr  a 


where  a  is  the  thickness  of  the  disc  being  analyzed.  The  time  deriva¬ 
tive  term  is  now: 


pV  C  =  p2nrAraC 

at  dt 


Combining  (46)  and  (59)  through  (64)  yields: 


(Kt  £)  -  2,u(-Kt  |f) 


2  T  rAr  V  f1  t  ) 
t _ o\ _ s  / 

N  R^SOR2  -  RIR2j 


=  p  2irrAraC 


Dividing  (65)  by  2itaAr  and  taking  the  limit  of  each  term  as  Ar  goes  to 
zero  leads  to: 


r  rv  (i-r) 

t  o  v _ s  ' 

N  R^a^SOR2  -  RIR2] 


_  9T 
prC  3F 


The  first  term  can  be  expanded  to  give: 


vhere  the  chain  rule  has  been  applied  in  the  case  of  the  last  term. 

With  the  aid  of  Equation  (67)  all  terms  in  Equation  (66)  can  be 
converted  into  finite  difference  approximate  expressions.  The  expres¬ 
sions  used  here  are: 


where  T(I)  refers  to  the  temperature  of  the  I1"*1  element  and  T'(I) 
refers  to  the  temperature  of  the  I1"*1  element  at  the  next  increment  in 
time.  The  heat  capacity  and  thermal  conductivity  used  were  obtained 
from  the  manufacturer  in  the  case  of  beryllium  and  were  obtained  by 
experiment  in  the  case  of  carbon.  Appendix  K  gives  the  properties  as  a 
function  of  temperature.  Incorporating  (67)  -  (70)  into  (66)  and 
solving  for  T'(I)  gives: 
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T' (I)  =  T(I)  + 


At 

rpC 


Kr 

[T (1+1)  -  2T(I)  +  T(I-l) 


I 


+  £  [T(I+D  -  T(I-l)  ]  +  r  §  J,Cm-rT(I-1)  2 


T  V  r  ( 1  ~  t~ ) 
t  o  \  s/ 

a  it  RrN^SOR2  -  RIR2j 


(71) 


This  is  the  equation  used  to  predict  the  interior  node  temperatures. 


Boundary  Condition  Equation  Development 
The  development  to  this  point  could  be  for  a  rotor  or  for  a 
stator.  Although  Figure  59  shows  typical  values  for  a  stator,  the 
interior  node  equation  development  applies  equally  well  to  a  stator  or 
a  rotor.  On  the  other  hand,  boundary  conditions  are  dependent  on  a 
particular  selection.  Figure  60  shows  a  rotor,  stator,  and  the  inter¬ 
face  between  them.  The  stator  was  selected  for  analysis  of  the  boundary 
conditions.  The  two  elements  of  interest  are  shown  in  Figure  60. 
Considering  the  lower  element  first  and  applying  the  first  law  of 
thermodynamics  as  given  in  Equation  (46)  yields: 

K  2 nra  -  h  2irra  (t  -  T  \  =  pf1  2Trra  C  (72) 

9r  \  °°  j  i  dt 

where  h  is  the  convective/radiative  heat  transfer  coefficient,  T  is  the 

temperature  of  the  element,  the  work  input  is  zero  and  T  is  the 

ambient  air  temperature.  Note  that  the  radial  thickness  of  this 

element  is  Ar/2  as  compared  with  Ar  for  the  interior  elements.  Dividing 
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ROTOR  OUTSIDE  RADIUS  =  7.2  IN. 


STATOR  OUTSIDE  RADIUS  =  6.75  IN. 


ROTOR  INSIDE  RADIUS  =  4.4  IN. 


STATOR  INSIDE  RADIUS  =  3.75  IN. 


Figure  60.  Rotor/Stator  Interface 
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(72)  by  2irra,  denoting  the  boundary  temperature  by  T(l),  denoting  the 
adjacent  element  temperature  by  T(2) ,  and  expanding  the  right  most  term 
of  the  equation  yields: 


T(2)  -  T 


^  -  h  [t(l)  -  T J  =  PC 


where  T'(l)  is  the  temperature  of  Element  1  at  the  next  increment  in 
time.  Equation  (73)  can  now  be  solved  for  T'(l)  to  give: 


r<1)  ■  I(1>  +^c  K 


T(l)  -  T 


This  is  the  equation  used  to  calculate  the  temperature  at  the  stator 
inside  radius. 

At  the  stator  outside  radius  the  boundary  condition  is  somewhat 
different  due  to  the  work  input.  Applying  the  first  law  of  thermo¬ 
dynamics  as  given  in  Equation  (46)  yields: 


-K  2irra  —  -  h2irra 


(T  - T.)  -  r  ■  *  r  2"a  c  £ 


where  is  the  expression  given  by  Equation  (59).  If  (75)  is  now 
divided  by  2irra,  the  node  at  the  stator  outside  radius  is  denoted  by  N 
and  the  finite  difference  formulation  of  each  term  is  indicated,  the 
following  equation  results: 


-  h[T(N)  -  T^- 
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where  T’ (N)  is  the  temperature  of  the  N  element  at  the  next  increment 
in  time.  Equation  (76)  can  now  be  solved  for  T'(N)  to  give: 

T'(N)  =  T(N)  +  ^777  -  [T(N)  -  T(N-l)  ] 

-  h[T(N)  -  T„]  -  (77) 

Equation  (77),  with  the  help  of  Equation  (59)  for  W^,  was  used  for  the 
calculation  of  the  temperature  at  the  stator  outside  radius.  The 
initial  condition  requirement  was  satisfied  by  initially  setting  all 
nodes  to  a  value  dictated  by  the  experiment  being  simulated. 

A  computer  program  employing  variable  heat  capacity,  variable 
thermal  conductivity,  variable  convective/radiative  heat  transfer 
coefficients  and  variable  work  input  was  written  to  calculate  the 
temperatures  from  the  equations  derived  above.  The  program  listing  is 
given  in  Appendix  L.  Although  employing  the  equations  above  for  the 
work  input  portion  of  the  brake  stop,  the  computer  program  employs  a 
cooling  scheme  beginning  at  1.25  times  the  stop  time  that  lumps  the 
entire  rotor  and  stator  into  one  mass  at  constant  temperature  and 
allows  that  mass  to  cool.  The  radiative/convective  heat  transfer 
coefficients  used  are  given  in  Appendix  J. 

One-Dimensional  Axial 

Interior  Node  Equation  Development 
The  model  selected  for  this  axial  analysis  is  shown  in  Figure  61. 
The  nodes  and  typical  elements  of  mass  are  shown  as  well  as  applicable 


ROTOR  STATOR 

~  HALF  - mJ  -  half  - 

thickness  thickness 
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Figure  61.  One-Dimensional  Axial  Model 


dimensions.  The  rotor/stator  interface  is  shown  separated  and  iden¬ 
tified  with  two  nodes  because  the  contact  resistance  of  the  interface 


is  considered  in  the  boundary  condition  analysis.  The  starting  point 


If  the  limit  of  this  equation  is  taken  as  Ax  approaches  zero,  the 
following  equation  results: 


(83) 


Expanding  Equation  (83)  and  applying  the  chain  rule  yields: 


(84) 


Equation  (84)  can  be  placed  in  finite  difference  form  by  employing 
Equations  (68)  through  (70)  and  simply  replacing  r  with  x.  This 
yields: 


R  T(I+1)  -  2  T(I)  +  T(I-l)  +  3K  P 
(Ax)2  3T  L 


T(I+1)  -  T(I-l)  2 
2  Ax 


=  PC 


T' (I)  -  T(I) 
At 


(85) 


Solving  Equation  (85)  for  T' (I)  gives: 


T'(I) 


T^>  + 


K  T(I+1)  -  2  T(I)  +  T(I-l) 
(Ax)2 


3K  P 
3T  [_ 


T(I+1)  -  T(I-l) 
2Ax 


(86) 


\ 


This  is  the  equation  used  to  predict  the  interior  node  temperatures. 
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Boundary  Condition  Equation  Development 


The  boundary  condition  at  the  centerlines  of  the  rotor  and  stator 
were  taken  to  be  adiabatic.  This  was  handled  within  the  finite  dif¬ 
ference  technique  by  creating  an  imaginary  node  just  adjacent  to  the 
centerline  node.  In  Figure  61  the  stator  imaginary  node  was  placed  to 
the  right  of  the  stator  centerline  while  the  rotor  imaginary  node  was 
placed  to  the  left  of  the  rotor  centerline.  At  each  time  step  the 
imaginary  node  of  each  centerline  was  set  equal  in  temperature  to  the 
temperature  of  the  node  just  on  the  opposite  side  of  the  applicable 
centerline.  As  can  be  seen  from  Equation  (68),  this  procedure  insures 
that  the  temperature  gradient  at  each  centerline  is  zero. 

In  order  to  consider  the  interface  boundary  condition,  it  is 
helpful  to  expand  the  rotor/stator  interface  area  of  Figure  61  as  is 
shown  in  Figure  62.  Similar  to  the  mathematical  development  by 
Schaaf,38  the  contact  resistance  at  the  interface  is  broken  into  two 
parts,  that  associated  with  the  stator  and  that  associated  with  the 
rotor.  The  actual  contact  resistance  values  and  the  calculation  pro¬ 
cedure  are  shown  in  Appendix  M.  The  heat  resulting  from  the  work  input 
is  required  to  pass  through  the  resistances  prior  to  flowing  into  the 
rotor  or  stator.  The  actual  contact  point  of  the  rotor  and  stator  is 
modeled  as  some  temperature,  T^,  that  is  necessarily  higher  than  either 
of  the  two  adjacent  nodal  temperatures.  The  resistances  and  Tffl  are 
shown  in  Figure  62.  Although  the  figure  shows  a  separation  between  the 
two  interface  nodes,  this  is  only  for  discussion  purposes.  The  nodes 
are  allowed  to  have  different  temperatures  but  are  required  to  have  no 
distance  between  them.  This  is  an  approximation  to  the  fact  that  the 

13 
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points  would  indeed  be  separated  by  a  small  amount  but  this  amount  is 
much  smaller  than  the  Ax  values  between  the  other  nodes.  To  assist  in 
the  equation  development,  numbers  were  assigned  to  the  nodes  as  shown 
in  Figure  62.  The  first  law  of  thermodynamics  as  shown  in  Equation  (46) 
is  repeated  below  for  reference: 


q.  -q  -  M  =  pV  C  J 
in  ^out  dt 


(86) 


Equation  (46)  as  applied  to  the  crosshatched  part  of  Figure  62  between 
Nodes  11  and  12  gives: 


-  K 


\ 


9T 

8x 


'^lef  t 


\  /Tn^  T\  .  DXR  3T(12)  (87) 

—  T(!2)  ~  \  —  C  — 


where  the  first  term  is  evaluated  at  the  left  face  of  the  element,  REf 
is  the  contact  resistance  associated  with  the  rotor  and  DXR  is  the 
incremental  distance  in  the  x  direction  associated  with  the  rotor. 

This  distance,  DXR,  is  allowed  to  vary  from  the  incremental  x  direction 
distance  within  the  stator,  DXS.  Dividing  by  A^,  expanding  the  right 
hand  side  and  putting  Equation  (87)  into  finite  difference  form  yields: 


-  (I(12)  - 1(11)1  -  h;  [i<12)  -  T»] 

-p2fclz  [*'  <12)a;  I(12>]  (88) 

where  the  12  used  as  a  subscript  indicates  that  the  subscripted  quan¬ 
tity  should  be  evaluated  at  T(12).  Similarly,  applying  Equation  (46) 
to  the  crosshatched  part  of  Figure  62  between  Nodes  13  and  14  gives: 
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.  DXS  3T(13) 

pAk  T  °13  "It 


where  the  first  term  is  evaluated  at  the  right  face  of  the  element  and 
REg  is  the  stator  contact  resistance.  Dividing  by  A^,  expanding  the 
right-hand  side  and  placing  (89)  into  finite  difference  form  yields: 


m  lT(14)  - 1(13)1  +i5~  [T»- 1(13)1 

s 


=  p  —•  C 


/  T' 
13  \— 


13)  -  T(13] 
At 


Since  the  material  of  the  rotor  is  the  same  as  that  of  the  stator,  the 
contact  resistance  of  each  is  also  identical,  i.e.: 


RE  =  RE  (91) 

s  r 

This  would  not  necessarily  be  true  if  the  two  materials  were  different, 
although  Equation  (91)  could  still  be  used  if  no  experimental  data  were 
available.  Subtracting  Equation  (88)  from  (90)  while  employing 
Equation  (91)  gives: 
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L 


—  [T(14)  -  T(13) ]  +  —  [T(12)  -  T(ll) ] 


+  —  [T (12)  -  T(13) ] 


DXS 

2 


P  c13 


T* (13)  -  T(13)  \ 


At 


/ 


-  p 


DXR  I  T'  (12)  -  T(12)  ^ 


At 


(.92) 


Applying  Equation  (46)  to  the  entire  element,  the  two  halves  of  which 
are  shown  crosshatched  in  Figure  62,  yields: 


-  K  A 


I?) 


k  3x 


. 


left 


3x  J. 


-  W 


right 


3T 


"  P  \  DXA  C 


(93) 


where  W  is  the  rate  of  work  input  due  to  the  friction  at  the  interface 
and  DXA  is  an  average  between  DXR  and  DXS.  Dividing  A^,  expanding  the 
right-hand  side  and  putting  in  finite  difference  form  gives: 


_  Klo  [T,(12)  T(ll)]  T£142  -  T(_13)_  _  ° 

Kl2  DXR  Kl3  DXS  W 


p  DXA  Ca  (Z^Tfet) 


(94) 


where  W"  is  the  rate  of  work  input  per  unit  area  and  the  subscript  "a" 
indicates  that  the  subscripted  term  must  be  evaluated  at  the  average 
temperature  of  Nodes  12  and  13,  i.e.,: 


T  (a)  =  [T(13)  +  T(12) ] 


(95) 
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With  the  inclusion  of  average  temperature  definition  of  (95) 
Equation  (94)  becomes: 


Equation  (98)  can  now  be  solved  for  T' (13)  to  give: 


T' (13) 


Ik  k 

I  [T(14) " T(13)]  +  ^  [t(12)  " T(11)1 


T(12)  -  T(13)  pDXR  Cl2 
RE  2At 


*13 

T(13)  +  DXS 

K 1 2 

[T(14)  -  T(13) ]  +  DXR  [T(ll) 

-  T (12) ]  -  W" 

* 

pDXA  Ca 

2At 

f 

pDXS  T(13) 

c13 

1 pDXS  C13  pDXR  C12 

(99) 

2At 

/  2At  +  2At 

Equations  (97)  and  (99)  contain  the  rate  of  work  input  term,  W".  The 
work  rate  is  the  force  pushing  the  rotor  past  the  stator  times  the 
distance  through  which  the  force  acts  divided  by  the  time  required  to 
travel  the  distance.  The  force  can  be  expressed  as  the  brake  torque, 
Tt>  divided  by  the  rolling  radius,  Rr.  The  distance  divided  by  time  is 
the  angular  velocity,  w,  given  by  Equation  (53)  multiplied  by  the 
rolling  radius.  Th js  the  work  rate  is: 


where  N  is  the  number  of  friction  pairs.  The  work  rate  per  unit  area 


w"  = 


_ 1 

ir  SOR2  -  RIR2 


(101) 


Equations  (97),  (99),  and  (101)  were  used  to  obtain  the  rotor/stator 
interface  temperatures  when  contact  resistance  was  a  finite  value.  To 
allow  comparisons  with  contact  resistance  equal  zero,  the  temperatures 
at  Nodes  12  and  13  were  set  equal  and  Equation  (97)  was  solved  to  give: 


A  computer  program  employing  variable  heat  capacity,  variable 
thermal  conductivity,  variable  work  input,  variable  x  direction  nodal 
spacing  from  rotor  to  stator,  and  contact  resistance  was  written  using 
the  equations  derived  above.  The  program  listing  is  given  in 
Appendix  N. 


Two  Dimensional  Radial/Axial 
(Rotor/Stator  Pair) 

Interior  Node  Equation  Development 
This  analysis  was  accomplished  assuming  that  each  rotor  was 
equivalent  to  every  other  in  temperature  and  that  each  stator  was 
equivalent  to  every  other  in  temperature.  Thus,  in  order  to  analyze 
the  total  heat  sink,  analysis  of  only  one  rotor  and  one  stator  was 
sufficient.  The  model  selected  for  this  analysis  is  shown  in 
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Figure  63.  The  nodes  are  shown  as  well  as  the  different  element  types 


ROTOR  STATOR 

HALF  HALF 

THICKNESS  THICKNESS 


Figure  63.  Two-dimensional  Radial/Axial  Model 
(Rotor /Stator  Pair) 


and  dimensions.  The  actual  numbers  shown  are  those  of  the  experimental 
brake  discs.  As  in  the  previous  analysis,  the  rotor/stator  interface 
is  shown  separated  and  identified  with  two  nodes  at  each  radial  posi¬ 
tion  because  the  contact  resistance  of  the  interface  is  considered  in 
the  boundary  condition  analysis.  The  starting  point  for  the  interior 
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node  equation  development  is  again  the  first  law  of  thermodynamics  as 
expressed  in  Equation  (46)  and  as  shown  below  with  the  work  term  of 
(46)  being  zero: 


ir  /I  dT 

q .  -  q  =  pV  C  t- 

m  out  dt 


(103) 


The  heat  rate  input  of  an  interior  element  is  expressed  as: 


4in  = 

'Kx  (2,rrAr)  S 

+ 

-K  (2irrAx) 

3T 

3r 

% 

x 

_ 

■J 

(104) 


where  K  and  K  are  the  thermal  conductivities  in  the  x  and  r  direc- 
x  r 


tions,  respectively.  Similarly  the  heat  rate  output  is: 


out 


■[- 


K  (2-rrrAr)  ~1 
x  3xJx+Ax 


h 


ail 


+  l-Kr  (2-nrAx)  3rJr+Ar 


(105) 


The  volume  of  the  element  is: 

V  =  2nr  Ax  Ar  (106) 

Substituting  (104)  through  (106)  into  (103)  and  rearranging  yields: 

(k  2urAr  -  / K  2irrAr  +  (k  2urAx  -jpl 

y  x  3x/x+Ax  lx  3x/x  y  r  3ryr+Ar 

-  2nrAx  =  p2nrAxAr  C  -|^  (107) 
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Dividing  both  sides  of  this  equation  by  2nrAxAr  and  taking  the  limit  as 
Ax  and  Ar  approach  zero  gives: 


(108) 


Expanding  this  equation  and  applying  the  chain  rule  leads  to: 


Kr  3T 
r  3r 


(109) 


Equation  (109)  can  be  placed  in  finite  difference  form  by  employing 
Equations  (68)  through  (70)  for  the  r-direction  terms  and  replacing  r 
with  x  for  the  x-direction  terms.  With  these  equations  and  the  Euler 
time  marching  technique,  Equation  (109)  converts  to: 
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Solving  this  equation  for  T'(I,  J)  gives: 


This  is  the  equation  used  to  predict  the  interior  node  temperature. 

Boundary  Condition  Equation  Development 
As  in  the  previous  development  the  boundary  condition  at  the 
centerlines  of  the  rotor  and  stator  was  taken  to  be  adiabatic.  This 
was  accomplished  within  the  finite  difference  technique  by  creating 
imaginary  nodes  adjacent  to  the  centerline  nodes  and  forcing  the 
temperature  of  these  nodes,  at  each  time  step,  to  be  equal  to  the 
temperature  of  the  node  on  the  opposite  side  of  the  centerline. 

In  order  to  consider  the  interface  boundary  condition,  the  rotor/ 
stator  Interface  is  expanded  as  shown  in  Figure  64.  As  in  the  previous 
development,  the  contact  resistance  at  the  interface  is  broken  up  into 
two  parts,  that  associated  with  the  rotor  and  that  associated  with  the 
stator.  Appendix  M  gives  the  contact  resistance  values  and  calculation 
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Figure  64.  Expanded  Rotor/Stator  Interface  (2D  Radial/Axial) 


I 

procedure.  The  actual  contact  point  of  the  rotor  and  stator  is  modeled 
as  some  temperature,  T^,  that  is  necessarily  higher  than  either  of  the 
two  adjacent  nodal  temperatures.  The  resistances  and  T^  are  shown  in 
Figure  64.  Adjacent  interface  nodes  are  shown  physically  separated  for 
clarity.  The  nodes  are  allowed  to  have  different  temperatures  but  are 
required  to  have  no  distance  between  them.  This  approximates  the  fact 
that  the  distance  between  the  nodes  is  actually  much  less  than  the 
Ax  value.  The  node  identification  shown  in  Figure  64  is  employed  to 
assist  in  equation  derivation.  The  first  law  of  thermodynamics  as 
shown  in  Equation  (46)  is  repeated  below  for  reference: 


4n 


-  q  -  W  =  pV  C  — ■ 
4out  ^  dt 


(112) 


Applying  Equation  (112)  to  the  crosshatched  section  of  Figure  64 
between  nodes  (3,  J)  and  (4,  J)  yields: 


T  DXR  „  StX  /  .  DXR  _  STX 

V«4J  T  2'r  -(KR4J  T  2»r  9?)r 

-  KX4J  T^'J^DialT^:i'J>  (DR2)  2"r  +  -"-R.-— ”  (0R2)  2»r 

r 

-pSf2„rC4  (113) 

where  all  terms  except  the  first  two  have  been  placed  in  finite  dif¬ 
ference  form,  C4  is  the  heat  capacity  of  the  (4,J)  element,  KX4J/KR4J 

t 

is  the  thermal  conductivity  in  the  axial/radial  direction  of  the  4J 
element  and  DR2  is  the  distance  between  nodes  in  the  radial  direction. 
DR2  is  a  function  of  the  distance  between  the  stator  outside  radius  and 
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the  rotor  inside  radius  and  of  the  number  of  nodes  in  this  area.  Thus, 
referring  to  Figure  63: 


DR2  =  .(soR-.  mi 


(114) 


It  is  possible  that  the  distance  between  nodes  in  the  radial  direction 
is  different  in  the  area  where  the  stator  and  rotor  make  contact,  as 
opposed  to  the  areas  where  the  stator  or  rotor  are  not  in  contact. 

Divide  Equation  (113)  by  x(DR2)  and  consider  the  case  when  DR2  approaches 
zero.  This  transforms  Equation  (113)  into: 


DXR  h  (KR4J  r  If)  - 

+  2  r  ^  J<4’  ~  -  P  (DXK)r  C4  jV.&.jLj:  (lls) 

The  derivative  in  the  first  terra  of  Equation  (115)  can  be  expanded  to 
give: 

3  .  3T\  32T  .  /dl\2  3  (KR4J) 

'  — ^KR4J  r  =  KR4J  r  —  +  r  [-) 


+  (KR4J)  U 


(116) 
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Substituting  Equation  (116)  into  Equation  (115)  and  dividing  the  result 


by  r  yields: 


DXR  (KR4J)  +  DXR  ^KR4J>  +  11 

ar2  \3r/  9T  r  3r 


-  2(KX4J)  - T^3>J)  +  2  -^2 - T^'J^ 

DXR  RF. 


=  pDXR  C4 


'(4,J)  -  T(4. 


(117) 


The  first  three  terms  can  now  be  placed  in  finite  difference  form. 
This  transforms  Equation  (117)  into: 


DXR  (KR4J) 


T(4,J+1)  -  2  T(4. J)  +  T(4, J-l) 


+  pyp Ft1  ( 4 > J~H )  -  T(4.J-1)]2  9 (KR4J) 

I  2  DR2  J  9T 


+  DXR  (KR4J)  |T(4. J+l)  -  T(4,J-1) 


p  (DXR)  C4 


-  2 (KX4J) 


T(4.J)  -  T(3.J) 


Similarly  to  the  above  development.  Equation  (112)  can  be  applied  to 
the  crosshatched  section  of  Figure  64  between  nodes  (5,J)  and  (6,J)  to 
give: 


DXS  (KR5J)  T15,J+1)  r  2  T(5,J)  +  T(5tJ-l) 

DR2 


It^J+I)  -  T(5,  J-l)l  2  9  (KR5J) 
s  |_  2  DR2  J  9T 


DXS  (KR5J)  |t(5,  J+l)  -  T(5,  J-l)~l  .  T(6,J)  -  T(5,J) 

r  |_  2  DR2  J  tKXi>J)  DXS 


+  2 


T  -  T(5, J) 
m 

RE 


=  pDXS  C5 


“ 


(5,J)  -  T(5, J) 
At 


(119) 


where  KR5J  is  the  thermal  conductivity  in  the  r-direction,  KX5J  is  the 
thermal  conductivity  in  the  x  direction,  and  C5  is  the  heat  capacity  of 
the  (5,J)  element.  Apply  Equation  (112)  to  the  entire  crosshatched 
area  of  Figure  64.  The  heat  transfer  rate  terms  will  be  similar  to 
those  in  Equations  (118)  and  (119)  but  contact  resistance  terms  will 
not  appear  and  a  work  rate  term  will  appear.  For  compactness,  desig¬ 
nate  all  terms  on  the  left-hand  side  of  the  equal  sign  of  Equation  (118) 
except  the  contact  resistance  term  by  S4J.  Similarly,  designate  all 
terms  on  the  left-hand  side  of  the  equal  sign  of  Equation  (119)  except 
the  contact  resistance  term  by  S5J.  Equation  (112)  now  becomes: 


S4J  +  S5J  -  -  p  (DXR)  C4  [T' T<4'J>j  +  p  (DXS)C5 

[r^,J\;T<5,J>]  (120, 

where  is  the  rate  at  which  work  crosses  the  element  boundary  and  is 
given  by  Equation  (59).  Subtracting  Equation  (118)  from  Equation  (119) 
and  solving  for  T'(5,J)  eliminates  T  and  yields: 


T'(5,J)  =  T(5, J)  + 


At 


p  (DXS)  C5 


p  DXR  C4 
At 


[ T * (4, J)  -  T(4, J) ] 


+  S5J  -  S4J  +  [T(4, J)  -  T(5, J) ] 


(121) 


where  it  has  been  noted  that  REs  =  RE^ .  Equations  (120)  and  (121)  are 
now  two  equations  with  two  unknowns,  namely  T'(4,J)  and  T'(5,J).  If 
Equation  (12)  is  solved  for  T*(4,J)  and  this  is  substituted  into 
Equation  (121),  T'(5,J)  becomes: 


T'(5,J)  =  T(5, J)  + 


At 


2p  (DXS)  C5 


W 


2  S5J  - 


irr  (DR2) 


+  —  [T(4, J)  -  T(5, J) ] 


(122) 


Either  Equation  (120)  or  (121)  can  be  used  to  obtain  the  T'(4,J)  value 
since  T' (5, J)  is  now  known. 


In  the  case  that  contact  resistance  is  zero,  the  previous  develop¬ 
ment  is  inadequate  since  the  contact  resistance  value  appears  in  the 
denominator.  For  this  case  the  rotor/stator  interface  shown  in 


' 


Figure  64  can  be  used  as  long  as  it  is  noted  that  adjacent  surface 

nodes  are  required  to  have  the  same  temperature  and  that  T  is  no 

m 

longer  required.  With  this  in  mind,  application  of  Equation  (112) 
yields: 


— K  2irr  DR2  T^’J)  +  K  2„r  DR2  ttM)  -  I&J) 

x  DXR  x  DXS 


+  (Kr2nr  ™L  +  JZ1  .  (K  2irr  P***™  |T\  _  * 

\  2  3r/r+Ar  \  r  2  9r/r  e 

.  2ut  p  SSLt-m  DR2  c  (123) 

Similarly  to  the  previous  development,  divide  Equation  (123)  by  2irrDR2 
and  pass  to  the  limit  as  DR2  approaches  zero.  Solving  the  result  for 
T’(4,J)  gives: 


2At 


T'(4,J)  -  T(4,J)  +  p  ^DXR  +  DXS)C 


-K 


T(4 , J)  -  T(3, J) 


DXR 


,  „  T(6, J)  -  T(5, J)  /DXR  +  DXS\ 

x  DXS  V  2  / 


T(4.  T+l)  -  2  T(4, J)  +  T(4, J-l) 
DR22 


9Kr  [t(4, J+l)  -  T(4,  J-l)~l2 
3T  [  2  DR2  J 


JL  T(4, J+l)  -  T(4, J-l)\  We 

r  2  DR2  /  2ur  DR2 


(124) 
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Equations  (120),  (12'’.)  and  (124)  were  used  to  find  the  interface 
boundary  temperatures. 

The  remaining  boundary  condition  equations  were  developed  similarly 
with  areas  and  masses  changed  as  applicable.  To  account  for  the  fact 
that  some  brake  rotors  and  stators  have  wear  pads  riveted  to  the  actual 
heat  sink  material,  a  set  of  wear  surfaces  boundary  conditions  was 
derived  that  included  a  mass  for  the  wear  pads  and  a  contact  resistance 
between  the  wear  pads  and  the  heat  sink  material.  These  equations  are 
shown  in  Appendix  0.  In  some  of  these  remaining  equations,  a  radiative/ 
convective  heat  transfer  coefficient  arises.  The  values  used  and  the 
calculation  procedure  are  contained  in  Appendix  J.  The  thermal  con¬ 
ductivity  values  and  heat  capacity  values  are  shown  in  Appendix  K. 

All  of  the  equations  developed  here  were  incorporated  into  a 
computer  program  employing  variable  heat  capacity,  variable  thermal 
conductivity,  variable  radiative /convective  heat  transfer  coefficient 
and  variable  work  input.  Although  employing  the  equations  derived 
above  for  the  work  input  portion  of  the  brake  stop,  the  program  employs 
a  cooling  scheme  beginning  at  1.25  times  the  stop  time  that  lumps  the 
entire  rotor  and  stator  into  one  mass  at  constant  temperature  and 
allows  that  mass  to  cool.  The  resulting  computer  program  is  given  in 
Appendix  P. 


Three  Dimensional 
(Rotor/Stator  Pair) 

Interior  Node  Equation  Development 
The  following  analysis  again  contains  the  assumption  that  each 
stator  is  the  same  temperature  as  every  other  and  that  each  rotor  Is 


the  same  temperature  as  every  other.  This  allows  analysis  of  only  one 
stator  and  rotor.  In  the  circumferential  direction  symmetry  is  also 
utilized  such  that  only  a  sector  of  a  rotor  and  stator  must  be  analyzed. 
A  typical  rotor  and  stator  are  shown  in  Figure  65.  The  area  of  each 
required  for  analysis  is  also  shown.  An  isometric  drawing  of  the 
sectors  selected  for  analysis  is  shown  in  Figure  66.  The  nodes  used 
for  analysis  are  also  shown.  The  equation  derivation  for  the  interior 
nodes  parallels  the  two-dimensional  development  given  earlier.  The 
basic  difference  is  the  inclusion  of  the  terms  in  the  circumferential 
or  theta  direction.  Considering  the  theta  direction  in  addition  to  the 
other  two: 

’in  ’  K  r  48  41  li)x  +  K  4x  r  16  f)r 

+  (-Ke  4x  4r  I  ff)a  <125> 

q  =  (-K  rAO  Ar  |^\  +  (-K  Ax  rA0  -jp) 

out  \  x  3x/x+Ax  \  r  3r/r+Ar 

+  (-K0  Ax  Ar  i  ||)e+A0  (126) 

The  volume  of  the  interior  element  is: 


V  =  rA0  Ar  Ax 


(127) 
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Substituting  Equations  (125)-(127)  into  Equation  (103),  dividing  by  the 
volume,  and  taking  the  limit  as  Ax,  Ar,  and  A0  approach  zero  yields: 


=  P  C 


3T 


3t 


(128) 


The  same  techniques  described  earlier  can  be  used  to  place  Equation  (128) 
into  finite  difference  form  and  solve  it  for  the  temperature  at  the  new 
time  increment.  This  results  in: 


At 


T'(I,J,K)  =  T(I,J,K)  + 


K  T(I+1,J,K)  -  2  T(I,J,K)  +  T(I-1,J,K) 
X  (Ax)2 


,  9Kx  [~T(I+1,J,K)  -  T(I-1,  J ,K)~| ‘ 
31  [  2  Ax  j 


+  K  T(I.J-H.K)  -  2  T(I,J,K)  +  [T(I, J-1,K)] 
r  (Ar)2 


3K 


3T 


v_  ] 


2Ar 


T(I,J+1,K)  -  T(I , J-1,K)  2 


K 


_r  T(I,J+1,K)  -  T(I,J-1,K) 
r  2Ar 


K 


+  _i  T(I,J,K+1)  -  2  T(I, J,K)  +  T(I, J,K-1) 


(A0)  ‘ 


1  8Ke  fT(I,J,K-H)  -  T(I,J,K-1)12  1 


Equation  (129)  was  used  to  predict  the  interior  node  temperatures. 

Boundary  Condition  Equation  Development 
The  boundary  condition  at  the  axial  centerline  of  the  stator  and 
rotor  was  taken  to  be  adiabatic.  This  boundary  condition  accounts  for 
32  of  the  total  of  196  boundary  points  shown  in  Figure  66.  The 
remaining  boundary  points  were  divided  into  two  basic  groups,  those 
containing  a  work  input  term  and  those  not  containing  a  work  input 
term. 

The  boundary  condition  group  without  a  work  input  term  comprises 
104  equations.  A  typical  one  of  this  group  will  be  derived  here. 
Consider  the  element  at  I  =  2,  J  =  8,  and  K  =  5.  The  heat  transfer 
rate  in  is: 


‘in 


_  /  Ax  Ar  1  3T  \  (  rA0_  Ax  3T\ 

\  0  2  2  r  30/0  \  r  2  2  3r/r 


(130) 


The  heat  transfer  rate  out  is: 


out 


-  K  r  tJr  S)x +  h  t- l1'2-8-5’  -  TJ  <m) 


The  volume  of  the  element  is: 


Ax  Ar  rAQ 
2  2  2 


(132) 


Substituting  Equations  (130)  -  (132)  into  Equation  (103)  and  dividing 
by  the  volume  yields: 
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(133) 


Expressing  Equation  (133)  in  finite  difference  form  gives: 


ov 

9  T(2,8,5)  -  1(2,8, 4)  _  _r  T(2,8,5)  -  T(2,7,5) 

r2A0  A6  r  Ar 


2K 

x  T(3,8, 5)  -  T(2, 8, 5)  _  2h 
Ax  Ax  Ar 


T(2,8,5) 


-  „  r  T' (2,8,5)  -  T(2, 8, 5) 
"  P  C  At 


(134) 


Solving  for  T' (2,8,5)  results  in: 


T 1  (2,8,5) 


T(2,8,5)  +  ~ 
P  C 


K 


6 


r2A0 


T(2,8,5)  -  T(2,8,4) 
A0 


^  T(2,8,5)  -  T(2, 7,5)  +  \  T(3,8,5)  -  1(2, 8, 5) 
Ar  Ar  Ax  Ax 


h_ 

Ar 


(135) 


Equation  (135)  was  used  to  predict  the  temperature  at  I  =  2,  J  =  8,  and 
K  =  5. 


The  boundary  condition  group  containing  a  work  input  term  com¬ 


prises  30  nodes  on  each  disc.  An  additional  complication  arises  in  the 
development  of  the  applicable  equations  due  to  the  rotation  of  the 


■r - r 


rotor  relative  to  the  stator.  Considering  one  node  on  the  stator,  at 
each  increment  of  time,  a  different  node  from  the  rotor  will  be  adja¬ 
cent  to  this  stator  node.  It  is  necessary  to  select  the  time  increment 
such  that  at  each  point  in  time  at  which  the  equations  are  used,  the 
nodes  on  the  rotor  exactly  line  up  with  the  nodes  on  the  stator.  It  is 
not  desirable  to  apply  a  boundary  condition  equation  when  the  rotor 
nodes  are  physically  between  the  stator  nodes.  This  lining-up  problem 
can  be  solved  by  using  a  variable  time  step  such  that  the  rotation 
increment  is  a  constant.  The  variable  time  step  is  necessitated  by  the 
fact  that  the  rotor  rotation  speed  is  variable  during  a  brake  stop. 

The  following  equation  shows  the  expression  employed  to  calculate  the 
variable  time  step  value: 

R  A0 

**■  -p,  (136) 

7 

Since  the  node  numbers  designated  in  Figure  66  are  fixed  to  each  disc, 
even  during  rotation,  it  is  necessary  to  develop  a  technique  for 
describing  which  rotor  wear-surface  node  lines  up  with  each  stator 
wear-surface  node.  The  circumferential  direction  is  the  only  direction 
requiring  this  analysis  since  lines  upon  which  K  is  constant  are  radial 
lines  and  since  relative  motion  is  only  in  the  circumferential  direc¬ 
tion.  A  pictorial  view  of  the  rotor  sector  rotating  past  a  stator 
sector  at  each  time  step  is  shown  in  Figure  67.  At  each  time  that  the 
set  of  rotor  nodes  is  directly  adjacent  to  the  stator  nodes,  an  equa¬ 
tion  is  desired  that  will  predict  the  temperature  of  the  nodes.  In 
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Figure  67.  Rotor  Sector  Positions  vs  Time  Increment 


order  to  avoid  writing  equations  for  eight  rotor  positions  as  well  as 
for  all  wear-surface  nodes,  an  equation  can  be  written  for  each  wear- 
surface  node  of  the  stator  and  a  separate  equation  can  be  written  for 
each  wear-surface  node  of  the  rotor.  Then,  when  a  wear-surface  node  of 
the  rotor  lines  up  with  a  wear-surface  node  of  the  stator,  the  equa¬ 
tions  can  be  added  and  then  solved  for  the  new  temperature.  Note  that 
the  new  temperature  of  the  rotor  wear-surface  node  will  be  identical  to 
the  new  temperature  of  the  stator  wear-surface  node  since  contact 
resistance  is  not  used  in  this  model.  A  typical  equation  for  a  rotor 
wear-surface  node  will  now  be  derived.  Consider  the  element  at  I  =  4, 

J  =  5,  K  =  2.  (See  Figure  66.)  The  heat  transfer  rate  into  this 
element  is: 


=  -  K  rA0Ar 
x 


T(4,5, 2)  -  T(3,5,2) 
Ax 


Ax 
r  2 


The  heat  transfer  rate  out  of  this  element  is: 


(137) 


1  aTN 

r  96  je 


(138) 


Applying  the  first  law  of  thermodynamics,  as  given  in  Equation  (46),  to 
the  element  being  considered  leads  to: 


+  (h  f  4r  7  Ifje+Ao  -  »»  -  p  ¥  "9  “  c  ||  (139) 

where  W4  is  the  rate  at  which  work  is  input  to  this  element.  This  work 
rate  term  is  unknown  but  it  will  shortly  be  shown  that  it  is  not  needed 
in  this  form. 

Ax 

If  Equation  (139)  is  now  divided  by  y-  rA0Ar  and  the  limit  of  each 
term  is  taken  as  A0  and  Ar  approach  zero,  the  following  results: 


ill 1  T(4,5,2)  -  T(3, 5 , 2)  13/  9T\  1  3  /  3T\ 

Ax  Ax  r  9r  y  r  9r/  ^0  \  0  90  / 


T’ (4,5,2)  -  T(4, 5, 2) 


(140) 


where  W't,  is  still  an  unknown  and  the  right-hand  side  of  the  equation 
has  been  expanded  as  described  previously. 

Consider  the  element  at  I  =  5,  J  =  5,  and  K  =  4.  (See  Figure  66.) 
The  node  of  this  element  is  a  stator  wear-surface  node  at  J  =  5  and 
will  at  various  times  during  the  brake  stop,  be  adjacent  to  the  rotor 
wear-surface  node  at  I  =  4,  J  =  5,  and  K  =  2  (previously  analyzed). 

For  this  stator  node  the  heat  transfer  rate  into  the  element  is: 
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The  heat  transfer  rate  out  of  this  element  is : 


q  =  -  K  rA0  Ar 
out  x 


T(6,5,4)  -  T(5,5,4) 
Ax 


-  K 


Ax 


rA0 


9T\ 

9r/r+Ar 


9t\ 

90/0+A0 


(142) 


As  before,  if  Equations  (141)  and  (142)  are  substituted  into 
Equation  (46)  the  following  results: 


K  rA0  Ar 
x 


T(6, 5,4)  -  T(5, 5, 4)  { 
Ax 


r  i) 


r+Ar 


L  Ax  4  1  dl\  L  Ax  AQ 

\K0  2  Ar  r  90/ 0+A0  \  r  2  r  9  9r/i 
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(143) 


where  W5  is  the  rate  at  which  work  is  input  into  this  element.  At  this 
point  W5  is  unknown  and  will  be  discussed  later.  Dividing  Equation  (143) 
by  rA0  Ar  and  taking  the  limit  of  each  term  as  A0  and  Ar  approach 
zero  yields: 
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2 _ x  T(6,5,4)  -  T(5. 5.4) 


I  JL  L  r  11\ 

r  3r  \  r  3r J 


13  /„  3T\  _  T’ 

r2  36  \K0  30/  W  5  pC 


(5,5,4)  -  1 
At 


(144) 


where  W' 5  is  an  unknown  and  the  right  hand  side  of  the  equation  has 
been  expanded. 

Considering  Figure  67,  it  can  be  seen  that  at  time  increment 
number  2  the  rotor  K  =  2  nodes  are  adjacent  to  the  stator  K  =  4  nodes. 
Thus  the  rotor  wear-surface  node  at  I  =  4,  J  =  5,  K  =  2  is  adjacent  to 
the  stator  wear-surface  node  at  I  =  5,  J  =  5,  K  =  4.  In  order  to 
determine  the  temperature  of  these  two  nodes  at  this  time  increment, 
Equations  (140)  and  (144)  are  added.  As  a  shorthand  notation  all  terms 
on  the  left-hand  side  of  Equation  (140)  except  the  work  rate  term  will 
be  designated  LHSR  while  all  terms  on  the  left-hand  side  of 
Equation  (144)  except  the  work  rate  term  will  be  designated  LHSS.  In 
addition,  the  following  definitions  are  used: 


RHSR  =  Cu  5  2 
At  z 

RHSS  =  2—  C5  5  4 

At 


(145) 


(146) 


where  the  heat  capacity  is  to  be  evaluated  at  the  temperature  at  the 
node  indicated  by  the  subscripts.  Using  the  above  shorthand  notations 
the  sum  of  Equations  (140)  and  (144)  becomes: 


LHSR  + 


LHSS  -  (w'4  +  W'5)  = 


RHSR  [T* (4,5,2)  -  T(4,5,2)] 


+  RHSS  [T T (5,5,4)  -  T(5,5,4)] 


(147) 


Note  that  the  sum  of  the  two  fractional  work  rate  terms  is  the  total 
work  rate  input  at  the  rubbing  interface  and  has  been  previously 
derived  as  Equation  (59).  It  is  not  necessary  to  know  which  portion  of 
this  term  is  associated  with  the  rotor  or  stator.  Note  also  that  the 
temperature  of  the  rotor  and  stator  must  be  identical  at  each  contact 
point,  or: 


T'(4,5,2)  =  T'(5,5,4)  (148) 

Equation  (147)  can  now  be  solved  for  T' (4,5,2)  [or  T1 (5,5,4)]  to  give: 

T' (4,5,2)  = 

[rHSR  T(4, 5, 2)  +  mss  T(5, 5,4)  +  LHSR  +  LHSS  -  (w\  +  W'5)] 

(RHSR  +  RHSS) 

(149) 

Equation  (149)  was  used  to  predict  the  temperature  at  time  step 
number  2  for  the  stator  node  at  1=5,  J=5,K=4,  and  for  the  rotor 
node  at  I  =  4,  J  =  5,  K  =  2. 

All  remaining  boundary  condition  equations  were  developed  simi¬ 
larly  with  areas  and  masses  changed  as  applicable.  These  equations  are 
given  in  Appendix  Q.  As  was  the  case  previously,  radiative/convective 
heat  transfer  coefficients  used  in  this  development  are  given  in 
Appendix  J  while  the  thermal  conductivity  and  heat  capacity  values  are 
given  in  Appendix  K. 

All  equations  developed  for  this  three-dimensional  model  were 
incorporated  into  a  computer  program  employing  variable  heat  capacity. 
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variable  thermal  conductivity,  variable  radiative/convective  heat 
transfer  coefficient  and  variable  work  input.  The  program  uses  the 
equations  derived  above  for  the  work  input  portion  of  the  brake  stop 
but  uses  a  cooling  scheme  beginning  at  1.25  times  the  stop  time  that 
lumps  the  entire  mass  of  the  rotor  and  stator  into  one  mass  at  constant 
temperature  and  allows  that  mass  to  cool.  The  resulting  computer 
program  is  given  in  Appendix  R. 

Three  Dimensional  Caliper 

Interior  Node  Equation  Development 
The  interior  equation  for  the  caliper  brake  is  valid  in  the  rotor, 
lining  and  housing.  The  model  selected  for  the  caliper  brake  is  shown 
in  Figure  68.  The  interior  equation  development  is  identical  to  that 
of  the  previous  section.  Thus  Equation  (129)  is  the  applicable  equa¬ 
tion  for  predicting  the  interior  node  temperatures  for  the  caliper 
brake. 


Boundary  Condition  Equation  Development 
The  boundary  conditions  are  split  into  three  groups;  those  at  the 
lining  and  housing  interface,  those  at  the  lining  and  rotor  interface, 
and  those  open  to  the  ambient  air.  Consider  a  node  at  the  left  lining 
and  housing  interface.  If  the  node  arrangement  is  selected  as  in 
Figure  69,  it  is  noted  that  some  nodes  lie  directly  on  the  lining  and 
housing  interface.  Considering  the  housing  as  Material  1  and  the 
lining  as  Material  2,  the  heat  transfer  rate  into  an  element  at  the 
lining  and  housing  interface  (small  crosshatched  element  in  Figure  69) 
is : 
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Figure  68.  Caliper  Brake  Model 
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(150) 


where  the  subscripts  indicate  the  applicable  material.  Similarly,  the 
heat  transfer  rate  out  of  an  element  at  the  lining  and  housing  inter¬ 
face  is: 


%ut  =  ~[Ke1  2  Ar  r  'e)e+A0-(K02^-Arl||)9+A0 


Kri^irAe|l\  _  /K  _^rA0 

1  2  3r/r+Ar  \  r2  2  3r/r+Ar 


K  Ar  rA0  .  . 

x2  3x/x+Ax 


(151) 


With  and  4out  as  shown  in  Equations  (150)  and  (151)  ,  the  first  law 
of  thermodynamics  as  given  in  Equation  (103)  yields: 


q.  -  q  _  =  pi  rA0  Ar  —x~  Ci 
out  1  9  1 


1  T'  _  T 


Ax2  t,  _  T 
+  p2  rA0  Ar  —  C2  — ^ — 


(152) 
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Combining  Equations  (150)  through  (152),  dividing  by  A0Ar,  and  taking 
the  limit  of  each  term  as  A0  and  Ar  approach  zero  gives: 


Ax^  3_ 

2r  80 


Ax2  3 
2r  80 


+ 


8_ 
2  8r 


Ax i  Ax2 

pi  -y-  Ci  +  p2  —  c2 
2 


(153) 


Placing  Equation  (153)  in  finite  difference  form,  using  techniques 
described  earlier,  and  solving  it  for  the  new  temperature  results  in: 


2At 


T'  (I,  J,K)  =  T(I,  J,K)  +  — - 1 ^ — t - —~r 

r  pi  Axi  Ci  +  p2  Ax2  C2 | 


'(4*>  %  +  S2) 

2r 


T(T,J,K+1)  -  2  T(I, J,K)  +  T(I, J,K-1) 
r2(A0)2 


+ 


T(I,J,K+1)  -  T(I, J,K-1) 
2rA0 


2 


Equation  (154)  is  the  equation  for  predicting  the  temperature  of  a  node 
between  the  left  housing  and  the  lining. 

The  crosshatched  element  at  the  interface  of  the  lining  and  rotor 
(Figure  69)  is  shown  in  an  expanded  view  in  Figure  70  and  will  be  used 
in  the  derivation  of  the  equation  predicting  temperature  for  the  two 
interface  nodes.  Two  nodes  are  used  in  this  case  to  allow  a  difference 
in  temperature  of  the  two  surfaces.  A  contact  resistance  term  will 
necessarily  appear  in  the  equations.  The  expression  for  the  heat 
transfer  rate  into  the  lining  element  is: 
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Figure  70.  Expanded  Lining/Rotor  Interface 
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where  the  subscript  2  refers  to  material  2.  The  heat  transfer  rate  out 
of  the  lining  element  is: 


out 


/  Ax2  k  StA 
lK62  ~T~  Ar  r  30/e+A0 


I  Ax2  3T  2\ 

Ik  rA0  — — 1  ,  . 

\r2  2  3r  J r+Ar 


(156) 


The  first  law  of  thermodynamics  (Equation  103)  as  applied  to  the  lining 
element  is: 


Ax2  T'2  -  T2 

^in  "  %ut  =  p2  rAS  - At 


(157) 


Combining  Equations  (155)  -  (157) ,  dividing  by  A0  and  Ar,  and  taking 
the  limit  as  A0  and  Ar  approach  zero  yields: 


AX2  ^ 
2r  30 


V  3T2\  +  AX2  3  U  3T2\ 
y  02  307  +  T  37  \  r2  r  -to) 


/  3T2\  /  \  Ax2  /T’2-T2> 

-  ^x2  r  Wk  +  K  "  r  P2  —  C2 


At 


(158) 
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Similarly  the  equation  for  the  rotor  element  of  Figure  70  is: 


2r  39  \  03  39  /  +  2  3r 

+  m  (vT3)  '  r  03  ~TT  °3 


r 


T'3-T3 

At 


3T3\ 

r  3x  /x+Ax 


(159) 


For  compactness  designate  all  terms  on  the  left-hand  side  of 
Equation  (158)  except  the  contact  resistance  terra  by  SHI.  Similarly 
designate  all  terms  on  the  left-hand  side  of  Equation  (159)  except  the 
contact  resistance  term  by  SH2.  Equation  (46)  or  equivalently 
Equation  (112)  can  now  be  applied  to  the  combination  of  the  lining  and 
rotor  elements  (Figure  70)  to  give: 


W 


SHI  +  SH2 


Ar  A0 


AX2 

=  P  2  ~T~  r  c2 


3T2 

3ir  +  p3 


C3 


3T3 

alT 


(160) 


where  the  fact  that  Ar  and  A9  have  approached  zero  is  no  problem  since 
Wg  is  derived  similarly  and  is: 


W 

e 


A9  Ar  T.  r  V 
_ t  o 

N  n  (rilo2  ' 


(161) 
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Denoting  the  work  term  of  Equation  (160)  by  W'  and  expanding  the  right- 
hand  side  leads  to: 

•  AX2  T  *  2  —  T2 

SHI  +  SH2  -  W'  =  p2  — 2 —  r  C2  - ^ - 

AX3  T 1 3  —  T3 

+  P 3  ~Y~  r  c3  - ^ -  (162) 

Subtract  Equation  (159)  from  Equation  (158)  to  get: 

v  AX2  T 1 2-T2 

T3-T2)  -  P2  —  r  c2  - XT 

Ax3  T' 3-T3 

“  P 3  —  r  C3  - —  (163) 

Equations  (162)  and  (163)  are  now  two  equations  with  two  unknowns,  T'3 
and  T'2*  Substituting  Equation  (163)  into  Equation  (162)  and  solving 
for  T'2  gives: 

T'*  *  T*  +  p2  4‘r  C2  {2SH1  -  +  k  (T3-T2)}  (164> 

T'3  can  be  calculated  using  either  Equation  (162)  or  Equation  (163). 
Equations  (162)  through  (164)  are  the  equations  for  predicting  the 
temperature  at  the  interface  of  the  left  lining  and  the  rotor.  The 
equations  for  the  right  side  of  the  brake  are  similar  except  for  sub¬ 
scripts  and  the  sign  of  the  x-direction  heat  transfer  rate  term. 
Equations  (162)  through  (164)  must  be  applied  to  the  brake  with  care 
since  there  exists  relative  motion  between  the  rotor  and  lining.  As  in 
the  previous  development,  a  variable  time  step  is  desired  so  that  the 
rotation  increment  during  one  time  step  will  allow  nodes  on  the  rotor 


SHI  -  SH2  +  —■  (' 
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to  exactly  line  up  with  nodes  on  the  lining.  This  variable  time  step 
is  calculated  by: 


R  A0 


In  addition  to  the  variable  time  step.  Equations  (162)  through  (164) 
must  also  contain  a  provision  for  determining  the  correct  circumferen¬ 
tial  position  when  using  the  old  (not  primed)  temperature  values.  This 
can  be  accomplished  by  subscripting  the  old  temperature  values  such 
that  the  subscripts  correspond  to  the  angle  traversed.  The  relation 
between  the  angle  traversed  and  the  time  increment  is  Equation  (165). 

It  should  be  noted  that  the  somewhat  involved  procedure  for 
equation  addition  within  a  computer  program  as  developed  in  the  previous 
section  is  not  needed  for  this  development  since  all  elements  are 
identical  in  the  rubbing  interface  area.  All  boundaries  not  discussed 
above  were  taken  to  be  adiabatic  since  previous  results  showed  little 
cooling  during  the  work  input  portion  of  the  brake  stop.  Cooling  of 
the  brake  after  rotation  is  stopped  can  again  be  handled  assuming  a 
lumped  mass  cooling  at  constant  temperature. 

No  computer  program  was  developed  for  the  caliper  brake  for  a 
number  of  reasons.  First,  no  tests  were  conducted  within  this  effort 
to  allow  direct  comparison  of  theory  to  experiment.  Second,  the 
experimental  data  from  other  authors  is  fairly  limited.  In  addition, 
brake  geometries  and  material  properties  of  other  authors  are  almost 
nonexistent.  Last,  but  certainly  not  least,  is  the  fact  that  the  same 
techniques  involved  in  the  caliper  brake  equation  development  were  used 
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in  the  multidisc  brake  equation  development.  This  implies  that  the 
conclusions  applying  to  the  multidisc  brake  theory  also  apply  to  the 
caliper  brake  theory. 

Two  Dimensional 
(Entire  Brake/Wheel/Tire) 

This  analysis  is  such  that  an  axial  temperature  gradient  can  exist 
from  stator  to  stator  or  rotor  to  rotor.  The  model  selected  for  this 
two-dimensional  multidisc  analysis  is  shown  in  Figure  71.  The  numbered 
areas  apply  to  the  nodes  assigned  to  that  area.  The  torque  tube  is 
divided  into  six  areas  with  node  numbers  1,  9,  21,  22,  23,  and  24.  The 
rim  is  divided  into  seven  areas  with  node  numbers  3,  7,  8,  11,  12,  13, 
and  14.  The  tire  is  shown  as  two  nodes,  4  and  6.  The  remaining  two 
nodes  are  Number  5,  the  contained  air  and  Number  2,  the  housing  and 
hydraulic  fluid  node.  The  heat  sink,  all  rotors  and  stators,  is 
divided  into  eight  areas  bounded  by  the  radial  center  lines  of  an 
adjacent  rotor  and  stator.  One  of  these  eight  areas  is  shown  cross- 
hatched  in  Figure  71.  Each  of  the  eight  heat  sink  areas  is  analyzed 
independently  using  primarily  the  equations  derived  in  the  two- 
dimensional  radial/axial  section  of  this  report.  The  nodes  utilized 
for  each  of  the  eight  areas  are  necessarily  the  same  as  those  used  in 
the  two-dimensional  radial/axial  section.  The  nodal  pattern  employed 
is  shown  in  Figure  63.  For  brakes  containing  wear  pads  at  the  rotor/ 
stator  interface,  the  nodal  pattern  is  shown  in  Appendix  0  as 
Figure  0-2.  Since  each  rotor  and  stator  is  in  contact  with  either  the 
rim  or  the  torque  tube,  a  heat  transfer  term  was  added  to  allow  heat 
flow  from  the  heat  sink  to  the  rim  and  torque  tube. 
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Figure  71.  Multi-Disc 


Each  numbered  node  shown  in  Figure  71  was  treated  as  a  constant 
temperature  lump  of  mass  with  thermal  properties  based  on  the  material 
being  depicted.  As  an  example  of  the  method  of  equation  development 
for  each  numbered  node,  consider  Node  12.  This  node  has  been  expanded 
and  the  heat  transfer  rates  for  the  node  have  been  identified  in 
Figure  72.  The  convection- radiation  heat  transfer  rate  from  Node  12 
was  obtained  by: 

QHR  =  h  x  A12  ^T12-T^  (166) 

where  h  is  the  convective- radiative  heat  transfer  coefficient  obtained 
as  shown  in  Appendix  J,  A12  is  the  area  of  Node  12  open  to  the  ambient 
air,  T12  is  the  temperature  of  Node  12  and  T^  is  the  ambient  air 
temperature.  The  conduction  heat  transfer  rate  from  Node  12  to  Node  4 
is  found  by  considering  the  internal  resistances  of  the  blocks  of  mass 
associated  with  Nodes  12  and  4.  The  resistance  associated  with  Node  12 
for  this  application  is: 


RE12  = 


LRIO 


VKR(ARIO) 


(167) 


where  RE12  is  the  resistance,  LRIO  is  the  length  from  the  center  of 
Node  12  to  the  contact  point  between  the  two  masses,  VKR  is  the  thermal 
conductivity  of  the  Node  12  (rim)  material,  and  ARIO  is  the  area  of 
contact  between  the  two  masses.  The  resistance  associated  with  Node  4 
for  this  application  is: 


RE4 


LTT4 

VKT(ARIO) 


(168) 
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Figure  72.  Node  12  Heat  Transfer  Rates 
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where  RE4  is  the  resistance,  LTT4  is  the  length  from  the  center  of 
Node  4  to  the  contact  point  between  the  two  masses,  VKT  is  the  thermal 
conductivity  of  the  Node  4  (tire)  material,  end  ARIO  is  again  the  area 
of  contact  between  the  two  masses.  The  heat  transfer  rate  between 
Nodes  12  and  4  can  now  be  written  as: 


Q124  =  (T12-T4) 


(169) 


Substituting  Equations  (167)  and  (168)  into  Equation  (169)  gives: 


Q124  =  (T12-T4) 


ARIO 

LRIO  LTT4 
VKR  VKT 


(170) 


The  other  heat  transfer  rates  shown  on  Figure  72  are  developed  simi¬ 
larly.  For  Node  12,  the  first  law  of  thermodynamics  can  be  applied  to 
give: 


Solving  for  the  temperature  at  the  time  t  +  At,  T12',  in  terms  of  the 
temperature  at  time  t,  T12,  yields: 

T12 '  =  T12  +  [QR(3)  +  QR(4)  +  Q1312  -  Q1211 

-  Q124  -  QHR]  (172) 

This  equation,  along  with  the  similar  equations  for  the  other  numbered 
nodes  and  the  heat  sink,  was  used  to  predict  the  temperatures  within 


the  entire  brake,  wheel,  and  tire. 


All  remaining  equations  for  the  numbered  nodes  are  shown  in 
Appendix  S.  The  radiative/convective  heat  transfer  coefficients  used 
are  given  in  Appendix  J  while  the  thermal  conductivity  and  the  heat 
capacity  values  are  given  in  Appendix  K. 

All  equations  developed  for  this  two-dimensional  model  of  the  heat 
sink,  tire,  rim,  torque  tube,  and  brake  housing  were  incorporated  into 
a  computer  program  employing  variable  heat  capacity,  variable  thermal 
conductivity,  variable  work  input  and  variable  geometry.  The  program 
uses  the  equations  derived  above  during  the  work  input  portion  of  the 
brake  stop  but  uses  a  cooling  scheme  beginning  at  1.25  times  the  stop 
time  that  lumps  the  entire  heat  sink  mass  into  one  mass  at  constant 
temperature  and  allows  that  mass  to  cool.  During  the  heat  sink  cooling 
period,  the  masses  of  the  numbered  nodes  in  Figure  71  simultaneously 
absorb  heat  from  the  heat  sink  and  lose  heat  to  the  surrounding  air. 

The  computer  program  resulting  from  the  above  equations  is  given  in 
Appendix  T. 
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CHAPTER  VI 

CORRELATION  OF  RESULTS 

Introduction 

Five  computer  programs  were  used  to  predict  the  temperature  at 
various  locations  within  the  brake  and  surrounding  hardware.  As  a 
shorthand  notation  the  programs  and  results  are  referred  to  as: 

One-Dimensional  Axial,  1-DA  or  1-D  Axial 

One-Dimensional  Radial,  1-DR  or  1-D  Radial 

Two-Dimensional  Radial/Axial  (R/S  Pair),  2-DRA  or  2-D  Radial/Axial 

Two-Dimensional  Entire  Brake,  2-D 

Three-Dimensional  (R/S  Pair) ,  3-D 

Many  of  the  required  computer  program  inputs  are  given  and  explained 
in  the  program  listings  in  the  appendices  and  in  the  form  of  graphs 
within  other  appendices.  For  reference.  Table  9  shows  some  selected 
program  input  values.  The  radii  of  the  carbon  heat  sink  are  shown  to 
be  different  than  those  of  the  beryllium  heat  sink  even  though  they  fit 
in  the  same  brake  housing  and  rim.  The  beryllium  disc  outside  radii 
are  smaller  than  the  corresponding  carbon  disc  outside  radii  and  the 
beryllium  disc  inside  radii  are  larger  than  the  corresponding  carbon 
disc  inside  radii.  These  values  were  used  in  the  theoretical  models  to 
account  for  the  fact  that  the  beryllium  discs  had  cutouts  at  various 
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SELECTED  COMPUTER  PROGRAM  INPUT  VALUES 
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TABLE  9 — continued 


locations  around  the:  circumference  (see  Figures  21  and  22);  thus,  in 
order  to  maintain  the  proper  thickness  and  mass  for  a  beryllium  disc, 
the  radii  were  modified. 

Contact  resistance  was  included  in  the  boundary  condition  equa¬ 
tions  developed  for  predicting  the  temperature  at  the  interface  of  a 
stator  and  rotor.  This  contact  resistance  term  was  included  in  develop¬ 
ments  1-DA,  2-DRA  and  in  the  development  for  the  caliper  brake.  No 
graphs  or  curves  are  presented  showing  this  contact  resistance  effect 
because  it  was  determined  that  less  than  a  1%  difference  in  peak 
temperature  occurred  at  the  contact  resistance  values  obtained  from  the 
disc  surface  measurements.  Even  when  this  contact  resistance  was 
multiplied  by  ten,  the  peak  temperature  changed  by  only  about  2%.  Thus 
for  the  models  2-D  and  3-D  the  rubbing  surface  contact  resistance  was 
not  used.  Since  the  rubbing  surface  contact  resistance  does  not  take 
into  account  contact  resistance  between  the  wear  pads  and  the  actual 
disc  material  (see  Figures  21  and  22) ,  a  separate  development  was 
employed  to  include  this  effect.  This  effect  is  significant  and  is 
discussed  in  more  detail  in  a  subsequent  section. 

Temperature-Time  History,  Experiment,  1-DA,  1-DR,  2-DRA 

For  the  carbon  brake,  the  thermocouple  at  a  radial  position  of 
5.5  in.  and  on  the  axial  centerline  of  Stator  2  (see  T/C  No.  3  on 
Figures  30  and  32)  was  chosen  as  a  typical  heat  sink  temperature. 

Figure  73  shows  a  plot  of  temperature  versus  time  comparing  experiment 
with  1-DR,  1-DA,  and  2-DRA.  The  time  that  the  brake  stops  rotating  is 
noted.  In  1-DA  there  is  no  radial  temperature  gradient  while  in  1-DR 
there  is  no  axial  temperature  gradient.  The  fact  that  a  particular 


Figure  73.  Temperature  vs  Time  at  R  =  5.5.  (Test  1.1-3,  ID  vs  2  DRA) 


radius  and  axial  position  are  specified  by  selecting  the  experimental 
location  is  neglected  in  the  above  two  cases.  The  exact  radial  and 
axial  locations  can  be  determined  for  the  2-DRA  model.  As  noted,  the 
axial  model  overpredicts  the  temperature.  This  is  the  case  since  no 
provision  for  cooling  is  incorporated  in  the  model.  Since  no  cooling 
takes  place.  Figure  73  contains  no  cooling  curve  for  the  1-DA  model. 

The  1-DR  model  also  overpredicts  the  temperature  but  is  closer  to 
experiment  than  the  1-DA  model.  For  the  1-DR  model  cooling  is  provided 
at  the  outside  and  inside  circumference  of  the  disc.  As  discussed  in 
the  theoretical  development  section  the  1-DR  model  predicts  tempera¬ 
tures  for  the  stator  only.  In  addition  to  the  circumferential  cooling, 
cooling  occurs  at  the  face  of  the  stator  toward  the  inside  radius  at 
the  area  that  is  not  in  contact  with  the  rotor.  The  reason  for  the 
temperature  over  prediction  in  the  case  of  the  1-DR  model  is  attributed 
to  two  things.  First,  the  radial  model  work  input  term  was  treated  as 
an  internal  heat  generation  since  no  boundary  was  present  in  the  axial 
direction  for  work  input.  This  insured  no  axial  temperature  gradient. 
On  the  other  hand,  the  thermocouple  was  placed  at  the  axial  center  of 
the  disc  which  is  the  coolest  axial  position.  Thus,  a  model  without 
the  axial  gradient  would  predict  a  higher  than  measured  temperature. 

The  second  reason  that  the  1-DR  model  could  predict  higher  temperatures 
is  because  the  1-DR  model  is  linked  to  the  1-DA  model.  This  is  true 
since  the  split  of  energy  into  the  rotor/stator  interface  for  the  1-DR 
model  was  taken  to  be  that  of  the  1-DA  model.  Also  shown  on  Figure  73 
is  the  heating  and  cooling  as  predicted  by  the  2-DRA  model.  This  model 
is  much  closer  to  experiment  as  it  takes  into  account  all  required 


convective/radiative  cooling  and  accounts  for  temperature  gradients  in 
the  radial  and  axial  directions. 

Figure  74  shows  the  predictions  from  the  three  models,  1-DA,  1-DR 
and  2-DRA  for  a  brake  stop  that  occurs  in  41  seconds  as  compared  to 
58  seconds  in  Figure  73.  Also  the  energy  input  into  the  brake  on 
Figure  74  is  100  percent  of  normal  kinetic  energy  (11.23  x  106  ft-lb) 
while  that  for  Figure  73  was  60  percent.  Again  the  1-DA  model  gives 
the  highest  temperatures  of  the  three  models,  followed  by  1-DR  and 
2-DRA.  The  2-DRA  model  again  gives  best  results  and  shows  a  peak 
temperature  difference  from  experiment  of  about  4%  as  compared  with 
about  4%  for  Figure  73.  The  2-DRA  model  does  predict  higher  tempera¬ 
tures  than  experiment  during  the  first  25  seconds  of  the  stop.  This 
lack  of  correlation  is  attributed  to  thermocouple  lag  time.  The 
temperature  gradient  during  the  heat  input  for  Figure  74  is  about  32°F 
per  second  while  that  for  Figure  73  is  about  100°F  per  second. 

Considering  the  beryllium  heat  sink.  Figure  75  gives  a  36-second 
stop  and  compares  1-DA,  1-DR,  and  2-DRA  to  experiment.  The  energy 
input  for  this  case  is  60  percent  of  the  normal  kinetic  energy.  For 
the  beryllium  heat  sink,  the  thermocouple  at  a  radial  position  of 
4.4  in.  and  on  the  axial  centerline  of  Stator  2  (see  T/C  No.  3  on 
Figure  29)  was  chosen  to  compare  to  theory.  Figure  76  shows  a  second 
stop  with  the  beryllium  heat  sink.  The  stop  time  in  this  case  is 
124  seconds  while  the  energy  input  is  100%  of  normal  kinetic  energy. 
Although  the  energy  input  levels  correspond  to  those  in  Figures  73  and 
74,  the  beryllium  heat  sink  remains  at  a  cooler  temperature  than  the 
carbon  due  to  the  difference  in  thermal  properties  (see  Appendix  K) 
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between  the  two.  Again  in  both  figures  the  1-DA  model  predicts  the 
highest  temperatures,  the  1-DR  model  predicts  the  next  highest  tempera¬ 
tures,  and  the  2-DRA  model  predicts  temperatures  closest  to  experiment. 
The  peak  temperature  predicted  by  the  2-DRA  model  for  Figure  75  is 
within  2%  of  experiment  while  for  Figure  76  the  peak  temperature 
predicted  is  within  8%  of  experiment.  The  temperature  gradient  during 
the  heat  input  is  about  15°F  per  second  in  Figure  75  and  about  8°F  per 
second  in  Figure  76.  The  test  with  the  highest  temperature  gradient 
during  heat  input,  Figure  75,  exhibits  the  same  trend  of  the  experi¬ 
mental  temperatures  lagging  the  2-DRA  predictions  as  was  noted  in 
Figure  74.  Although  the  same  trend  is  noted  in  Figures  74  and  75,  the 
lag  shown  in  Figure  75  cannot  be  solely  attributed  to  thermocouple  lag 
because  of  several  differences  between  the  2-DRA  model  inputs  for  the 
two  cases.  The  beryllium  brake  disc  uses  thin  brake  pads  as  the 
wearing  surface  and  the  use  of  these  pads  creates  a  contact  resistance 
between  the  pad  and  the  beryllium  material.  This  contact  resistance 
was  assumed  to  be  constant  and  was  selected  to  be  in  the  range  normally 
found  in  this  type  application  (.005  HR-Ft2-°F/BTU) .  This  contact 
resistance  value  used  probably  varies  with  temperature  as  the  contact 
area  between  the  wear  pad  and  the  beryllium  material  varies  with 
temperature  due  to  pad  warpage.  Another  model  input  difference  in  the 
2-DRA  model  for  beryllium  versus  carbon  is  the  thermal  properties 
(Appendix  K) .  Heat  capacity  of  the  beryllium  material  is  about  one- 
half  that  of  the  carbon.  This  means  the  beryllium  is  not  able  to  store 
as  much  thermal  energy  as  the  carbon.  Also  the  thermal  conductivity  of 
the  beryllium  is  higher  than  that  of  he  carbon  so  that  the  thermal 
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resistance  is  lower  and  the  beryllium  will  therefore  lose  heat  faster 


than  the  carbon. 

The  most  important  point  to  be  made  in  this  section  is  that  for  the 
2-DRA  model,  the  maximum  temperature  was  predicted  within  4%  of  experi¬ 
ment  for  three  of  the  four  stops  studied  and  that  the  fourth  was  within 
8%.  It  is  this  maximum  temperature  that  must  be  accurate  to  be  able  to 
accurately  predict  the  temperatures  of  any  surrounding  hardware.  The 
precise  temperature  gradients  during  heating  are  not  as  important  since 
these  take  place  over  about  one  minute  compared  with  the  approximately 
ten  minutes  required  for  the  surrounding  hardware  to  reach  its  maximum 
temperature. 


Temperature-Time  History,  Constant/Variable 
Torque  and  Pads/No  Pads 

Recall  from  the  theoretical  development  chapter  that  the  work 
input  term  for  each  model  was  proportional  to  the  torque  generated  by 
the  brake  during  the  brake  stop.  The  torque  can  be  calculated  if  brake 
pressure,  deceleration,  brake  geometry,  and  coefficient  of  friction  are 
known.  This  calculation  results  in  a  constant,  average  torque  since 
brake  pressure,  deceleration  and  coefficient  of  friction  are  assumed 
constant.  Unfortunately,  brake  test  experience  shows  that  the  brake 
torque  is  not  constant  with  time.  Not  only  does  torque  vary  with  time 
but  also  it  varies  in  different  ways  for  different  brake  materials. 
Figure  77  is  a  comparison  of  average  torque  and  typical  torque  curves 
for  the  carbon  and  beryllium  brake  (for  actual  torque  curves  of  carbon 
and  beryllium,  see  Appendix  E) .  Also  in  Figure  77  with  the  same  time 
scale,  is  a  plot  of  flywheel  energy  versus  time.  It  can  be  seen  that 
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typical  beryllium  brake  torque  curve 


the  highest  energy  of  the  flywheel  and  therefore  the  highest  energy 
input  into  the  brake  occurs  during  the  first  portion  of  the  stop  and 
that  after  the  midway  point  of  the  stop  (at  the  stop  time  divided  by 
2),  the  vast  majority  of  energy  to  be  absorbed  by  the  brake  has  already 
been  absorbed.  Thus,  the  shape  of  the  torque  curve  is  most  important 
during  the  first  half  of  the  stop  and  use  of  an  average  torque  for 
input  into  a  thermal  model  would  not  produce  an  adequate  simulation. 

The  difference  in  temperatures  resulting  from  the  use  of  an  average 
versus  an  exact  torque  curve  in  the  2-DRA  model  is  depicted  in  Figure  78 
for  carbon.  The  constant  torque  maximum  temperature  is  7%  lower  than 
the  variable  torque  maximum  temperature  of  the  2-DRA  model.  Looking  at 
Figure  77,  it  would  be  expected  that  use  of  the  average  torque  curve 
for  carbon  would  result  in  a  prediction  of  lower  temperatures  but  due 
to  the  actual  torque  beginning  at  a  lower  value  than  the  average,  the 
lower  temperatures  (7%)  are  not  drastic.  A  more  drastic  difference 
should  be  noted  for  the  case  of  beryllium  since  the  actual  torque  curve 
is  always  lower  than  the  average  torque  curve  during  the  first  half  of 
the  stop  (note  Figure  77).  Figure  79  shows  a  beryllium  experimental 
curve  along  with  three  curves  from  the  2-DRA  model.  The  upper  curve  is 
the  constant  torque  curve  and  gives  a  peak  temperature  23%  higher  than 
the  actual  torque  curve,  the  curve  closest  to  the  experimental  curve. 

Also  shown  in  Figure  79  is  a  curve  generated  from  the  2-DRA  model 
assuming  no  wear  pads  on  the  beryllium  disc  but  using  the  same  equiva¬ 
lent  disc  thickness.  Thus,  the  total  disc  thickness  in  this  case  is 
found  by  adding  the  wear  pad  thickness  to  the  actual  beryllium  thick¬ 
ness.  The  entire  resulting  disc  is  assumed  to  be  beryllium.  The  input 
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Figure  78.  Temperature  vs  Time  at  R  =  5.5  in.  (Test  1.1-3  Constant  Torque  Effect) 
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torque  for  this  case  was  the  actual  torque  curve.  As  can  be  seen  in 
Figure  79,  this  type  of  2-DRA  model  is  not  very  accurate  in  that  the 
temperature  lag  evident  in  the  experimental  data  is  not  predicted  and 
in  that  the  peak  temperature  is  26°F  different  from  experiment  while 
the  2-DRA  model  including  wear  pads  was  only  10°F  different  from 
experiment . 


Temperature-Time  History,  2-DRA  and  3-D 
The  three-dimensional  temperature  prediction  model,  3-D,  of  a 
rotor  and  stator  sector  was  shown  in  Figure  66.  For  the  purposes  of 
comparing  experiment,  2-DRA  and  3-D  at  the  same  physical  location  in 
the  disc,  and  adequately  representing  the  3-D  model,  the  3-D  model  was 
modified  somewhat.  In  order  to  simulate  the  conduction  heat  transfer 
out  of  the  lug  positions  into  the  rim  and  torque  tube,  the  nodes  at 
(2,10,4),  (2,9,4),  (3,10,4),  (3,9,4),  (4,10,4),  (4,9,4),  for  the  rotor 
and  at  (5,1,4),  (5,2,4),  (6,1,4),  (6,2,4),  (7,1,4),  (7,2,4)  for  the 
stator  were  forced  at  each  time  step  to  be  a  given  percentage  lower 
than  the  actual  calculated  value  for  that  node.  Also,  during  this  time 
period  the  radiation  convection  heat  transfer  coefficient  was  taken  to 
be  zero.  This  coefficient  was  taken  as  zero  since  the  primary  heat 
loss  mechanism  for  a  disc  during  a  stop  is  conduction  through  the  lug 
areas.  The  actual  percentage  reduction  in  the  temperatures  of  the 
nodes  above  was  chosen  so  that  the  percentage  reduction  was  the  same 
for  simulations  of  the  carbon  1.1-3  test  and  the  beryllium  7.2-3  test 
and  so  that  the  resultant  peak  temperatures  of  the  simulations  corre¬ 
sponded  to  the  experimental  values  of  the  two  tests.  This  temperature 
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reduction  technique  was  used  until  1.25  times  the  stop  time  at  which 
time  the  conduction  simulation  was  terminated  and  the  disc  was  lumped 
into  one  mass  and  allowed  to  cool  by  convection/radiation.  The  result 
of  this  effort  is  shown  in  Figure  80  for  the  carbon  disc.  The  2-DRA 
and  3-D  curves  are  within  1.5  percent  of  each  other  at  the  stop  time 
and  never  are  farther  apart.  Figure  81  shows  the  results  for  a  beryl¬ 
lium  disc.  Since  the  primary  objective  of  the  3-D  model  development 
was  to  look  at  circumferential  temperature  gradients,  wear  pads  were 
not  simulated.  Thus,  for  the  comparison  of  2-DRA  and  3-D  in  Figure  81, 
neither  of  the  theoretical  models  used  wear  pads.  The  two  theoretical 
curves  are  within  5%  of  each  other  at  about  40  seconds  and  this  is  the 
largest  deviation  noted.  This  deviation  is  larger  than  that  of  the 
carbon  and  is  attributed  to  the  fact  that  the  beryllium  disc  is  not 
simulated  as  well  as  the  carbon  in  the  area  of  the  wear  pads  or  in  the 
area  of  the  radial  cutouts  (Figures  21  and  22). 

Circumferential  Direction  Temperature  Gradient,  3-D 
To  determine  the  relative  importance  of  the  circumferential 
temperature  gradient  in  the  brake  application,  the  3-D  Model  was  used 
to  generate  temperatures  in  the  circumferential  direction.  The  con¬ 
duction  heat  transfer  out  of  the  lugs  was  treated  as  previously  dis¬ 
cussed.  The  resulting  temperatures  were  plotted  for  a  carbon  brake 
stop  in  Figure  82.  The  temperatures  plotted  occurred  on  the  axial 
centerline  of  the  stator  and  at  the  radial  and  circumferential  posi¬ 
tions  shown  in  Figure  82.  Three  times  during  the  brake  stop  were 
selected  for  plotting,  namely  .25,  .5,  and  .75  times  the  stop  time.  At 
each  time,  a  curve  is  shown  of  temperature  versus  circumferential 
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Figure  81.  Temperature  vs  Time  at  4.4  in.  (Test  7.2-3,  2  DRA  vs  3D) 
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Figure  82.  Temperature  vs  Circumferential  Direction  (Test  1.1-3) 
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direction  at  the  inside  and  outside  radii  of  the  stator  The  curves  at 
the  outside  radius  show  no  temperature  variation  in  the  circumferential 
direction  but  the  curves  at  the  inside  radius  do  show  a  variation  since 
the  conduction  occurs  at  the  inside  radius.  Although  the  variation  in 
temperature  increases  with  time,  the  variation  is  at  a  constant  7%. 

The  same  type  of  analysis  except  for  the  beryllium  brake,  results  in 
curves  as  shown  in  Figure  83.  The  variation  at  the  inside  radius  is 
again  a  constant  percentage  but  the  percentage  has  increased  to  14%. 
This  occurs  since  the  heat  capacity  of  beryllium  is  about  twice  that  of 
carbon.  Although  a  temperature  gradient  does  exist  in  the  circum¬ 
ferential  direction,  it  is  not  recommended  that  each  disc  of  the  brake 
be  modeled  in  three  directions  since  two-dimensional  disc  analysis  is 
close  to  experiment  and  since  computer  time  and  storage  dramatically 
increase  when  the  third-dimension  is  added.  If,  on  the  other  hand, 
information  is  desired  in  the  lug  area,  then  the  three-dimensional 
model  can  be  used  for  one  disc  pair  only  and  the  required  information 
can  be  found. 

Disc,  Torque  Tube,  Rim  and  Contained 
Air  Temperatures,  2-D 

For  disc  temperature  comparisons  in  the  carbon,  the  thermocouple 
on  the  axial  centerline  of  the  stator  at  a  radius  of  5.5  inches  was 
used,  as  before.  The  thermocouple  in  the  beryllium  stator  was  also 
selected  as  before,  with  a  radius  of  4.4  inches  and  on  the  axial 
centerline.  Figures  84-87  show  temperature  versus  time  plots  for  the 
stator  and  torque  tube  of  two  tests  using  carbon  and  two  tests  using 
beryllium  heat  sink  materials.  The  stator  heating  curves  are  similar 
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Figure  83.  Temperature  vs  Circumferential  Direction  (Test  7.2-3) 


BERYLLIUM  8.1- 


Figure  87.  Temperature  vs  Time,  Stator-Torque  Tube  (Test  8.1-2) 


to  the  heating  curves  of  the  2-DRA  model  shown  in  Figures  73-76.  This 
is  necessarily  so  since  the  2-DRA  model  is  essentially  a  subroutine 
within  the  2-D  model.  The  difference  between  the  two  models  is  that 
cooling  of  the  disc  for  the  2-DRA  model  is  solely  by  convection/ 
radiation  and  cooling  for  the  2-D  model  also  includes  heat  conduction. 
The  fact  that  the  heating  curves  in  all  four  cases  do  not  change  from 
model  to  model  shows  the  relative  unimportance  of  the  cooling  phe¬ 
nomenon  during  the  work  input  portion  of  the  stop.  For  the  two  carbon 
tests,  the  theory  shows  a  higher  temperature  during  cooling  than  is 
shown  by  the  experimental  curve.  However  the  theoretical  and  experi¬ 
mental  curves  converge  as  time  increases.  In  the  case  of  beryllium  the 
theoretical  cooling  curves  are  almost  coincident  with  the  experimental 
curves.  A  comparison  of  the  torque  tube  temperature  versus  time  curves 
in  Figures  84-87  shows  that  the  predicted  peak  temperatures  are  within 
about  6%  of  experimental  peak  temperatures  for  the  carbon  tests  and 
within  10%  for  the  two  beryllium  tests.  The  time  at  which  the  peak 
temperatures  occur  within  the  torque  tube  is  predicted  quite  well  for 
all  four  tests. 

The  relation  between  theory  and  experiment  for  the  fuse  plug  and 
the  drive  lug  for  the  four  cases  discussed  above  is  depicted  in 
Figures  88-91.  It  is  to  be  expected  that  the  thermal  model  will 
deviate  from  experiment  more  at  greater  distances  from  the  work  input 
than  at  the  actual  work  input  location.  This  is  true  since  any  model 
inaccuracies  will  be  additive  over  distance  and  since  contact  at 
material  interfaces  is  modeled  as  the  same  temperature  for  both  mate¬ 
rials  (no  contact  resistance).  For  the  carbon  tests.  Figures  88  and 
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Figure  91.  Temperature  vs  Time,  Fuse  Pluj 


89,  the  general  shape  of  the  predicted  curves  is  the  same  as  experi¬ 
ment.  The  maximum  peak  temperature  difference  from  theory  to  experi¬ 
ment  is  less  than  2%.  The  beryllium  predictions,  Figures  90  and  91, 
are  not  as  close  to  experiment  as  are  the  carbon.  One  contributor  to 
this  fact  is  that  the  beryllium  heat  sink  is  not  modeled,  as  exactly  as 
is  the  carbon  heat  sink.  The  actual  mass  and  volume  of  the  beryllium 
was  modeled,  but  the  holes  for  the  wear  and  attachment  and  other  holes 
to  allow  for  expansion  were  not  modeled.  Thus  the  fuse  plug  is  within 
7%  of  experiment  while  for  the  case  in  Figure  90,  the  drive  lug  is 
29%  higher  than  experiment.  In  three  of  the  four  figures  shown,  the 
drive  lug  predicted  temperature  is  within  2%  of  the  experiment. 

The  relation  between  the  2-D  predictions  for  temperature  in  the 
bead  seat  and  contained  air  is  shown  for  four  cases  in  Figures  92 
through  95.  Figures  92  and  93  are  for  carbon  while  Figures  94  and  95 
are  for  beryllium.  The  bead  seat  peak  temperature  is  within  about  8% 
for  three  out  of  four  cases  and  is  16%  off  for  the  fourth  case. 

Figure  92.  The  contained  air  temperature  is  the  node  farthest  from  the 
work  input  and  this  fact  helps  explain  why  in  three  out  of  the  four 
cases  considered  the  predicted  peak  temperature  is  between  14%  and 
35%  different  from  experiment.  Another  inaccuracy  in  the  prediction  of 
temperature  within  the  contained  air  is  that  the  interaction  with  the 
portion  of  the  tire  adjacent  to  the  flywheel  is  neglected.  The  free 
convection  effect  of  the  hotter  air  rising  within  the  tire  cavity  is 
also  neglected. 

As  a  summary  of  the  comparison  between  experiment  and  the  2-D 
theory  peak  temperatures.  Table  10  gives  the  percent  variation  of  the 
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2-D  model  from  the  experimental  peak  values.  The  six  locations  at 
which  temperature  was  predicted  are  listed  in  order  of  expected  accu¬ 
racy  beginning  with  the  most  accurate.  As  can  be  seen  in  Table  10, 
three  out  of  four  variations  in  the  first  five  lines  are  less  than  10%. 
Other  than  the  contained  air  variation,  only  two  of  the  20  variation 
entries  are  substantially  greater  than  10%. 

TABLE  10 

EXPERIMENT  VS  THEORY  PEAK  TEMPERATURE  COMPARISON* 


Location 

Stop 

1.1-3 

Stop 

4.2-3 

Stop 

7.2-3 

Stop 

8.1-2 

Stator  Axial  Centerline 

4.9 

.9 

-2.2 

7.9 

Torque  Tube 

-  5.9 

4.5 

-10.5 

9.6 

Fuse  Plug 

1.3 

.4 

-1.2 

-6.8 

Drive  Lug 

-.8 

-1.6 

28.9 

1.6 

Bead  Seat 

16.2 

6.8 

-8.4 

6.7 

Contained  Air 

-5.4 

14.2 

14.8 

-35. 

*Table  values  are  percent  variation  of  the  2-D  model  from  experimental 
peak  values. 


Two  areas  for  which  plots  are  not  shown  are  in  the  hydraulic  fluid 
location  and  in  the  axial  direction  from  disc  to  disc.  The  theoretical 
model  has  the  capability  to  predict  the  brake  housing  temperature  which 
is  not  dramatically  different  from  the  fluid  temperature  but  the 
capability  to  predict  specifically  the  fluid  temperature  does  not 
exist.  This  capability  would  require  a  more  sophisticated  model  in  the 
housing/fluid  area  and  could  require  a  specific  model  for  each  brake 
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simulated.  On  the  other  hand,  the  two-dimensional  complete  brake  model 
does  contain  the  capability  of  predicting  the  axial  temperature  varia¬ 
tion  from  disc  to  disc.  Although  the  capability  exists,  it  was  noted 
that  the  variation  from  stator  to  stator  or  from  rotor  to  rotor  was  not 
predicted  by  the  model.  This  phenomenon  is  not  considered  very  impor¬ 
tant  for  this  four  rotor  brake,  since,  of  the  ten  temperature  versus 
axial  position  plots  in  Appendix  G,  only  one  shows  a  definite  axial 
temperature  gradient  and  that  is  less  than  a  20%  gradient  through  seven 
brake  discs. 


Conclusions 

1.  The  one-dimensional  axial  model  predicts  temperatures  up  to 
31%  higher  than  experiment  and  no  radial  temperature  gradient  can  be 
determined.  No  cooling  exists  and  thus  no  cooldown  curve  could  be 
predicted. 

2.  The  one-dimensional  radial  model  predicts  temperatures  up  to 
16%  higher  than  experiment  and  no  axial  temperature  gradient  can  be 
determined.  The  cooling  curve  is  close  to  experiment. 

3.  Contact  resistance  at  the  rotor/stator  interface  produces  a 
minimal  effect  in  surface  temperature  difference,  rotor  to  stator,  even 
when  10  times  the  measured  contact  resistance  is  used. 

4.  The  two-dimensional  model  predicted  disc  peak  temperatures 
within  less  than  5%  of  experiment  in  three  out  of  four  cases,  with  the 
fourth  case  less  than  8%. 
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5.  For  brake  discs  with  wear  pads,  it  is  important  to  explicitly 
include  the  wear  pad  in  the  thermal  analysis  and  to  include  the  contact 
resistance  between  the  wear  pad  and  the  disc  heat  sink  material. 

6.  In  a  thermal  analysis  it  is  necessary  to  use  the  variable 
torque  versus  time  curve  in  lieu  of  an  average  torque  curve.  Failure 
to  do  this  can  result  in  deviations  in  the  peak  temperature  of  about  7% 
in  carbon  and  about  23%  in  beryllium. 

7.  The  three-dimensional  program  shows  a  circumferential  tempera¬ 
ture  gradient  of  about  7%  for  carbon  and  about  14%  for  beryllium.  It 

is  not  recommended  that  each  disc  <.  7  the  entire  brake  be  modeled  in 
three  directions  since  two-dimensional  disc  analysis  is  close  to 
experiment  and  since  computer  time  and  storage  dramatically  increase 
with  addition  of  the  third  dimension.  For  the  same  rotor/stator  pair 
analysis  section,  48  nodes  are  required  for  the  two-dimensional  analysis 
and  228  nodes  are  required  for  the  thre  —dimensional  analysis.  If  the 
exact  temperatures  and  gradients  are  required  at  keyways,  the  use  of 
the  three-dimensional  model  for  one  rot^r/stator  becomes  necessary. 

8.  Two-dimensional  analysis  of  a  rotor/stator  pair  gives  good 
agreement  with  experiment  as  far  as  disc  temperature,  including  cool¬ 
down,  is  concerned.  A  variable  convective/radiative  heat  transfer 
coefficient  is  used  in  this  case  and  compensates  for  the  lack  of  heat 


conduction  from  the  discs. 


9. 


A  two-dimensional  analysis  of  the  stators,  rotors,  torque 


tube,  fuse  plug,  drive  lug,  bead  seat  and  contained  air  gives  good 
(within  10%)  correlation  with  experiment  in  most  cases.  A  constant 
convective/radiative  heat  transfer  coefficient  is  used  in  this  case  and 
doc  not  compensate  for  the  conduction  from  the  discs.  This  is  handled 
explicitly. 
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CHAPTER  VII 


UNIQUENESS  AND  OVERALL  CONCLUSIONS 

Uniqueness 

The  following  items  are  unique  to  this  effort: 

1.  The  theoretical  development  includes  six  brake  temperature 
prediction  models  ranging  from  one-dimensional  to  three- 
dimensional  and  the  computer  results  are  compared  between 
models  and  to  experiment.  An  extensive  development  such  as 
this  is  not  reported  in  the  literature. 

2.  The  thermal  analysis  includes  caliper  as  well  as  multi-disc 
developments  and  thus  is  applicable  to  any  vehicle  employing 
disc  type  brakes  (i.e.,  car,  truck,  aircraft,  bus,  etc.). 

3.  Contact  resistance  is  considered  not  only  at  the  sliding 
interface  between  a  stator  and  rotor  but  also  between  the 
wear  pads  and  primary  heat  sink  material. 

A.  The  capability  exists  not  only  to  predict  the  temperature 

within  the  axial  and  radial  directions  of  a  stator  and  rotor 
but  also  to  predict  the  axial  temperature  variations  from 
stator  to  stator  and  rotor  to  rotor. 

5.  The  heat  transfer  coefficient  used  in  the  boundary  conditions 
is  obtained  from  experimental  data  in  two  ways.  The  first 
assumes  only  convective  and  radiative  heat  transfer  as  the 


cooling  mechanisms  (no  conduction) ,  while  the  second 
explicitly  accounts  for  the  heat  conduction  leaving  a 
particular  node. 

6.  The*  thermal  conductivity  is  treated  as  a  variable  with 
temperature  as  well  as  with  direction.  The  equations 
developed,  therefore,  contain  the  derivatives  of  thermal 
conductivity  with  respect  to  temperature  and  direction. 

7.  The  estimate  often  made  of  the  fraction  of  generated  heat 
that  flows  into  a  rotor  or  stator  is  shown  to  be  unnecessary. 
The  particular  fraction  is  a  direct  result  of  the  correct 
application  of  the  first  law  of  thermodynamics  and  is  not  an 
input  into  the  first  law  equation. 

8.  The  importance  of  using  a  variable  brake  torque  input  instead 
of  average  torque  to  the  computer  model  is  clearly  demon¬ 
strated  and  verified  by  comparison  to  experimental  data. 

9.  In  addition  to  the  theoretical  analysis  uniqueness,  extensive 
testing  was  accomplished  in  the  area  of  thermocouple  selec¬ 
tion  and  in  the  area  of  full  scale  multi-disc  brakes. 

Overall  Conclusions 

The  overall  objective  of  this  research  is  to  develop  a  general 
model  for  accurately  predicting  temperatures  within  the  stators, 
rotors,  and  brake  housing  of  a  typical  disc  brake.  This  objective  is 
accomplished  through  the  systematic  development  of  six  brake  temperature 
prediction  computer  models.  The  validity  of  the  models  is  demonstrated 
through  comparison  ta  experimental  data  collected  specifically  for  this 
research. 
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Specific  conclusions  regarding  a  particular  portion  of  the  research 
can  be  found  on  pages  60,  133,  and  237.  These  pertain  to  the  tem¬ 
perature  measurement  analysis,  the  theoretical  development  and  the 
correlation  of  results  chapters,  respectively. 


CHAPTER  VIII 


MODEL  LIMITATIONS  AND  RECOMMENDATIONS 

Model  Limitations 

The  model  selected  from  the  six  developed  in  this  research  as 
being  the  best  design  tool  is  the  two-dimensional  model  with  the 
capability  of  predicting  rotor  and  stator  temperatures  as  well  as  the 
surrounding  brake  hardware  temperatures.  This  model  has  some  limita¬ 
tions  and  those  are  listed  below: 

1.  The  farther  the  distance  from  the  heat  generation,  the  less 
accurate  the  model  becomes.  Thus,  the  predicted  contained 
air  temperature  deviates  from  experiment  in  three  out  of  the 
four  cases  considered  by  15  to  35%. 

2.  The  axle  temperatures  are  not  predicted  by  this  model  in  that 
the  torque  tube/axle  and  rim/axle  interfaces  are  treated  as 
adiabatic. 

3.  The  equations  employed  are  programmed  in  such  a  way  that  the 
task  of  using  experimental  temperature  data  as  inputs  to  the 
program  in  an  attempt  to  calculate  radiation/convection  heat 
transfer  coefficients  cannot  easily  be  accomplished. 

4.  The  brake  housing  containing  the  hydraulic  fluid  is  modeled 

as  a  lump  of  mass  with  the  properties  of  the  metallic  housing. 
Since  this  area  contains  cavities  of  hydraulic  fluid,  retaining 


springs  and  pistons,  the  model  is  not  as  accurate  here  as  in 
other  areas  such  as  the  rotors,  stators,  rim,  and  torque 
tube. 

5.  The  model  is  two-dimensional  and  thus  does  not  give  any 

information  concerning  circumferential  temperature  gradients. 
(This  is  obtainable  from  the  three-dimensional  model,  however.) 

Recommendations 

The  following  items  are  recommended  for  follow-on  efforts: 

1.  If  more  accurate  information  is  desired  in  the  brake  housing 
and  axle,  the  model  should  be  modified  to  contain  more  nodes 
in  these  areas. 

2.  The  present  model  employs  experimentally  determined  heat 
transfer  coefficients.  It  would  be  beneficial  to  conduct  a 
research  program  using  theoretical  techniques  to  predict 
these  coefficients  and  subsequently  to  compare  the  results  to 
experimentally  determined  coefficients. 

3.  In  this  research,  the  finite  difference  numerical  solution 
technique  is  used  along  with  the  Euler  time  marching  tech¬ 
nique.  This  is  based  on  results  from  a  finite  difference, 
finite  element,  and  exact  solution  to  a  one-dimensional 
transient  heat  transfer  problem.  Research  should  be  con¬ 
ducted  to  compare  the  finite  element  technique  to  the  finite 
difference  technique  in  a  transient,  variable-property  heat 
transfer  problem  in  two  and  three  dimensions.  Other  time 
marching  techniques  should  also  be  employed.  Emphasis  should 
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THERMOCOUPLE  CALIBRATION,  CALIBRATION 
CERTIFICATE,  AND  CHECKLIST  FOR 
RUNNING  THERMOCOUPLE  TESTS 
OF  CHAPTER  III 


THERMOCOUPLE  CALIBRATION 


T/C 

Code 

%  Deviation 

At  750°F3 

%  Deviation 

At  1100°F3 

Tests  During 
Which  Used4 

U.1 

-  .173 

+  .054 

24,  27 

U2 

-  .08 

+  .009 

El 

+  .173 

+  .227 

18,  19,  23 

E2 

+  .107 

+  .136 

G1 

0 

-  .036 

Baseline 

G2 

+  .12 

-  .182 

29,  30 

2 

-  .413 

-  .600 

3 

-  .267 

-  .618 

32,  34 

4 

-  .333 

-  .682 

5 

-  .213 

-  .582 

1Manufacturer  tolerance  is  +75%  which  is  ±5.6°F  at  750°F  and 
8. 25°F  at  1100°F. 

2A11  data  using  the  above  thermocouples  were  corrected  using  the 
750°F  t  •  i.ibration  point  for  temperatures  below  925°F  and  using 
the  1100°F  calibration  point  for  temperatures  above  925°F. 

3A  negative/positive  deviation  means  that  the  deviation  must  be 
added/subtracted  to/from  the  thermocouple  reading  to  get  the 
actual  temperature. 

^Thermocouples  not  used  during  testing  were  backup  thermocouples. 


INSTRULAB,  INC. 

1205  Lamar  Street  THERMOCOUPLE 

Dayton,  Ohio  45404 

513-223-2241  CALIBRATION  CERTIFICATE 


Date  11-16-77 
W.O.#  14790 

Cust.  P0 //  04622 


Systems  Research  Labs,  Inc. 

WPAFB,  Bldg.  31,  Area  B 
Dayton,  Ohio  45433 

CALIBRATION  OF  TEN  (JO)  TYPE  "K"  THERMOCOUPLES  (SUPPLIED  BY  SRL) 
750°F  1100°F 


°F 

mV 

°F 

mV 

U1 

-1.3 

(.030) 

+0.6 

(.015) 

U2 

-0.6 

(.014) 

+0.1 

(.002) 

El 

+1.3 

(.032) 

+2.5 

(.060) 

E2 

+0. 8 

(.018) 

+1.5 

(.037) 

G1 

0.0 

(.000) 

-0.4 

(.010) 

G2 

0.9 

(.022) 

-2.0 

(.048) 

2 

-3.1 

(.075) 

-6.6 

(.159) 

3 

-2.0 

(.049) 

-6.8 

(.131) 

4 

-2.5 

(.059) 

-7.5 

(.180) 

5 

-1.6 

(.038) 

-6.4 

(.154) 

AMBIENT 

TEMP: 

80°  ±  5°F 

RELATIVE  HUMIDITY: 

45% 

To  find  correct  temperatures  when  using  above  thermocouples,  subtract 
plus  (+)  deviations  and  add  negative  (-)  deviations  from/ to  readings 
obtained. 

Accuracy  of  results: 

100  to  1600°F  —  ±1.5°F 

Above  1600°F  —  ±2.5°F 

Equipment  Used:  (All  traceable  to  National  Bureau  of  Standards) 

Leeds  &  Northrup  K-3  Potentiometer  Ser.  No.  1526510 

Calibrated  6-2-77 _ 

Traceable  to  NBS  Test  No.  214726,215102 

Eppley  Standard  Cell  Cat.  100  Ser.  No.  841565 _ 

Calibrated  5-27-77 _ 

Traceable  to  NBS  Test  No.  214726,215102 


Plat.  10%  Rhod.  Vs.  Plat,  Working  Standard  T/C 

Calibrated  10-13-77 _ 

Standard  Plat.  10%  Rhod.  Vs.  Plat.  T/C 

Calibrated  10-14-76 _ 

Traceable  to  NBS  Test  No.  215470 

INSTRULAB,  INC. 

(signed  by) :  Robert  A 


.  Dudek 


PROCEDURE  FOR  THERMOCOUPLE  EXPERIMENT 


1.  Mount  carbon  specimen  between  lamps. 

2.  Insert  thermocouples  as  required. 

3.  Turn  on  main  circuit  breaker  to  power  panel. 

A.  Turn  on  secondary  circuit  breaker  to  test  machine. 

5.  Insert  ice  in  thermocouple  reference  junction  container. 

6.  Plug  in  lamp  power  "ready"  circuit. 

7.  Plug  in  and  turn  on  x-y-y'  plotter  and  time  base.  Allow  to  warm  up 
one  half  hour  prior  to  use. 

8.  Zero  both  temperature  channels  on  the  plotter. 

9.  Turn  lamp  control  to  full  power  on  the  manual  setting. 

10.  Move  test  machine  control  lever  to  left  and  wait  for  ready  light  to 
activate. 

11.  Turn  on  compressor  in  machine  shop. 

12.  Turn  on  lamp  cooling  air  making  sure  not  to  exceed  12  psig  down¬ 
stream  of  the  regulator. 

13.  Turn  on  lamp  cooling  water  until  valve  is  fully  open. 

1A.  Energize  lamps  for  about  three  seconds,  then  push  stop  button. 

This  allows  the  carbon  sample  to  heat  to  about  150°F  so  that  a 
constant  starting  temperature  can  be  reached.  Allow  baseline 
thermocouple  temperature  to  reach  100°F. 

15.  Immediately  start  time  base  on  x-y-y'  plotter  and  energize  lamps 
until  the  baseline  thermocouple  reaches  1500°F. 

16.  Depress  stop  button  on  test  machine. 

17.  Allow  plotter  to  record  cooling  curves  until  one  sheet  of  plotter 
paper  is  used,  then  turn  off  time  base  and  move  to  new  sheet. 

18.  Change  test  thermocouples  when  carbon  is  sufficiently  cool  and 
repeat  1A-18  until  desired  thermocouples  are  tested. 

19.  Turn  off  all  equipment,  cooling  air,  cooling  water,  and  the 
compressor  making  sure  to  turn  off  the  cooling  air  prior  to 
turning  off  the  compressor. 
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Parameters 


Figure  E-2.  Test  Parameters  vs  Time  (.Test  2.1-2) 


Figure  E-3.  Test  Parameters  vs  Time  (Test  2.2-3) 


Figure  E-4.  Test  Parameters  vs  Time  (Test  3.1-2) 


■ 

'  f 

. 

L.  .  .1 

--  .  — 

-l 

. 

4 

Figure  E-5.  Test  Parameters  vs  Time  (Test  4.2-3) 
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F-8.  Temperature  vs  Time  for  Stator  2,  Rotor  2,  Torque  Tube  (Test  7.2-3) 
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Stator  2,  Rotor  2,  Torque  Tube  (Test  9.2-2) 
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TEMPERATURE  VERSUS  AXIAL 
POSITION  PLOTS  AT 
VARIOUS  TIMES 
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Figure  G-l.  Temperature  vs  Axial  Position  (Test  1.1-3) 
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Figure  G-5.  Temperature  vs  Axial  Position  (Test  4.2-3) 
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Figure  G-6.  Temperature  vs  Axial  Position  (Test  6.1-3) 
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Figure  G-8.  Temperature  vs  Axial  Position  (Test  7.2-3) 


on/, 


TEMPERATURE  (°F) 


APPENDIX  H 


TEMPERATURE  VERSUS  TIME  PLOTS 
FOR  RADIALLY  LOCATED 
THERMOCOUPLES 
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Figure  H-2.  Temperature  vs  Time  for  Radial  Thermocouples  (Test  2.1-2) 


H-4.  Temperature  vs  Time  for  Radial  Thermocouples 
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Figure  1-1.  Temperature  vs  Radial  Position  (Test  1.1-3) 
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Figure  1-3.  Temperature  vs  Radial  Position  (Test  2.2-3) 
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Figure  1-4.  Temperature  vs  Radial  Position  (Test  3.1-2) 
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Two  methods  were  used  to  calculate  convective/radiative  heat 


transfer  coefficients.  The  first  involved  using  the  classical  approach 
of  assuming  no  heat  conduction  from  the  stators  or  rotors  to  the  rim, 
torque  tube  or  any  other  surrounding  brake  hardware.  The  following 
equation  is  applied  to  the  heat  sink,  all  rotors  and  stators: 

h  A  (T  -  T„)  =  p c  v  §  0-1) 


Using  experimental  cooling  data  in  the  form  of  temperature  versus  time 
plots,  the  temperature  and  time  increment  can  be  found.  This  informa¬ 
tion  along  with  the  material  properties  and  disc  geometry  allow  all 
quantities  in  Equation  (J-l)  to  be  known  except  for  the  h  value.  Using 
the  Euler  time  marching  technique  and  solving  Equation  (J-l)  for  h 
results  in: 


pC  V  (T’-T) 
A^T-T  \ At 


(J-2) 


where  T'  is  the  temperature  of  the  heat  sink  at  a  time  increment.  At, 
from  the  temperature,  T.  Evaluation  of  the  convection/radiation  heat 
transfer  coefficient  in  this  way  assumes  only  convective/radiative 
cooling,  thus  the  resulting  values  were  used  only  in  theoretical  models 
that  did  not  contain  conduction  from  the  heat  sink.  These  models  were 
one-dimensional  radial,  two-dimensional  radial/axial  and  the  cooldown 
portion  of  the  three-dimensional  model.  The  resulting  convective/ 
radiative  heat  transfer  coefficients  plotted  versus  temperature  are 
shown  in  Figure  J-l  for  carbon  and  in  Figure  J-2  for  beryllium.  The  h- 
values  are  different  for  beryllium  and  carbon  due  to  disc  geometry 
differences  and  due  to  the  fact  that  heat  shields  were  used  during 
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Beryllium  Combined  Heat  Transfer  Coefficient 


carbon  tests  but  not  during  beryllium  tests.  On  each  figure  two  curves 
are  shown,  one  for  which  no  exterior  cooling  was  provided  and  one  for 
which  the  cooling  blowers  on  the  test  dynamometer  were  used. 

The  second  method  employed  for  calculating  convective/radiative 
heat  transfer  coefficients  involved  including  the  conduction  heat 
transfer  from  the  disc  to  surrounding  brake  hardware.  Thus  the  first 
law  equation  for  the  heat  sink  would  appear  as: 


’SR  +  ’SI  +  ’SH  +  h  A  (T  "  T»)  "  »  C 


(J-3) 


where  q  terms  are  the  conduction  heat  transfer  rate  terms  and  the 
subscripts  are  defined  as:  SR  -  sink  to  rim,  ST  -  sink  to  torque  tube, 
and  SH  -  sink  to  housing.  The  first  law  was  also  applied  to  the 
housing,  rim,  wheel,  torque  tube,  and  tire  including  the  conduction 
terms  from  adjacent  brake  hardware  (see  Figure  33).  In  this  way  a 
different  value  of  h  could  be  found  not  only  at  different  temperatures 
but  also  for  each  separate  piece  of  brake  hardware.  To  obtain  h  for 
the  heat  sink.  Equation  (J-3)  was  solved  for  h  to  give: 


P  C  V  At  ^SR  ^ST  +  ^SH 


(J-4) 


The  q  terms  were  calculated  by  setting  them  equal  to  a  temperature 
difference  divided  by  a  thermal  resistance.  The  temperature  difference 


was  obtained  from  experimental  data,  and  the  resistance  was  calculated 
based  on  the  thermal  conductivity,  area  through  which  heat  was  trans¬ 
ferred,  and  the  distance  from  the  area  through  which  heat  was  conducted 
to  the  geometrical  center  of  the  heat  sink.  The  results  of  calculating 
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h  in  this  fashion  was  that  each  piece  of  brake  hardware  had  almost  the 
same  coefficient  as  any  other  piece  and  that  this  value  did  not  change 
significantly  with  temperature.  The  radiative/convective  heat  transfer 
coefficients  calculated  in  this  way  are: 


Carbon  with  Blowers, 


h  =  2. 


BTU 


HR-FT-FT-  F 


Carbon  without  Blowers, 


h  =  1. 


BTU 


HR-FT-FT- °F 


Beryllium  with  Blowers, 


h  = 


2.2 


BTU 

HR-FT-FT- °F 


h  = 


Beryllium  without  Blowers 


.4 


BTU 

HR-FT-FT-°F 
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The  following  three  graphs  show  thermal  conductivity,  heat 
capacity,  and  thermal  diffusivity  for  beryllium  and  carbon  as  a  func¬ 
tion  of  temperature.  The  carbon  thermal  conductivity  is  shown  in  both 
the  radial  (with  ply)  direction  and  the  axial  (across  ply)  direction. 
The  beryllium  data  was  obtained  from  Brush-Beryllium  Company  while  the 
carbon  data  was  measured  specifically  for  this  effort  by  Southern 
Research  Institute. 


TEMPERATURE  (°F) 


TEMPERATURE  (°F) 

Figure  K-2.  Thermal  Conductivity  vs  Temperature 


(UH/gld)  AllAISDJJia  1Vl/\IU3HJL 


Thermal  Diffusivity  vs  Temperature 
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PROGRAM  LISTING  FOR  ONE-DIMENSIONAL 
RADIAL  (IDR)  TEMPERATURE 
PREDICTION  MODEL 
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APPENDIX  M 

CONTACT  RESISTANCE  VALUES 
AND  CALCULATION 
PROCEDURE 


The  procedure  used  to  calculate  the  contact  resistance  at  the 
rotor/stator  sliding  interface  was  to  employ  the  empirical  formula 
developed  by  Yu  P.  Shlykov  and  Ye  A.  Ganin. ^  The  formula  developed 
was  compared  to  experimental  data  from  other  authors  as  well  as  to 
experimental  data  on  four  materials  as  tested  by  Shlykov  and  Ganin. 
The  empirical  formula  is: 


1_ 

RE 


.7  P  10** 


a 


B 


(Ml) 


where  RE  is  the  interface  contact  resistance,  K.  is  the  thermal  con- 

A 

ductivity  of  the  air  trapped  between  the  two  surfaces  in  contact,  h  is 

the  roughness  (arithmetic  average)  of  the  contacting  material  surface 

(assuming  identical  materials) ,  P  is  the  pressure  between  the  two 

surfaces,  is  the  thermal  conductivity  of  the  contacting  material, 

o  is  the  yield  strength  of  the  contacting  material  and  .7  is  an 
15 

.  .  ,  .  .  -1.  The  air  thermal  conductivity 

empirical  constant  given  in  meters 

was  taken  at  the  average  interface  temperature  and  the  actual  values 
were  taken  from  the  work  of  J.  E.  Fontenot,  Jr.40  The  graph  given  by 
Fontenot  was  linearized  to  give: 

-6 

Ka  =  (.00045  T  +  .335)  x  10  BTU/ (IN-SEC-°F)  (M2) 

The  thermal  conductivity  of  the  contacting  material,  K^,  was  also 
selected  at  the  average  interface  temperature  and  was  permitted  to 
vary  with  temperature.  The  yield  strength,  was  taken  to  be 
14,000  psi  for  carbon  and  10,000  psi  for  beryllium.  The  value  of  the 
roughness,  h,  was  measured  for  both  the  carbon  and  the  beryllium  using 
a  profilometer .  The  arithmetic  roughness  value  for  the  carbon  was 
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.00137  in.  while  the  value  for  beryllium  was  .00101  in.  The  final 
expression  used  to  calculate  the  contact  resistance  of  the  beryllium 
heat  sink  was: 


RE  = 


513,270,000  x  Ka  +  .2134  x  p  x 


(M3) 


where  is  found  from  Equation  (M2),  ^  is  found  from  an  interpolating 
subroutine  and  P,  an  input  parameter,  is  the  pressure  between  the  two 
surfaces.  Similarly,  the  final  expression  for  calculating  the  contact 
resistance  of  the  carbon  material  was: 


RE  = 


374,000,000  x  Ka  +  .1524  *  P  *  ^ 
Equations  (M3)  and  (M4)  require  the  following  units: 


(M4) 


K.  - 


BTU 


A  (IN-SEC-°F) 

_  BTU 
^1  (FT-HR°F) 


P  -  psig 


_  (HR  FT-FT-°F) 
\  BTU 


i 
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PROGRAM  LISTING  FOR  ONE-DIMENSIONAL 
AXIAL  (IDA)  TEMPERATURE 
PREDICTION  MODEL 


PROGRAM  IDA  ( INPUT , OUTPUT) 
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BOUNDARY  CONDITION  EQUATIONS  FOR  IIDRA 
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Boundary  condition  equations  derived  for  taking  into  account  the 
wear  pads  that  exist  in  some  crakes  are  based  on  Figure  0-2.  The 
temperature  of  the  pad  is  designated  TP(J)  while  the  pad  thickness  is 
designated  PT.  The  resulting  equations  for  this  interface  condition 
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PROGRAM  LISTING  FOR  TWO-DIMENSIONAL 
RADIAL/AXIAL  (R/S  PAIR)  (TIDRA) 
TEMPERATURE  PREDICTION  MODEL 
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THIS  PROGRAM  ASSUMES  THAT  WEAR  PADS  EXIST 


SES  APPLICABLE  EQUATIONS . IF  IPAO  =  0  NO  WEAR  PADS  ARE  ASSUMED 
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CONTINUE 

RETURN 

END 


APPENDIX  Q 

BOUNDARY  CONDITION  EQUATIONS 
FOR  THREE-DIMENSIONAL 
(R/S  PAIR)  MODEL 


BOUNDARY  CONDITIONS  WITHOUT  A  WORK  RATE  TERM 


For  the  boundary  conditions  in  this  section  the  following  13 
equations  are  used  to  allow  the  remaining  equations  to  be  condensed. 


THO  = 


K0  [T(I, J ,K+1)  -  T(I, J,K) ] 
(r&0)2 


(Q— 1) 


THI  = 


-K0  [T(I,J,K)  -  T(I,J,K-1) ] 
(rA0)2 


THIO  = 


K  [T(I, J,K+1)  -  2  T(I,J,K)  +  T(I, J,K-1) ] 


dh 

9T 


(A0): 


[T(I, J,K+1)  -  T(I,J,K-1)1 


(2A0) 


RO  = 


Kr  [T(I, J+1,K)  -  T(I,J,K)] 
(Ar)2 


(Q-2) 


(Q— 3) 


(Q-4) 


RI  = 


-  Kr[T(I,J,K)  -  T(I, J-1,K) ] 

(Ar)2 


(Q-5) 


RIO  = 


Kr[T(I,J+l,K)  -  2  T(I,J,K)  +  T(I, J-1,K) ] 

(Ar)2 


Kr[T(I,J+l,K)  -  T(I, J-1,K) ] 
(2rArj 


3K 
_ r 

3T 


[T(I,J+1,K)  -  T(I ,  J-.1  ,K)  ] 


(2Ar) 


(Q-6) 
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_ _ _ _ _  .  ..  L 


'  .  -  PIP 


xo 


Kx  [T(I+1,  J,K)  -  T (I , J,K) ] 
(Ax)2 


XI 


XIO 


Kx  [T(I, J,K)  -  T(I-1, J,K) ] 
(Ax)2 


Kx  [T(I+1, J,K)  -  2  T(I, J,K)  +  T(I-1, J,K) ] 
(Ax)2 


3K 
_ x 

8T 


[T(I+1, J ,K)  -  T(I-1, J ,K) ] 


(2Ax) 


-  h 


THH  = 


RH  = 


XH  = 


[T(I,J,K)  -  TJ 
(rA6) 

h  [t(I,J,K)  -  tJ 
Ar 

-  h  [t(i,j,k)  -  tJ 


Ax 


T' (2,10,4)  =  T(2,10,4)  +  2  W(THI+RI+THH+RH+XO) 

T' (3,10,4)  =  T(3,10,4)  +  +  THI+RI+THH+Rh) 

T' (4,10,4)  =  T(4,10,4)  +  2  W(THI+THH+RI+RH+XI+XH) 
T' (4,10,2-3)  =  T(4, 10,2-3)  +W  (THIO+RI+RH+XI+XH) 


(Q-7) 


(Q-8) 


(Q-9) 

(Q-10) 

(Q-ll) 

(Q-12) 

(Q-13) 

(Q-14) 

(Q-15) 

(0-16) 

(0-17) 


373 


(Q— 18 ) 


T'  (2, 4-7, 5)  =  T(2,4-7,5)  +  2W  +THT+Xoj 


T' (2,10,2-3)  =  T(2, 10,2-3)  +  2W  (RI+RH  +  +  XO 


(* 


T'  (3, 4-7, 5)  =  T(3,4-7,5)  +  W(XIO+RIO+2THI) 


T' (3,10,2-3)  =  T( 3, 10, 2-3)  +  W(XIO+THICH-2RI+2RH) 


T’ (2,3,5)  =  T(2, 3,5)  +  2W(ROI-RH+XCH-THI) 


T' (2,10,1)  =  T(2, 10,1)  +  2 W ( R I+RH+THO+XO ) 


T' (3,10,1)  =  T(3, 10,1)  +  2W (RI+RH+XICH-THO) 


T'  (4,10,1)  =  T(4, 10,1)  +  2W (XI+XH+RI+RH+THO) 


T' (2,9,4)  =  T(2, 9, 4)  +  2W  +  THI+THH+XOj 


T' (3,9,4)  =  T(3,9,4)  +  W(RIO+XIO+2THI+2THH) 


T'  (4,9,4)  =  T(4, 9, 4)  +  2W  ( +  XI+XH+THI+THHj 


T' (4, 9, 2-3)  =  T(4, 9, 2-3)  +  W(RIO+THIO+2XI+2XH) 


T'  (3,3,5)  =  T(3, 3, 5)  +  2W  +  RO+RH+THI^ 


T'  (4,9,1)  =  T(4, 9, 1)  +  2W  (  +  Xl+XH+THO 


( 


RIO 


T' (2,8,4)  =  T(2,8,4)  +  -W  (2TH I+THO+THH+2RI+RO+RH+3XO) 


(Q-19) 

(Q-20) 

(Q-21) 

(Q-22) 

(Q-23) 

(Q-24) 

(Q-25) 

(Q-26) 

(Q-27) 

(Q-28) 

(Q-29) 

(Q-30) 

(Q-31) 


(Q-32) 


(Q-33) 


T' (3,8,4)  = 


T' (2,8,5)  = 

T' (3,8,5)  = 

T' (6,8,1)  = 


T' (6, 8, 2-4) 


T' (7 ,8,1)  = 


T' (7,8, 2-4) 


T' (6,8,5)  = 

T’ (7,8,5)  = 


T'(6,4-7,5) 


T' (7,4-7, 5) 


T(3,8,4)  +  -  W(l. 5XIO+2THI+THO+THH+2RI+RO+RH) 


T(2,8,5)  +  2W(RI+RH+THI+X0) 


T(3, 8, 5)  +  2W  +  RI+RH+THlj 


T(6, 8, 1)  +  2W  +  RI+RH+THoj 


=  T(6, 8, 2-4)  + 


"(¥ 


2  +  2RI+2RH+THIO J 


T(7,8,l)  +  2W(RI+RH+XI+THI) 


=  T(7,8,2-4)  +  2W  p  +  RI+RH+XI 


T(6,  8, 5)  +  2W  +  RI+RH+THI  j 


T(7, 8, 5)  +  2W(RI+RH+XI+THI) 


=  T(6,4-7,5)  +  W (XIO+RIO+2THI) 


=  T(7,4-7,5)  +  2W  +  XI+THlj 


T' (6,3,5) 

T’ (7,3,5) 
T'  (7,1,4) 

T' (7,2,4) 


T(6, 3, 5)  +  2W  +  RO+RH+THI 

T(7, 3, 5)  +  2W (RO+RH+XI+THI ) 
T(7,l,4)  +  2  W(RO+RH+XI+THI) 

T(7, 2 ,4)  +  2W  +  XI+THlj 


(Q-34) 

(Q-35) 

(Q-36) 

(Q-37) 

(Q-38) 

(Q-39) 

(Q-40) 

(Q— 41) 
(Q— 42) 

(Q— 43) 

(Q— 44 ) 

(Q-45) 

(Q-46) 

(Q-47) 
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T* (7,3,4)  = 

T(7, 3, 4)  +  -  W(RI+3XI+2RO+RH+2THI+THO+THH) 

(Q-48) 

T’  (6,3,4)  = 

T(6, 3, 4)  +  W(l. 5XIO+2THI+THO+THH+RI+RH+2RO) 

(Q-49) 

T1 (6,1,4)  = 

T(6, 1,4)  +  2W  +  THI+THH+RO+RH^ 

(Q-50) 

T’ (7, 1,2-3) 

=  T(7, 1,2-3)  +  2W  ^RO+RH+  +  XI  j 

(Q-51) 

T' (7,1,1)  = 

T(7,l,l)  +  2  W (RO+RH+THO+XI) 

(Q— 52) 

T' (7, 2-7,1) 

=  T(7, 2-7,1)  +  2W  +  THO+Xlj 

(Q— 53 ) 

T' (5,1,4)  = 

T(5,l,4)  +  2  W ( XO+XH+RO+RH+TH I+THH ) 

(Q-54) 

T' (5,2,4)  = 

T(5, 2,4)  +  2W  (— •  +  XO+XH+THI+THH^ 

(Q-55) 

T'  (6,2,4)  = 

T(6, 2,4)  +  2W  ^  +  THI+THH^ 

(Q— 56 ) 

T' (6, 2-7,1) 

=  T(6, 2-7,1)  +  W (XIO+RIO+2THO) 

(Q-57) 

T' (6,1,1)  = 

T(6, 1, 1)  +  2W  y  +  RO+RH+THO^ 

(Q-58) 

T'(5, 1,2-3) 

=  T(5, 1,2-3)  +  2W  +  RO+RH+XO+XH^ 

(Q-59) 

T’ (5,1,1)  = 

T(5, 1,1)  +  2  W (XO+XH+RO+RH+THO) 

(Q-60) 
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r'  (5,2,1)  =  T(5,2,l)  +  2U  +  XO+XH+THoj 

T' (5, 2, 2-3)  =  T(5, 2, 2-3)  +  W(2XO+2XH+THIO+RIO) 

T*  (3,3,1)  =  T(3, 3,1)  +  2W  +  RO+RH+THO^ 

T' (3, 4-9,1)  =  T(3,4-9) ,1)  +  W(XIO+RIO+2THO) 

T' (2,3,1)  =  T(2, 3,1)  +  2W  (THO+RCH-RH+XO) 

T' (2, 4-9,1)  =  T(2, 4-9,1)  +  2W  +  XO+THO^ 

T'  (2, 3, 2-4)  =  T(2, 3, 2-4)  +  2W  +  ro+rh+XO 

T'  (3, 3, 2-4)  =  T(3,3,2-4)  +  W(THIO+XIO+2ROh2RH) 
T' (6, 1,2-3)  =  T(6, 1,2-3)  +  W(THIO+XIO+2RO+2RH) 


(Q-61) 

(Q-62) 

(Q-63) 

(Q-64) 

(Q-65) 

(Q-66) 

(Q— 6  7 ) 

(Q-68) 

(Q-69) 


BOUNDARY  CONDITIONS  CONTAINING  A  WORK  RATE  TERM 
Figure  Q-l  shows  the  way  in  which  the  rotor  and  stator  wear- 
surface  elements  were  identified.  The  rotdr  and  the  stator  each  have 
10  different  kinds  of  elements.  For  instance,  all  rotor  elements 
numbered  6  have  five  heat  flow  paths  (r-direction  in,  r-direction  out, 
0-direction  in,  0-direction  out,  x-direction  in).  The  rotor  elements 
numbered  5  are  different  since  they  have  four  heat  flow  paths  (r- 
direction  in,  r-direction  out,  0-direction  in,  x-direction  in).  As 
discussed  in  Chapter  V,  equations  were  developed  for  the  rotor  inde¬ 
pendently  from  the  stator  and  the  two  added  at  the  applicable  time 
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increment.  The  following  20  equations  are  those  developed  for  the 
different  numbered  elements  as  shown  on  Figure  Q-l.  Equations  (Q-l) 
through  (Q-13)  are  again  utilized  so  that  the  equations  will  be 
shorter.  The  stator  elements  are  identified  with  the  letter  S  while 
the  rotor  elements  are  identified  with  the  letter  R. 


SI;  RI  +  THI  +  RH  +  XO  -  W5  =  ~  TX5?1>K)1 

3  W 


S2 ;  RI  +  THIO  +  RH  +  X0  -  W5=  f T-^5 ?  J > 5)-_ 


S3;  RI  +  THO  +  RH  +  XO  -  W5  =  -(5->  -  T-^-5? ^K— 


S4;  m  +  THO  -t  XO  -  W5  =  .^'C.5^K)  -T^JJC)! 
z  w 


S5;  +  thi  +  xo  -  W5  = 


„r  RIO  THIO  *  FT* (5, J,K)  -  T(5, J,K) ] 

S6;  ~2~  +  ~~2 —  +  xo  -  w5  - - - - 


S7;  RO  +  THI  +  XO  +  RH  -  W5  =  ~  T(j>  J-ijQ-L 


CQ-70) 

(Q-71) 

(Q-72) 

CQ-73) 

(Q-74) 

(Q-75) 

(Q-76) 


S8;  RO  +  THI  +  1.5  XO 


+  ( RH+RI+THH+THO+XH 5 )  _  ^  1.5  [T*(5,J,K)  -  T(5,J,K)] 


W 


S9;  2X0  +  THIO  +  RIO  +  XH  -  W5  =  -2  -tT-  ~  T(5>J*KU 


S10;  RIO  +  2 THO  +  2X0  +  XH5  -  W5 


y  _  2  [T1 (5, J,K)  -  T(5, J.K) ] 
W 


(Q— 77) 

CQ— 78) 

(Q-79) 
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APPENDIX  R 

PROGRAM  LISTING  FOR  THREE-DIMENSIONAL 
(R/S  PAIR)  TEMPERATURE 
PREDICTION  MODEL 
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APPENDIX  S 


EQUATIONS  FOR  TWO-DIMENSIONAL 
(ENTIRE  BRAKE/WHEEL/ 

TIRE)  MODEL 


The  heat  transfer  rates  from  one  node  of  Figure  71  to  another  node  are 
(The  definition  of  all  symbols  is  shown  as  comment  statements  in 
Appendix  T) . 


Q78  = 


(T7-T8)  A7T8 
L7T8  LRT7 
VKR  VKR 


(S-l) 


n?  (T7-T6)  A6T7 
L7T6  L6T7 
VKR  VKT 


(S-2) 


075  =  (T7-T5)  A5T7 
^  L7T5  L5T7 
VKR  VK5 


(S-3) 


065  =  (T6-T5)  A5T6 
^  LTT5  L5TT 
VKT  VK5 


(S-4) 


045  =  (T4-T5)  A5T6 

'  7  rprp  r  T  C  rrm 


LTT5  L5TT 
VKT  VK5 


(S-5) 


034  =  (T3-T4)  A3T4 
L3T4  L4T3 

trim  xnrrn 


(S-6) 


0]n=  (T11-T3)  ARITI 
^  L3T11  (RT+ST) 

VKR  VKR 


(S-7) 


_  (T9-T1)  A1T9 
W  ~  L9T1  L1T9 
VKTT  VKTT 


(S-8) 


Q92  = 


(T9-T2)  A2T9 
L9T2  L2T9 
VKTT  VK2 


(S-9) 
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Q1211  =  -(T12-T11)  ARITI 
/2  LllT3\ 

\  VKR  J 


(S-10) 


Q1312  =  -(T13-T12)  ARITI 
(2  L11T3\ 

\  VKR  ) 


(S-ll) 


Q1413  = 


(T14-T13)  ARTTT 
/ 2L11T3\ 
y  VKR  J 


(S-12) 


0147  =  .(T14-T7)  ARITI 
L11T3  L7T14 

VKR  VKR 


(S-13) 


Q145  =  .(T14-T5)  ARIQ 
LRIO  ,  L5TR 
VKR  +  Vk5 


(S-14) 


Q135  =  _(T13-T5)  ARIQ 


LRIO  L5TR 
VKR  VK5 


(S-15) 


Q124  =  (TL2-T14)  ARIQ 


LRIO  LTTR 
VKR  VKT 


(S-16) 


Q114  =  (TH-T4)  ARIQ 


LRIO  LTTR 
VKR  VKT 


(S-17) 


Q219  =  XT21-T9)  ATITI 


L11T3  L9T21 

VKTT  VKTT 


(S-18) 


Q2221  =  1T22-T21)  ATITI 
(2  L11T3\ 
y  VKTT  / 


(S-19) 


mooo  _  (T23-T22)  ATITI 
^  / 2L11T3\ 

( VKTT  ) 

(S-20) 

no/.  oo  _  (T24-T23) ATITI 
/2L11T3  \ 

\VKTT  ) 

(S-21) 

_  (T23-T8)  A23T8 
(L23T8  L8T23 

VKTT  VKR 

(S-22) 

QH1  =  H  A1  (Tl-TINF) 

(S-23) 

QH2  =  H  A2 (T2-TINF) 

(S-24) 

QH3  =  H  A3 (T3-TINF) 

(S-25) 

QH4  =  H  A4 (T4-TINF) 

(S-26) 

QH5  =  H  A3 (T5-TINF) 

(S-27) 

QH6  =  H  A6  (T6-TINF) 

(S-28) 

QH7  =  H  A7 (T7-TINF) 

(S-29) 

QH8  =  H  A8 (T8-TINF) 

(S-30) 

QHR  =  H  AR(TR(I)-TINF) 

(S-31) 

QHT  =  H  AT (TSS (I)-TINF) 

(S-32) 

OR  (I)  =  *  ~*  t  TRA  TR(I)]  ARTR 

LRITR  LRTRI  (S-33) 

. 5VKR  . 5VKKR 
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QS  (i) 


. 5 [TSA-TSS (I) 1  ATTS 
LTTS  LSTT 
VKTT  +  VKKS 


(S-34) 


The  temperature  at  the  new  time  step  in  terms  of  the  temperature  at  the 
previous  time  stp  is: 


TNI  =  T1  +  (Q91-QH1) 

P  VL* 


(S-35) 


TN2  =  T2  +  ^  (Q92-QH2) 


(S-36) 


TN3  =  T3  +  (Q113-Q34-QH3) 


(S-37) 


TN4  =  T4  +  (Q34+Q114+Q124-Q45-QH4) 


(S-38) 


TN5  =  T3  +  (Q45+Q135+Q65+Q145+Q75-QH5) 


(S-39) 


TN6  =  T6  +  ^  (Q76-Q65-QH6) 


(S-40) 


TN7  =  T7  +  ^  (Q147-Q76-Q78-QH7-Q75) 


(S-41) 


TN8  =  T8  +  ^  (Q78-QH8+Q238) 


(S-42) 


TN9  =  T9  +  ^  [QS(1)-Q92-Q91-QHT(1) ] 


(S-43) 


TN11  =  Til  +  [QR(1)+QR(2)-Q113+Q1211-QHR(1)-QHR(2)-Q114 ]  (S-44) 
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TN13  =  T13  +  —  [QR(5)+QR(6)+Q1413-Q1312-Q135-QHR(5)-QHR(6) ]  (S-46) 
TN14  =  T14  +  -^  [QR(7)+QR(8)-Q1413-Q147-Q145-QHR(7)-QHR(8) ]  (S-47) 
TN21  =  T21  +  [QS (2)+QS (3)+Q2221-Q219-QHT(2)-QHT(3) ]  (S-48) 
TN22  =  T22  +  ~  [QS(4)+QS(5)+Q2322-Q2221-QHT(4)-QHT(5) ]  (S-49) 
TN23  =  T23  +  ~  [QS(6)+QS(5)+Q2423-Q2322-QHT(6)-QHT(7)-Q238]  (S-50) 
TO24  =  T24  +  ~  [QS(8)-Q2423-QHT(8) ]  (S-51) 
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